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Chemistry. — Phase determination in direct Fourier-synthesis of crystal 
structures, By J. M. BIyvoeET. 


(Communicated at the meeting of March 26, 1949,) 


A direct Fourier synthesis of the electrondensity in crystals is based on 
the fact that the amplitude of a density wave is proportional to that of 
the X-ray beam reflected by the corresponding netplane. The phases of 
the reflected waves being lost in ordinary X-ray methods 1), the phase 
relations between the different density components remain unknown at 
this stage. 

The method of isomorphous substitution A,R- AgQR introduces, along 
familiar lines, known reference waves, viz. the density components of the 
A-configuration (for atomic scattering power A,—A,; we suppose this 
configuration to be centrosymmetrical and its determination accomplished). 

The reference waves being known, their superposition effect on the 
density waves of the A,R structure reveals their mutual phase differences, 
ie. the unknown phases ax: of the density waves of the AiR structure 
relative to the centre of symmetry of the A-configuration: 


2 2 z 
Pakts, p = Fakig, p + Paki, 4 1 2 |Paetla, p Paris 4, 608 Phkla, p+ 


Here |F4.r| is derived from the diffraction-intensities of A,R; Fa,p | 
idem for crystal Aj:R; F4,-4, is calculated for the A-configuration of 
atomic scattering power A,—Ay. 

Now this procedure gives us gxxi except for its sign, thus performing 
the phase (sign) determination only for the case of a centrosymmetrical 
A;R structure (projection). Otherwise a synthesis with both max: and 
— nxt can be resorted to, resulting in a duplicated structure-model 2). 

Now we wish to call attention to the fact that in this non-symmetrical 
case there is, in principle, a general way of determining the sign of pax.’ 
We can use the abnormal scattering of an atom for a wavelength just 
beyond its K-absorption limit. This effect is made use of already in X-ray 
analysis to discriminate between atoms of scattering power nearly equal 
under normal conditions 3). Even the abnormal phase shift has been used, 
as long ago as 1930, by COSTER in a manner similar to that proposed 
here 4). Let the abnormal phase shift introduced at atom A be 6. Then the 


1) Cf. those methods in which the scattered beam is made to interfere with the 
drect beam, ia. H. OTT, Ann. Physik 31, 264 (1938), D. GABOR, Nature 161, 777 (1948). 

2) Cf. C, BOKHOVEN, J. C. SCHOONE, J. M. BIJVOET, Proc. Kon. Ned. Akad, v. 
Wetensch., Amsterdam, 52, 120 (1949). 

3) See R. W. JAMES, The Cryst. State II (1948) Chapter IV (1). 

4) D. COSTER, K. S. KNOL and J. A. PRINS, Z. £. Physik, 63, 345 (1930). 
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phase differences + pari between the density waves concerned result in 
the phase differences yaxi + 6 resp. — paki + 0 for the scattered radiation, 
and so will become distinguishable. At present we are testing the applica- 
bility of his method in an actual analysis. 

It will be clear that by the above method it also becomes possible to 
attribute the d or / structure to an optically active compound on actual 
grounds and not merely by basic convention as in organic chemistry. 


Van ‘t Hoff Laboratorium der 
Rijksuniversiteit, Utrecht. 


Mathematics. — Contradictoriteit der elementaire meetkunde. By L, E. J. 
BROUWER. 


(Communicated at the meeting of March 26, 1949.) 


In een vroegere mededeeling 1) is de contradictoriteit der aequivalentie 
der relaties ,,.> 0" en ,,o> 0” aangetoond 2). Hieruit volgt tevens de con- 
tradictoriteit der aequivalentie der relaties ,,= 0" en ,,of = 0 of o>0”, 
Immers laatstgenoemde aequivalentie zou eerstgenoemde impliceeren. 


Beschouwen we de snijpuntsstelling der Euclidische planimetrie, luidende 
dat voor elke twee lijnen a en / van het Euclidische vlak, die néch kunnen 
samenvallen noch evenwijdig kunnen zijn, een gemeenschappelijk punt kan 
worden aangewezen. Nemen we een rechthoekig codrdinatenstelsel aan, 
en kiezen we voor a de X-as en voor / uitsluitend lijnen door het punt 
P(0,1) met richtingscoéfficiént = 0 en = —1, dan leert de snijpuntsstel- 
ling in het bijzonder, dat voor elke lijn 1 door P met richtingscoéfficiént 
<0 en = —I een snijpunt S met de X-as kan worden aangewézen, der- 
halve een natuurlijk getal n(/) kan worden aangegeven zoodanig dat 
x, <2", derhalve voor de richtingscoéfficiént o(1) van 1 de relatie 
e(!) <—2-", dus de relatie o(1) <o 0 geldt. Daar derhalve uit de snij- 
puntsstelling der planimetrie de contradictoor gebleken aequivalentie der 
relaties ,,< 0" en ,,<° 0" volgt, is de snijpuntsstelling der Euclidische plani- 
metrie eveneens contradictoor gebleken. 

De hieruit voortvloeiende contradictoriteit der Euclidische planimetrie 
is evenwel door deze uitsluitend in de verte plaatsvindende manifestatie 
nauwelijks afdoende gedemonstreerd, terwijl de snijpuntsstelling der pro- 
jectieve planimetrie, luidende dat voor elke twee lijnen van het projectieve 
vlak, die niet kunnen samenvallen, een gemeenschappelijk punt kan worden 
aangewezen, hierdoor nog niet wordt aangetast. 


Beschouwen we wederom1) den met het eenheidscontinuum samenval- 
lenden puntwaaier J. Zij f een willekeurig punt van J, verstaan we onder 
az de assertie, die f rationaal verklaart, en beschouwen we de soort o der 
sequenties S(y, ay) en de soort 7 der correspondeerende tweezijdige dem- 
pingsgetallen E(y, ay), waarbij y steeds voorstelt de tweevleugelige aan- 
stuiving met kern 0 en telgetallen (—1)"2-" (n= 1,2,...), doch f vrij 
varieert binnen J, met dien verstande dat het nde element k,(“n) van f steeds 
na Cn(y, af) doch véér cn4; (y, ay) wordt geschapen. 


1) Proc. Kon. Ned. Akad. v. Wetensch., Amsterdam, 52, p. 122 (1949). 

2) Als resultaat der mededeeling is de non-aequivalentie van > en o > geformuleerd. 
De bewijsgang levert echter het verder strekkende resultaat der non-aequivalentie van 
o> Uncen +o 0, 
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Zouden nu te eeniger tijd de relaties , 4.0" en ,,6f <0 of >0” voor 
het continuum aequivalent blijken, dan zou daarmee in het bijzonder voor 
elk element e van 7 kunnen worden aangetoond: 

hetzij e <0, in welk geval de onmogelijkheid van e °>0, dus de absur- 
diteit der rationaliteit van f zou zijn gebleken, . 

hetzij e > 0, in welk geval de onmogelijkheid van e<o 0, dus de niet- 
contradictoriteit der rationaliteit van f zou zijn gebleken. 

Derhalve zou het eenheidscontinuum zijn gesplitst in twee puntkernsoor- 
ten, in elk van welke een element zou kunnen worden aangewezen, het- 
geen onmogelijk is 3). Zoodat van de aequivalentie der relaties ,,54 0" en 
»Of <0 6f >0” voor het continuum de absurditeit is gebleken. 

Hieruit volgt tevens de contradictoriteit der aequivalentie van ,,a = a” 
en ,,of a= 0 of a = 0” voor het continuum. Immers laatstgenoemde aequi- 
valentie zou eerstgenoemde impliceeren. 


Beschouwen we een Euclidisch vlak met een rechthoekig coérdinaten- 
stelsel. Zij a de X-as, m de lijn x = 1, a een willekeurig reéel getal + 0, 
f de absolute waarde van a, P het punt (—1, 28), Q het punt (1,a),/ de 
verbindingslijn der punten P en Q, y een zoowel van + als van 3 verwijderd 
liggend reéel getal. Dan snijdt de verbindingslijn van P met het punt (y, 0) 


de lijn m in het punt (. 287 | dat met Q slechts kan samenvallen, als 


hetzij 2 6 oa == 6, hetzij 2.8 me = —, d.w.z. als hetzij B(y—3) = 0, 
hetzij 6(3y—1)—= 0, welke vergelijkingen beide valsch zijn, zoodat het 
onmogelijk is, dat de verbindingslijn 1 van P en Q de X-as a snijdt in het 
punt (y, 0). 

Mocht derhalve voor de lijnen / en a een snijpunt kunnen worden aan- 
gewezen, dan is dit of het punt (4,0), of het punt (3,0). In het eerste 
geval is a noodzakelijk <0, in het laatste geval > 0. 

Gold dus de snijpuntsstelling der projectieve planimetrie, kon derhalve 
voor iedere a 4 0 een gemeenschappelijk punt van a en I worden aange- 
wezen, dan zouden de relaties ,,=40" en ,,6f <0 6f> 0” voor het conti- 
nuum aequivalent zijn. Zoodat wegens de contradictoriteit dezer aequi- 
valentie eveneens de snijpuntsstelling der projectieve planimetrie contra- 
dictoor is gebleken, waarmee zoowel van de Euclidische als van de 
projectieve planimetrie de contradictoriteit afdoende is gedemonstreerd. 

Verwacht mag worden dat, met behulp van uit vrije variatie van f voort- 
vloeiende soorten van dempingsgetallen door af, ook van andere reeds 


onjuist gebleken theorieén der klassieke wiskunde de contradictoriteit zal 
kunnen worden vastgesteld. 


%) Vgl. Mathem. Annalen 97, p. 66, noot 10), 


Mathematics. — Approximate division of an angle into equal parts. By 
J. G. VAN DER CorpuT and H. Mooy. 


(Communicated at the meeting of March 26, 1949.) 


The editors 1) of Mathematica A published an approximate construction 
of the trisection, due to M. MARTENS, a former head-master. Although 
this construction is very accurate for acute angles (with an error smaller 
than 21’ 24”), it is surpassed in this respect by the very simple construction 
given by S. C. VAN VEEN 2), of which the error for acute angles is less 
than 2’ 36”. 

These articles led H. Moolj to deal in his thesis3) with the problem of 
dividing a given angle approximately into a number of equal parts and 
even of the approximate construction of ra, where a is a given angle and r 
denotes a number between zero and one, which can be constructed with a 
pair of compasses and a ruler. 

Moolj’s construction is based on the following theorem. 


First theorem. 


Describe a circle, the centre of which coincides with the vertex O of the 


given angle AO B= 2a < 180°. 
Suppose that the circle intersects the legs of the angle at A and B. 
Let C be the middle of the segment AB. Produce AB to D in such a 


way, that 


2, 
ep Ac: 
4r 


3—3r? ; boss 
Describe a circle with D as centre and 4 AC as radius, which inter- 
cr 


sects the smaller arc AB at E. Then ~COE is approximately equal to 
ra and the difference is smaller than 
4r(1—r’) (tg¢a+1sin a— a). 
Thus the trisection (r= 4) gives the following construction, which is 


identical to VAN VEEN’s. 

Describe a circle with the unit as radius, the centre O of which is the 
vertex of the given angle AOB= 2a< 180°. Assume that the circle 
intersects the legs of the given angle 2a at A and B. Let C be the middle 


of AB. Produce AB to D such that 
CD=ZAC. 


1) Trisectie, Mathematica A 6 (1937—38), p. 1—4. 

2) §, C. VAN VEEN, Benaderde trisectie, Mathematica A 7 (1938—39), p. 229—237. 

3) H. Mooy, Over de didactiek van de meetkunde benevens benaderingsconstructies 
ter verdeling van een hoek in gelijke delen, Thesis Amsterdam 1948. 
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Describe a circle with D as centre and AB as radius, which intersects 
the smaller arc AB at E. Then “COE is approximately equal to } a and 
the difference is smaller than 


$i (tgza+7sina— Ja). 
In this article the proof of the results found by Moor will again be 
given. Further we shall give a second approximate construction of ra, 


which is a little more complicated, but gives a more accurate approximation. 
The new construction is based on the following theorem. 


Second theorem. 


Describe a circle with the unit as radius, the centre O of which coincides 
with the vertex of the given angle AOB= 2a < 180°, 

Assume the circle intersects the legs of the given angle at A and B; 
let C be the middle of AB. Produce AB to D such that 


1 1 
ness 2 ya en 3 
CD=7,84r)AC + 7407 r?)(9—r?) AC?. . . (1) 
Describe a circle with D as centre and the radius 
3 ae: ny Se ea 3 
act! r’) AC + 549; r).(9-—-£7) AC*, 5. = = aield) 


If this circle intersects the smaller are AB of circle O at E, then 
ZCOBE is approximately equal to ra and the difference is smaller than 


4r(tg4a+4sin a+ 5 sin? a—,3, asin? a—a), 
As for the trisection (r= 4) we get the following approximate con- 
struction. 
Describe a circle with the unit as radius, the centre O of which is the 
vertex of the given angle OA B= 2a< 180°. 
Assume that the circle intersects the legs of the given angle 2a at A and 
B and that C is the centre of AB. Produce AB to D, such that 


CD=fAC +, AC*. 
Describe a circle with D as centre and with 
2AC+ 8, AC? 
as radius. This circle intersects the smaller arc AB at E. Than “COE 
is approximately equal to! a and the difference is smaller than 
$ (tga + 4sin a+ 5% sin'a— 8 a sin? a— $0), 
This construction becomes somewhat simpler, if the factor (9—r2) in the 


second term of (1) en (2) is replaced by 9. 
If r= 4, we get 


CD=tAC+ AG? ee 
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and the radius of the circle with D as centre becomes 


VW ANCR SF es re Ps 


This approximation is extraordinarily accurate. In order to have a 
survey of the accuracy of the approximation of the trisection, we give the 
following tables, referring to the three constructions I, II and III; here I 
denotes VAN VEEN’s construction (which also occurs in Moolj’s thesis); 
II is the construction according to the formulae (1) en (2) and finally III 
is the simpler construction in which the formulae (3) and (4) are applied. 


Difference from 1/3 a, 


180° Lae Pde ee 9 Leo ery 
120° ) ive 19" ee CAL NG; 9 
90° zZ 138" a Fee Seyi 
60° 18” 229. Zao 


For the proof of the above assertions we need the following lemma 


Lemma. 


Assume the real numbers A, B, C and D satisfy the inequalities 
A Ba CO; A?+B?> D*; CHD 
and the real number t satisfies 
A sint+ Bcost=D. 
Then the equation 
max GO cog ces Cee 6! wp a a aK) 
has at least one real solution x for which the inequality 
yA 1 ee 
VA2+ B?—C? + VA? + B?— D? 
holds and for which x—t has the same sign as 


woe 
Acost—Bsint 


Jx—th< 


Moreover 


A cos t—B sin t= + YA? + B?—D? +0. 
For the proof we substitute x = t + u in (5) and we introduce z determined 
by 
tg4u=z (—a<u<a). 
Equation (5) becomes 


A sin (t+ u)+ Bcos(t+u)=C, 
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whence 


I ees 2z L—2 225 as 
a(t eentt py pont) +B (Tg cea a Oo 


consequently 
2?(A sin t + Bcos t+ C)—2z(A cos t—B sin t)—(A sin t+ Bcost—C)= 0, 
hence 
(D+ C)z2+2(Bsint—Acost)z+C—D=0. . . (6) 
The discriminant of this quadratic equation is 
A =(B sin t— A cos #)? —(C?— D’) 
and therefore 
A=A’?+ B?—-C’?=0. 
The roots of (6) are 
tab ye A OSE TB sin t + \A? + B?—C 
G2 u— ie a D 
A? + B?*— D?— A?—B? + C? 
~ (C+ D)(Acost—B sin t= JA?+ B?—C?)’ 
The sign of +) A?2+B?—C? can be taken equal to that of A cos t—B sin t 
and then the sign of tg 4u is the same as that of 


C=) 
A cost—Bsint’ 


According to our convention u lies between —z and a, so that }u has 
‘the same sign as tg 3u, thus 


Wibece pee eee | 
VA?+ B?—C?+|Acost—Bsint| 
CDI 


TfAt + Bia Cha yA Ae De 


where u = x—t. This establishes the proof. 
For the proof of the first theorem we put 


AQ= 1; /COE=x; CD=psina; DE=qsina and “ BDE=y. 
By projecting OED on OC and on AD, we get the relations 


cos x—q sina sin y — cos a 
and 
sin x + q sinacosy = psina, 
By eliminating y we obtain 


psinasin x + cosacosx=1-+ 4(p?—q?—1) sin? a, 
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thus 
A sinx + Bcosx=C, 


where 
A=psina; B=cosa and C=1+4(p?—q?—1)sin?a. . (7) 


We have to determine A and C in such a way, that x is approximately 
equal to ra; therefore we put 


Msinte--Beosra=D,.. .. .. . (8) 
Then, according to the above lemma, 
Yo 
jra—x| < seer) (9) 


VA?+ B?—C a VA?+ B?—D?° 
Expanding cos ra and sin ra in powers of sin a, we get 


Ce a r? (2? in* : 
ents T AL 6] 


(11?) G21 


(10) 
2.2 
sinra—rsina+ Ai 2 sin? a + 2 


3! 
From the relations (7) and (10) it follows that (8) becomes 
r (17—r?) 


D=psina)rsina + OT sin? a+...} + cosa cosra, 


sinfa +... 


Differentiation of (10) gives 


r? (1?—r?) (32—r?) 


sin? a + 4] Taintat...t 


17— 
reosra=reosa} so Mi 


By substituting this result in the above value of D, we get 


a 4 (tet 
D=psina} rsina+ C57") sinta + sink ata 


a eT Yee 
oA stat ey. : 


(11) 
1752 
+ (1—sin? a) 1 + —— sin 
2! 
Putting 
an = {(2h—1)’—@} ans: ap = 1; y(h)=—(2h—1 + 0c) +2prh 
and replacing sin a by s and r?2 by 9, we find 
@, $7” an-2 {(2h—3)—0} 
eth. s 402 
In order to get a good approximation of ra, we choose p and q in such a 
way, that the expansion of C—D begins with sin6 a. This gives 


+ (p’—q’—1) =pr—1 + 3 (1—2’) and 


pre) (AG) 
ga 4! 5 eae 


B22 


hence 


or+3 __3—3r? 
Paget ee 


(13) 


In what follows an upperbound is given for the error, provided that r 
lies between 0 and 1. Substituting in (12) the value of p given by (13), 
we get 
r+3 

2r 


a (hh) = 2p leo hr = § (h—2) (e—1) 


hence 


@ shapes sy al 


p= an (A201): 


h=0 
In C—D the terms with s°; s2 and s4 disappear, hence 
nos Sees thd) ee ee 


h=3 


= (1—e)? (3?—g) U's®, 


where 


eo ark aS 
By virtue of 0S @ <1 the inequality U < Up is evident, consequently 
C—D& s° (1—e)’ (3—@) Up 
where Up, denotes the value of U at the point 0 = 0. Hence 
C—D = (l1—e)? (1—4o)(C—D), . . . . . (14) 


where (C—D), denotes the value which C—D takes for 0 = 0. 
The denominator of the right hand side of (9) contains the terms 


V A2 + B2—C2 and yA? + B2—D?. Here - 


A’ + B’—C? = p? s? + cos? a—{1 + } (p?—q?—1) s?}? 
90s ot (pie gal se 
Substitution of the values of p and q, found in (13), gives 


oil—p}? (lpi? 
2 2 2 2 
A sate! 7p s*— i6 


st, 


hence 


YER HCO = MEe 5/2 ot 5 EP OTS 


By putting t= ra in the last formula of our lemma and substituting (7), 
we obtain 


+ VA? + B°—D?=psinacosra—cosasinra. . . (16) 
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. * . 2 
With the aid of the series (10), by substituting p = oe , and, putting 
r 


Co = 1 and crn = (4h2—@) ca_y (h>1), we get 


VAR pF 3, 8 | Pe ae ais, 


erii—¢) > (2h— —5) ch—2 athe. 
4 h=1 (2h—1)! 


The right hand side is > 0, consequently 


A? a Tp = 3(l—e) Q 
2 2 2 
jA?+ BP- D?= z s(1 igs 


and a fortiori 
AEB KD > 9 6/2 9 (17) 
a i 9 eee 
The relations (9), (14), (15) en (17) give 


2 (l—g)? (1 —4) (C—D), 
3(1— _——— — ——— 
a2! vim tewt + 2s yay Xe 


From Cy = 1 + 352 and Dy = 2asina+ cosa it follows that 


|ra—x|< 


(C—D),=1+ }sin? a—2asina—cosa 


a 
= sina (‘95 + 4 in a—ta). 


Hence 
2(1—e)? (l—$e)s (9344 gata} 
|ra—x| ee ge 
2.3(1—o) s aia oer 
4r 71 — os 
Consequently 


|ra—x| < $r(1—g) (tg 4a +4 sin a— fa). 


This establishes the proof of the first theorem. For the proof of the second 
theorem we put 


2242's (18) 

Then 
C=1+4 (p?—q?—1)s?=14+Ks?+Lst, , . . (19) 
where P, Q, K and L are properly chosen functions of r. Formula (8) gives 


7 ene) ssinra+cosacosra, 


r 
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hence 


where 
apn = (12—o) (32—@) .. . ((2h—1)?—@) and ay) = 1, so an = {(2h—1)?—e} ani, 


consequently 


_ yy an an-1 Pani Q an_2 

D=3s) @h)I ~ @h—2)! + @h—1)I * Qh—3)! 
(il) (2h—3)? —o 
= BP" an 2) OG Oh—3—0) — Roy t 
(2h—3) =e Q 
(Opens kaa: 
This gives 
oe Capes 
D=2 op ? 4) 

where 


p(h) = —{(2h—1) + @ } | (2h—3)?—e} + 2h} (2h—3)—e} P + 
+. 2h (2h—1) (2h—2) Q. 


The functions P,Q, K and L are chosen such that the expansion of C—D 
begins with s8. By virtue of 


2 : 4 
D=1+ 5 (—I-0 + 2P) +7 {—G+e)ll—2)+ 4(1—e) P-++4.3.2.Q) + 


ded IO) 'y oop Mien ae taek | pees 
ra 0) (9—0).-+ 6 (9-9) P.-+- 655.40} +... 


we get in connection with (19) 


2K=—12p+2Po eee 
24L=— (3+ e)(1—e) +4(1—e) P+ 4.3.2.Q; . . (21) 
—(5+0)(9—0) + 6(9—e)P+6.5.4.Q=0. . . (22) 


Substitution in (19) of the values of p and 2 K, found respectively in (18) 
and (20) furnishes 


qo=(P—oF + 2s5(PO—[p)-+ s* QO? es) 
To simplify we put the right hand side of this equation equal to a perfect 
square, by choosing - 
PQ—Le =(P— 0) Q, hence L=Q. 
From (21) follows 


419) Psi eye ae poeees . . (24) 


25 


Formula (22) gives 
120 Q=6 +0)(9—e)—6 Oo). 2 +8, 


hence 


Now (23) becomes 
eq? =(P—e + Qs’)? and rq=+(P—e4+ Qs?) 


The left hand side of the last relation and also both P—o and Qs? are 
positive, so that the plus sign holds good. This gives 


rq= ae 1180+ (9— -@) s?} 
and by (18) 
; 1—o) (9— 
p= pes ie a AS) 


In order to obtain an upper bound of the error, we deduce from our 
lemma that this error is 


|ra—x|< Pe ee 
) A? + B?—C? + yA?+ B?— D? 
here 
Ces Seti 
fins (2h)! 
and 


p(h)=—{2h—1 + @} {(2h—3)’—9} + 
+ 2h}(2h—3)?—o} P + 2h (2h—1) (2h—2) Q. 
Consequently 
yp (h)=— {(2h—1) + Q} {(2h—3)—e} + 
e) (9—e) 


ee 
Ee te tg bo = + 2ni2h— anes Yee, 
which gives after reduction 


y (h)= ae {—102h? + 267 h—180—o (2h? + 3h—20)}, 


c-D=—° pte CE Th Tit? {102 h?—267 h + 180+ @ (2h? + 3h—20)}. 
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In this relation we have for h24 

an—-2 = (1—e) (1 —4 @) baz (Q), 
where bi_2(@) is a polynomial in ge, which has in the interval 0S oe <1 its 
maximum value at the point e — 0, therefore 


= a2 (1-0) (14.0) br-2 (0) = (1-0) (1-4 0) an-2 (0). 
Further 
102 h?—267 h + 180+ 0 (2h? + 3h—20) S (1 + 3; e) (102 h?—267 h + 180); 
for this relation is equivalent to the inequality 
62 h? + 93 h—620 = 102 h? — 267 h + 180, 
which is evident in virtue of h2—9h + 2020. Hence 


CaD= 
<= (oe) 4 9) & (h—3) © ano(0) (1 + pro) (10242 —267h + 180) 
bm 60 as (2h)! ce a 


This furnishes 
C—D3(1—e)(1—4Fo)(l+yra)\(C—D)o. . - - (26) 
From (7), (8) and (25) we deduce 


LP ce Oe ee =p} (99) = st. 
‘Dz re: sin asin ra + 74) sin asinra-+cosacosra, 
hence 
D)=asina + 3; asin? a+ cosa. 
Further 
C,=1+4s? +; a 
hence 


(C—D) =1-+1s7?+ 5 st—fas—,5as?—cosa, 
consequently 
(C—D), = sin a(tg}a+4sin a + 3; sin? a— 3, asin?a—2a). (27) 
Further : 
A? + B?—C? = p’ sin? a + cos? a—(1 + K sin?a+ L sin‘* a)? 
PQ 
=(F ae na 
where Q = I, hence 
2 
A+ BOOa a +2F 0, +S s+ 
+1—s?—1—K? si—Q? it ee eee 


2PQ _ 
g 


2 
so} 12K +3 k2-2) +51(S—2aK) —arstt. 


327 


Substitution of the values of P, Q and K gives 


yay eg ee le 9 , 9-11 
Pe 0g Se gs Pn = 23 
VA?+ r sy ; + 100 5 -+ 


(9—e) (9—1210) , O-o} « 


me eescones. aeons 
In order to deduce from this relation the inequality 
nee ee 
C= yf tal tae ve eae 126) 


we remark that this inequality is equivalent to 


9I—lle (9—@)(9—121e) ,  (9—e)? 


2 A Be et F oe, 
= 4 ue sie. 36000 S 3600.2 
a eee sa 4 
Sy tetire) (It get spe’ 


and therefore equivalent to 


Lo A pt Ee a (os 
( f Jers 10 * 200 ** 3600 * 


121 get Sogo ules «ay Seal | 
3 ee Bi 3500 a var) ° see : cone 9.961 & =° 
To prove this inequality it is sufficient to consider the most unfavourable 


case, viz. s = 1; in this case the inequality becomes 


369 1101, 163179 18561 ina eel 
400 0 12400 ' 961.3600 © — 961.3600° 9.9612 = 


which is true for every value of 9 between 0 and 1. 
By substituting in (16) the value of p found in (25) we get 


ae aay es Ure (9-11. 9){1—o) r{l—p) 
2 ; as ee rf ee 3. 
ee ag a a 
im s r(2?—o)...((2n—6)?—0) (1—0) (2n—7) { 102n?-369n +339 +(n—3)(2n+7)o } s?”—! 
n=4 120 (2n—1)! ; 
This gives 
aa a), Vo) l—e) | 
\A+B—D= 4r is 240r 
s (29) 
— 1-— 
ee Ida 0) (1+ wre). 
ee Q 
In fact this relation becomes after reduction 
dy 3? + Pfs 0 — dats 08? + ety ? —Tees* =0, 


Pups 
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that is 
dy 8? +0 (SPs — reo S°— Tes sei ee az 0 
and this inequality holds because for 0s <1 the factor of @ is at least 
equal to 
drs — tes — 7 > 0. 
Now by virtue of (26), (27), (28) and (29) the relation (9) becomes 
|ra—x|< 
< 2(1—o)(1—4e)(1+ 31.) sina (tgs Bog esas sin? a— 3 asin? a— #a) 
pea Deiat AV 


which gives after reduction 


|ra—x|<#r(tg}a+}sina+ ¥ sin? a—, asin?a— fa). 
This becomes in the special case r= 4 
°  |ra—x|<$(tg4a+}sina + 3; sin? a—,% asin? a—'Z a). 


This establishes the proof of the second theorem. 


Mathematics. — Een iteratieproces voor de oplossing van een rand- 
waardeprobleem bij een lineaire partiéle differentiaalvergelijking van 


de tweede orde. I. By J. J. DRonKERs. (Communicated by Prof. 
W. VAN DER Woube.) 


(Communicated at the meeting of February 26, 1949.) 


§ 1, 
In dit artikel wordt een iteratieproces behandeld, ter verkrijging van een 
oplossing z(x, y) van een lineaire partiéle differentiaalvergelijking van de 
tweede orde: 


07z 07z 07z Oz dz te 
Bet Page arog font lag tt f= 0 oF es (1) 


waarbij de randwaarden gegeven zijn voor x = 0: 


0 
ia CRNA re Ae ass ah ee ace to] 


Hierbij zijn a, b, c, d, e en f functies van x en y, waarbij in het ver- 
volg die waarden van x en y beschouwd worden, die gelegen zijn in de 
intervallen 0O=x=AenO0Xy=B. 

Nu wordt het volgende verondersteld: 

Van de functies b en c bestaan haar afgeleiden naar x van de eerste 
orde. Deze zijn tevens continu in x, hetgeen ook het geval is met de overige 
functies a, d, e en f. : 

In het bijzonder zijn alle functies a t.m. f, de genoemde afgeleide functies 


Ox 


Oz 
van b en c naar x van de eerste orde, benevens de randwaarden & 
x=0 


en z(0,y) analytisch in y en begrensd. 
Voor de afgeleide naar y van de k* orde van ieder der genoemde functies 
geldt dan analoog als voor de functie a: 


Se 6 a ga ae ee i a a 


waarbij M, en @ onafhankelijk zijn van x, y en k, terwijl M, zo gekozen 
kan worden, dat ze ook voor alle andere functies toegepast kan worden. 
Voor alle functies geldt dus, dat de afgeleiden naar y van de ke orde 
na deling door k! en o* begrensd zijn. 
De iteratiemethode, die in dit artikel wordt behandeld, is verschillend 
van die, welke PICARD1) toepast om een integraal te verkrijgen van de 
hyperbolische partiéle differentiaalvergelijking van de tweede orde, indien 


1) EE. PICARD: Mémoire sur la théorie des équations aux dérivées partielles et la 
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deze — eventueel na een codrdinatentransformatie — van de vorm is: 


07z 
ty es +45 BP 


Een differentiaalvergelijking van deze vorm wordt zelfs in het vervolg 
uitgezonderd. 

Ook is van PICARD een iteratiemethode afkomstig om een oplossing te 
verkrijgen van de elliptische partiéle differentiaalvergelijking van de 
tweede orde. 

In dit artikel wordt de partiéle differentiaalvergelijking van de tweede 
orde (1), die van het hyperbolische, elliptische of parabolische type kan 
zijn, vervangen door 2 lineaire partiéle differentiaalvergelijkingen van de 


fade pete, 


eerste orde. Daarna wordt een iteratiemethode aangegeven om een oplos- 
sing te verkrijgen van deze twee vergelijkingen met randwaarden, die af- 
hankelijk zijn van de randwaarden, welke in (2) zijn aangegeven. Dan 
wordt tevens de integraal z(x,y) verkregen, die aan (1) en de gestelde 
randwaarden voldoet. 

Alle iteratiemethoden betreffende partiéle differentiaalvergelijkingen 
vertonen veel overeenkomst. De in dit artikel weergegeven oplossing 
onderscheidt zich door de voorwaarden, die aan de verschillende functies 
van (1) en (2) zijn opgelegd in verband met het convergentiebewijs van 
de iteratiemethode. Dit bewijs is van een andere structuur, dan gewoonlijk 
bij overeenkomstige iteratiemethoden voorkomt. 

Om de convergentie van de iteratiemethoden van PICARD te bewijzen 
moet o.a. verondersteld worden, dat het gebied waarvoor de oplossing 
geldt, zo klein is, dat hiervoor de te bepalen functie | z(x, y)|, benevens 
Oz 
dx! ly 

Een dergelijke veronderstelling wordt in dit artikel niet gemaakt, zoals 
uit de hiervoren aangegeven voorwaarden blijkt. Er worden alleen voor- 
waarden opgelegd aan de bekende functies a, b, enz. van (1) en aan de 
randwaarden. 

Betreffende de variabele x wordt alleen continuiteit verondersteld, ter- 
wijl de functies in y analytisch aangenomen worden. 

Bij vele in de literatuur der toegepaste wiskunde voorkomende verhan- 
delingen betreffende de convergentie van iteratiemethoden, is het 0.a. nodig 
te veronderstellen, dat alle afgeleiden naar x of y van de bekende functies 
bestaan en tevens begrensd zijn. In verband hiermede moet vaak voor x 
of y of beiden een bepaalde lengte-eenheid gekozen worden (bv. sin 2x 
voldoet niet aan deze voorwaarden, sin x wel). 


en 


kleiner zijn dan een constante c. 


méthode des approximations successives. Journal de mathématiques pures et appliquées. 
Paris 1890. 


E, GOURSAT: Cours d’analyse mathématique, troisiéme édition, tome III, pages 133 et 
229 etc. Paris 1923. 


COURANT und HILBERT: Mathemathische Physik II, pag. 317 ev. 
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Doordat nu echter aangenomen wordt, dat de bekende functies analytisch 
in y moeten zijn, worden de mogelijkheden aanzienlijk uitgebreid, zo kan 


nu bijv. ook de functie Ti beschouwd worden. Eveneens is het niet 


meer nodig om een bepaalde lengte-eenheid te kiezen. 
Betreffende de gemaakte veronderstellingen worden nog enkele opmer- 
kingen gemaakt: 


2 
le. In (1) is de coéfficiént van a3 gelijk aan één genomen. Indien 
x 
2 
deze term ontbreekt en dit niet het geval is met a, wordt de coéfficiént 
y 
hiervan gelijk aan één gesteld. In de randgegevens (2) moeten dan x en y 
verwisseld worden, terwijl ook, hiermede corresponderend, hetgeen van de 
functies a t.m. f gegeven is, gewijzigd zal moeten worden. 


Wij zullen echter wel veronderstellen, dat zeker één der coéfficiénten 

2 2 
van See en ea gelijk aan én gemaakt kan worden. Zoals reeds gezegd, 
beschouwde PIcarD speciaal de gevallen, waarbij of de beide coéfficiénten 
nul zijn, 6f aan elkaar gelijk en dus gelijk aan één gemaakt kunnen worden. 


2e. In het voorgaande is als randwaarde z(0, y) aangenomen. Laat een 
dergelijke kromme algemeen voorgesteld worden door: 


x=9(u); y=y2(u)i z=—q3(u).... . (3) 
terwijl verder 92 @4(u) en = v5(u) bekende functieaszifa, mits het 
verband bestaat: 

Ps (u) YY (4) + @s (u) pf) (a) = oY) (u). 


Wordt in plaats van x, de codrdinaat x, ingevoerd, zodat x = y,(u) + x,, 
dan is de kromme (3) weer in het vlak x, = 0 gelegen. 


3e. Zoals bekend volgen de karakeristicke krommen van de partiéle 
differentiaalvergelijking (1), die in het xy vlak gelegen zijn, uit de ver- 
gelijking: 
dy’?’—bdxdy+adx*?=0. 


Wordt nu een kromme volgens (3) gegeven, dan mogen de functies g, 
en $2 voor geen enkele waarde van u voldoen aan de betrekking: 


2 
pl? — bol!) of + agi? =0 
dat wil zeggen de raaklijn in een punt P van de projectie van de kromme 
(3) op het xy vlak mag niet samenvallen met de raaklijn in P aan een 
karakteristieke kromme, die door het punt P gaat. 
2 
Daar de coéfficiént va:a a3 in (1) ongelijk aan nul is, zal de projectie 
x 


van de kromme z(0,y) op het xy vlak nimmer samenvallen met een 
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karakteristieke kromme of een deel daarvan. Wel kan in één of meer 
punten P de raaklijn aan de projectie samenvallen met de raaklijn aan de 
karakteristieke kromme in P. 

Zoals bekend, worden in verband met het al of niet reéel zijn van de 
karakteristicke krommen, de partiéle differentiaalvergelijkingen van de 
tweede orde onderscheiden in hyperbolische, parabolische en elliptische 
differentiaalvergelijkingen. Alleen in het laatst genoemde geval zijn de 
karakteristieke krommen imaginair. Dan zal dus een gesloten kromme in 
het xy vlak in geen enkel punt aangeraakt kunnen worden door een karak- 
teristicke kromme, hetgeen wel het geval kan zijn bij de hyperbolische en 
parabolische differentiaalvergelijkingen. Deze raakpunten zullen dan uit- 
gezonderd moeten worden. 

In het vervolg wordt niet verondersteld dat de differentiaalvergelijking 
noodzakelijk elliptisch moet zijn. Voor dit geval is de in dit artikel behan- 
delde iteratiemethode uit praktisch oogpunt wel het meest aangewezen. 

4e. In vele praktische vraagstukken betreffende de elliptische partiéle 
differentiaalvergelijkingen is de kromme x — 0, waarop z(0, y) gegeven is, 
gesloten. 

Laten wij weer aannemen, dat binnen het gebied G door x = 0 omsloten, 
de functies a, b, c, d, e en f voldoen aan de eisen in het voorgaande ge- 
noemd. 

Als de integraal z(x, y) op G bepaald is door haar waarden voor x = 0 
en de eis dat zij en haar partiéle afgeleiden tot de tweede orde binnen G 


doorlopend moeten zijn, kan de functie = niet willekeurig gegeven 
x=0 


zijn, maar moet ze zo bepaald worden, dat een dergelijke integraal z(x, y) 
kan worden verkregen, 


$2. 


In plaats van (1) beschouwen wij de twee partiéle differentiaalverge- 
lijkingen van de eerste orde: 


r) é é 
(a) ge tay thay tea 


Ook nu zijn a, 1, 71, 6, By en yg functies van x en y. 
Elimineren wij u uit (4), dan wordt gevonden: 


é¢2 16 é p, 6 
jaa hig ess + a sat (n+ ese Pit ay,— = 


adx a Ox 


0x adx/)’ 
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Nu zullen (1) en (5) identiek zijn, indien de zes functies a, By, 71, 6, Bo 
en 72 voldoen aan de zes vergelijkingen: 


B,=—b; af,=—a; 7, +-—=-c; 

aox 
dp, P, da 6B, Oy, 7 Oa O72 
Rete ae rae 3g 
88d da 
dx adx 


Wij bepalen hiervan een particulier stelsel oplossingen. 

De functie #; is onmiddellijk bekend, terwijl volgens de tweede en derde 
betrekking y; en #2 in a zijn uit te drukken. Elimineren wij daarna y, Bo 
en 7g uit de vierde en vijfde vergelijking, dan wordt een ingewikkelde 
partiéle differentiaalvergelijking gevonden voor a, namelijk: 


2q da 62a 2(éda\? 2a(da\* 2bédada 
pa tags thee 2 (22) ay 


a 
(7) 
6b\ da/ , db 26a id de) ab Sa 


Beye ey) aoa 


Volgens de gegevens (zie § 1) bestaan alle coéfficiénten voor Ox = A; 
0 = y = B en zijn in ieder geval in x continu. Ten opzichte van y zijn alle 
coéfficiénten analytisch. 

Deze partiéle differentiaalvergelijking is ingewikkelder dan de partiéle 
differentiaalvergelijking (1). Nu is het echter alleen nodig om een parti- 


culiere oplossing voor a te bepalen, zodanig dat a(x,y) en x in y 


analytisch zijn en ten opzichte van x continu. Dan behoeft niet aan be- 
paalde randwaarden voldaan te zijn. Daarna kunnen de overige functies 
By, v1 enz. direct bepaald worden en ook hun partiéle afgeleiden naar y 
van de k4 orde. 

Zodra a bepaald is, is 6 volgens de laatste betrekking van (6) onmid- 
dellijk te bepalen als een oplossing van een gewone lineaire differentiaal- 
vergelijking. 

Wij kunnen uit (6), in verband met hetgeen in § 1 is verondersteld, 
a, By, 71, 6, Bz en 72 steeds zo bepalen, dat haar partiéle afgeleiden naar y 
aan de voorwaarden voldoen, dat de afgeleiden naar y van de k* orde na 
deling door k! en een factor o* begrensd zijn en dus ook analytische 
functies in y zijn. 

In het bijzondere geval dat: 


dd , bc 2b da 
ee ef Ya 8 
Re  Sacdy. dy? : ) 


zal a = constant = k een particuliere oplossing zijn. 
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Dan is volgens (6): 


ra 1db 1a, 66 _ 
bj =—); a yt ee ee k dy’ bx 


Ook zal er vereenvoudiging optreden, als de functies a, b, c, d, e en f 
alleen van x of y afhankelijk zijn. 

Het stelsel (6) kan dan vereenvoudigd worden door bijv. te veronder- 
stellen, dat bij afhankelijkheid van x ook de functies a, By, 71, 9, Bo en 72 
slechts van x afhankelijk zijn, zodat 


Bs _ 872 _ 
oy oy 
Dan zal a een particuliere oplossing moeten zijn van de gewone differen- 
tiaalvergelijking: 
2 
a@s—2 ie eee a (ge .] =: 

Ten slotte nog een opmerking over het geval, dat in (1) de functie a 
gelijk is aan nul. Dan kan het mogelijk zijn om (1) door één partiéle 
lineaire differentiaalvergelijking van de eerste orde, nl. van de vorm (4a) 
te vervangen, waarbij dan a=0. Volgens de tweede betrekking van 
(6) is nl. afo== 0, dus of fs == 0 of ¢ =— 0. 

Als Bg = 0 en a 40, is (1) weder door twee lineaire partiéle differen- 
tiaalvergelijkingen van de eerste orde te vervangen. Zodra echter a = 0 
kan zijn, kan (1) worden vereenvoudigd tot één lineaire partiéle diff.verg. 
van de eerste orde, waarin een willekeurige functie van y voorkomt, die 
met behulp van de randwaarden (2) nader kan worden bepaald. 

Dan moet verder nog deze lineaire partiéle differentiaalvergelijking van 
de eerste orde worden opgelost. 

In plaats van de partiéle differentiaalvergelijking (1) en de gegeven 
randwaarden te beschouwen, trachten wij deze te vervangen door het 
stelsel vergelijkingen (4) met daarbij behorende randwaarden, op de wijze 
die in het voorgaande beschreven is. 

Deze randwaarden worden dan als volgt bepaald: 


Daar z(0,y) en fe gegeven zijn, volgt uit (4a): 
Ox } = 


1 Oz 1 1 
ON =p (Sz) _ PHO 20 0,4) -210-2 00,5) —20.9). 0) 


Alle termen van het rechterlid zijn bekend, zodat dit ook het geval is met 
u)(0,y) en dus u(0,y), waarbij aan de integratieconstante een wille- 
keurige getallenwaarde gegeven kan worden. Tevens zijn alle afgeleiden 
naar y van u(0, y) te berekenen. 

Wij bepalen nu eerst de functies z(x,y) en u(x, y) zodanig, dat zij 
aan (4) voldoen, terwijl deze functies voor x = 0 overgaan in de bekende 
functies z(0, y) en u(0, y). 


i i) 


De functie z(x,y) zal dan tevens aan (1) en de gestelde randwaarden 
voldoen. 


$3. 


In deze paragraaf wordt aangegeven hoe de functies z(x,y) en u(x, y), 
die aan (4) en gegeven randwaarden moeten voldoen, volgens een iteratie- 
methode kunnen worden bepaald. 

De randwaarden z(0, y) en u(0, y) geven wij dan kortweg aan met zo 
en Wo en beschouwen deze als gegeven functies. 

Verder wordt onder z\*) de afgeleide naar y van de ke orde verstaan enz. 

Van al deze ke afgeleiden naar y (k= 1, ...... hi Sap roe ) wordt weer 
verondersteld, dat zij na deling door k! en of begrensd zijn en dan kleiner 
zijn dan een bepaald getal M. Het verband met een bepaalde lineaire 
partiéle differentiaalvergelijking van de 2¢ orde, dat in de vorige paragraaf 
is aangenomen, wordt echter voorlopig buiten beschouwing gelaten. 

Als eerste benadering van z(x, y) en u(x, y), Z, en &, genaamd, wordt 
gesteld: 


== x x x 
2, =2 +2,;= 24 ul) fadx+ 2) [ Bidx + 2 [y,dx+ [ddx 
0 0 0 0 
a i (10) 
Z;=u,+u; =u + 2!) [ Bodx + 29 f yd x. 
i) 


Zij geldt voorO =x=A;0=y=B. 

Voor x = 0 is z,; = Zo, zodat z, aan de gegeven randwaarden voldoet. 

In het vervolg zullen wij in de formules de integratiegrenzen 0 en x 
niet meer aangeven. 

Laat in het algemeen: 


Ze-1 29 t--- + 2n—1 CM By uy +a, +... + o,-1 


berekend zijn. Dan worden Zn en un als volgt gevonden: 


8Zn Sin 
Voor x" Ox wordt gesteld: 


Az n—-1 n-1 n-1 
2=a Jo+ 2, SAM+y, Sats. 
Ox k=0 k=0 k=0 


Oiin n-1 n—1 
=p, Fw + 72 Dz. 
Ox k=0 k=0 


Hieruit volgt: 


n—-1 n—-1 n—-1 
25 + fa (22?) dx + /B (2,2) dx oh (2%) dx + [ddx, 


k=0 


B= n—1 
fin = Uy + f[ Bo ( 5 21) dx+ fy. (2+) dx 
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Dus is: 
z= / (aul) + Biz) + 2 + 6) dx, 
z,=f(au + 6,22 +rz)dx, hae eee ie) 
zn =| (e a), + B, 20) , + y,Zn-1) dx 

en: 


uy = { (Bo zy) + 2 Zo) dx 
t= /(hr2 traldxe }...- . + (12) 
ttn = | (B2 2) + 72 Zn-1) dx. | 


In het bijzondere geval, dat zp; en up-1 identiek gelijk aan nul zijn, © 
is dat eveneens het geval met Zn en Un enz. 

In § 6 wordt aangetoond, dat als de functies a, b, enz. voldoen aan de 
voorwaarden in § 1 genoemd, voor de variabelen x en y intervallen be- 
staan, waarvoor de genoemde reeksen gelijkmatig convergeren en z en u 
oplossingen zijn van de differentiaalvergelijkingen (4), die tevens aan de 
gestelde randwaarden voldoen. Worden daarna in de formule voor z, 
de functies ug, u, enz. uitgedrukt in ts ~ en z(0, y) en haar afgeleiden 

=0 : 
naar y, dan wordt de functie gevonden, die aan (1) en de gegeven rand- 
waarden voldoet. 

Elimineren wij u@), met behulp van (12) uit de uitdrukking voor Zn 
(zie (11)), dan wordt gevonden: : 


Zn== fa Wee dx?+fa [B22 jdxt+fa ee Sask 
oie [rPzn-2dx? + [ By 2) ;dx + [71 2nd x. . (13) 


_ In verband met de gegevens is het nl. geoorloofd om de differentiatie 
naar y en de integratie naar x te verwisselen. 
Door tweemaal achtereenvolgens naar x te differentiéren, blijkt: 


Us EO? ase Ts 1 da dzn OZa-~i 


Ox Sa Ue py tae a ae er reeare most es 


gad, de 5 eee 
flee oe (leon) fe 


a 0x oy 
dy, 71 8a OY2 
a be ier ee Zn=1 -P mabe fie 
Volgens het gegeven is ook: 
6? ZnS Tee ra 
dydx Ox dy" 


Wordt ten slotte toch weer aangenomen, dat er tussen de vergelijkingen 
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(1) en (4) het verband bestaat, dat in § 2 is verondersteld, dan kan voor 


(14) geschreven worden: 


6? of fae 67 Z5-5 6? Zn-2 O2n 
ey ae CT re + 


é }é n- ) 6 t= 


é }é 
—(e+0%%) tnt a3? zn. 


Deze formule geldt voor n = 2. 
Voor n = 1 luidt zij als volgt: 


5x2 y} oy 


) 
sta—letn)24—(d+e Yata 2 20 (e442) Zo. 


(15) 


(15a) 


Mathematics. — A non-homogeneous inequality for integers in a special 
cubic field. By A. V. PRASAD. (Second communication.) (Com- 
municated by Prof. J. G. VAN DER CORPUT.) 


(Communicated at the meeting of March 26, 1949). 
Lemma 9. 
ad OO Se ae ee 


Proof. Suppose, if possible, that a = 04— 6-5. Now by (16), (19) 
and (12) 
a<yfa|p|267} <A V5 +8) (I-9} <3 < 
so 
la—O*|< 6. 
Hence, using (18) with é = @-4, 


| (8 — o*) (B—$4)| > (1—2) 4, 


|B—*| > Ay (1—2) > 2-019, 
if ¢ is sufficiently small. Thus 
| B| = | p—o*| — Ae > 2-019— 4-2 > 2:019—0-570 = 1-449. (49) 
But, by (16) and (12) and our supposition about a, 


vey i" our yf oa an 


> 001 5-001 
<V bree tocal < VER <0 


contrary to (49). This contradiction proves the lemma. 


so that 


6. The object of the next eight lemmas is to prove that a must be close 
to 02 + 1 and that 6 must be close to d? + 1. 
For this purpose it is convenient to work in polar coordinates. We write 
Pare! Ge Se woe © a ee) 
Peadew, 1.0. we ee cee 
where r>0, d>0 and y and y are real. 


We suppose that y is chosen so that 0 = y <2a, and that y is chosen 
so that 


SL AO ees i ee ee 
Then, as9 OO =1 andd0+6+ =O, we have 


r= Olas e=t 


and 


= —16'=—0-7623512..., 
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so that 

Be PO Se ei. (83) 
We note that 

2y == 279° 20’-6..., 

21 Toma 9 a OS Sa 

4y = 558° 41’-2..., a) 

5 y = 698° 21’-6 
We also write 

ee are ae TG. 8 (55) 


where R>0 and OS w< 2a, Taking x = 1, y=0,z=1 in (31) and 
(33) we have 


€+67(1+9)(1+¢07=5......, (56) 


Hence 


R=y{(l+ (1 +39} =} gif =1- 3472054... 


and 
—~ 2794+ _ 244+ 9+4+0)—6 
cos © = 2R = 2R 
4—6? 
= gp 0° 833251... 
Thus 
eee Don, uel, 5 gone = (57h 


Lemma 10. For n = 2,3,4 or 5, 
ee 1 l—e 


a as (hkl eae 
cos (x ny)< 2 ee é a Qn |a—6"| 
— 5 or0| (22)! 4) 

Ye, 


Proof. For any rational integer n, we have 
(6) (8-8) = (delz—r" einv) (de-tx—p" e~ inv) 
= d? + r"—2r"dcos(y—ny), 


so that 
cos (y—n py) = {d? + r"—(B—$") (B—$")}/2 r"d. 
Using (18) with = 0-" we obtain, for all rational integers n, 


l—e 


ae weiss 2 Ot Ni Oe te eld es 
cos (7—ny) S 2rd }¢ <2 |a—O" | 


340 


Since, by Lemma 9, 
a< 64—05=20?+0-2< 27 
we have 
\ |a—6"| < an 
for n — 2,3,4o0r5. Thus if ¢ is sufficiently small 
1 l—e 
6" |a—6"| = 


But by (51) and (16) 


for n==2, 3,4; or 5. «1. = {60} 


5+6 


a 


d= pps 
Now (58) follows from this result together with (59) and (60). 
Lemma 11. Jf a= 062+ 1, then 


3220 < p< BIS, ee 
Proof. By Lemma 8, 
a> 0 +6 = 2 01 ee 
Thus 
5+6 _5-+0001 
a ee 
oe a rae 


and by Lemma 10 and (60) 


Bes 1 Ime 
5-0001/ (267-1 ' 6m |a—on|\’ 


for n = 2,3 or 4. Now numerically 


1 ee roe ee. 


cos (y—ny) < ant 


% 
8 . 4 3 
2\5-0001/ ~2\ 5-0001 <4 (0- 50194)? < 0: 35424, 
and 


5:0001  5-0001 
267—1 = 2-50975... 


< 1 399205: 


Also we have 
ja—47| <1, 
ja—635|< 2? 4+1—-AF=67, 
lea—6*|< 6*—-247+1=06~. 
Thus, if « is sufficiently small, 


cos (7 —2y) < (1+32472) (035424) {1.99263 -+ 0- 56985 — 1} 
< 0:73323 < cos 42° 50’, 

cos (y—3y) < (1+ 52471) (0-35424) {199263 + 0-43016 —2-32471} 
< 0-05298 < cos 86° 57’, 

cos (x —4y) < (1+75488) (0+ 35424) {199263 + 0-32472—1-75487} 
< 0-34967 < cos 69° 31’. 
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Using the numerical values for 2y, 3 and 4y we see that it is not 
possible for y to satisfy any of the inequalities 


|x—279° 20’-6|< 42° 49’-9, 
lx—59° 0’- 9| < 86° 56’-9, 
|z—419° 0’-9|< 86° 56’.9, 
|z—198° 41’- 2|<69° 30’-9, 
Hence, by (52), 
322° << x < 333°, 
. Lemma 12. Jf a = 62 +1, then 
Boeck Oar te Ae. a 68) 
Proof. If a= 62+ 1 we have 
e<ott <a. 


Thus by Lemma 10 and (60), 


FEN | }5 000 nS ea 
a [2 oe 
cos (y—ny) << 6” eae, G24] gn la= 


for n = 3, 4 or 5. Now, numerically 


at = ae 75487 . 


5 - 0001 5- 0001 a <4{0+ 55097}! < 0- 37114, 


1) 6241 1$2- 75487. : 
. at : - 37113, 
eee 225-0001 $ >#10 Pome hh OC ILI? 


and 
5-0001  5-0001 
67+] a 1 75487... 


Also, since 62+ 1=a< 64—6-5, 
\a—@| <6 —6>— BF =28?—3 < 0+ 50976 < (1: 96170)" 
la—6*|< 64*-@?—-1= 6-4, 
la—F? |< O—-G—-1=6. 


Binet. 81534. 


Thus, if ¢ is sufficiently small, 


cos (y—3y) < (1-52471) (0-37114) {1-81534 + 0-43016—1-96170} 
<0: 16060 < cos 80°45’, 

cos (y— 4 y) < (1+ 75487) (0- 37113) {1-81534 + 0-32472 —3-07959} 
<—0-61190 < cos 127° 43’, 

cos (y—5 y) < (2-01981) (0: 37114) {1-81534 + 0: 24513 —0-75487} 
< 0+97872 < cos 11° 50’. 
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Using the numerical values for 3y, 4y and 5 y, we see that it is not 
possible y to satisfy any of the inequalities 
la— 59° 0’-9|< 80° 44-9, 
|x—419° 0’-9|= 80° 44’-9, 
|y—198° 41’-2| <= 127° 42°-9, 
|x—338° 21’-6|< 11° 49’-9. 
Hence, by (52) 
326° <4 <327 


Lemma 13. For n = 2, 3,4 or 5, 


1 62 \3( 5 1 1 
cos (r—nv) <> (“5 2 i ane =: =a. 


1(/5+6\3 ee 
tare zn —a}. 
Z a 
Proof. It is clear from (12) that 
5+ 6\3 5 \3 
(Cea ae 
a \t at. 
sf)> (3) 
Using these inequalities together with (60) in (58) we obtain 
1 (/ad"\3( 5 1 1 
cos (y—ny) < 2 ( 5 | a T Ga i |a—6"| 
ere uiew ee 
fol. Poa ene (65) 
fag 3 é 5+ 6\3 
<= Rigs sees: he 
a5) ae oe 
for n= 2, 3, 4 or 5; Now-sitce 
PLO <a< Oo, 
we have 
agr\t(5 4 l—e 
oe Aenean 
ae NS 1 
Saye 6 
<a ee petiele - 6) 
3-1 4rd AS 
<| 5 p02 45-54 


== (2»542)§ (7 > 3) < (1.6) (7-3) < 12. 
Now provided «¢ is sufficiently small (65) and (66) prove (64). 
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Lemma 14. If a= 62 + 1 and n= 2 or 3, or ifa= 62+ 1 andn=4 
or 5, then 


cos (y—nwy) < cos(w—ny) + 1-36 
; 
—4|a—6—1|-4} (344) —a} Ae (67) 


Proof. We suppose that n = 2, 3, 4 or 5 and write 


_ f/x \"45 1 a) 
fa 2) = ( 5 (x? OA” [x—oA]§" (68) 
for x satisfying 
Pee er gs 169) 


Then 
d oh Mea Sa eae: 1 
Eee |=3,(2) +3,(3) apareeat 
x \} 1 
+(2) ame | ve) 


‘aides pred 2x ae De 
=10m (2 Vt 9 | 5 (2-5 + ga) 8) 


We consider two cases. First suppose that n = 2 or 3 and 
yA ged ides wall ited oo» Ge nr 75 2 


Then 
0% < x—O6" S11, 
and by (70) 


d ' / Nas | OY (eck Le Oe 
gp brle)> 2OW! ary X| 208 1) rs A 
1,/ 7 , 
>a) seen} | 22 aime a1 | 


ae 7 emer! Sey, ie ee 
=z) suaoe 2)8~ 


Consequently 


fa (x) S fn (0? + 1)—4|x-—@7—1| . . . . . (72) 
if x satisfies (71). Thus, if n = 2 or 3 and a= 62 +1, it follows from 
(64), (68) and (72) that 


(BP41\8( 5 1 1 
cos (z—ny) < 44" (=; Vig Ber TET 
+ 1-36—4|a—07—1| . (73) 


“4-4 


23 


344 


Now suppose that n = 4 or 5 and 


21s rc = 20 6 2 (74) 
Then 
O05 < O"—x< O—-GP—1=8, 
and, by (70) 
d - } = 1 
qx ln) <4 + 6? (4) ){64— = MA) 
2(67+1 5 1 
«| : a hye (o-9 +1 
—1) 52640) 411 
2(1 340-741 5 —u 
<——4(2-3)(1-4) i288 3 
Consequently 
fn (x) < fa (O? + 1)—4|x—@?—1] . . .°- - (75) 


if x satisfies (74). Thus, if n = 4 or 5 and a= 62 +1, it follows from 
(64), (68) and (75) that (73) is satisfied. 

We now show that (67) follows from (73) in the required cases by 
proving that 


g?+41\3 1 
ak —l1 n|2 ‘ — / 
cos (o—ny)= 4% | 5 eae | 6? + 1—6"| ‘ (76) 
for n = 2, 3, 4 or 5. Since 
O=relv and ¢?+1=Re” 


jeter 


we have, as in the proof of Lemma 10, 
cos (w —ny) = {R? + 2" — |g? + 1—9" 7 }/2e"7R 
ee chase ok ph 1 (77) 
= 90( 24) Smee 


+1—-@=1, 
2+1—-A=6-, 
674+ 1— 6*#=— 6-4, 
C7+1—-A=——6, 
and @ and © satisfy similar identities. Thus 
|?-+1—6|-|g?-+1—g"2=1, 


for n= 2, 3, 4 and 5. Consequently when n has one of these values the- 


where 


Now 
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expression on the right hand side of (77) reduces to that on the right hand 
side of (76). This proves (76) and completes the proof of the lemma. 


Lemma 15. [} a < 62+ 1, then 


snd (x—@)|<6,. . 2... serie (78) 

eed Lai 5 Sera me as aera 2) 
and 

‘5+ )\! 

(242) So 


Proof. If a= 62 + 1, it follows by Lemma 14, that 


3 
cos (y—ny) < cos(w—ny) + 1-3d—1 Ja—o—1~3} (548) —di, 
for n = 2 and 3. Thus 
- — 2sin}$ (z+ @)—ny} sin 4 (y— a) 


\z 
<1:30—4|a—6—1/—4}(544) ah - (81) 
for n = 2 and 3. 


By Lemma 11 and the numerical values of w and wy we have 
324° << $ (y+ 0) < 330°, 
so that 
44° <4 (7+ 0)—2y< 51°, 
—96° << 4 (z7 + 0) —3y< — 89°. 

Thus 

2 sin {4 (y+) —2y} >2sin 44° > 1-3, 

2 sin {¢ (7 + @) —3y} <—2sin 96° < — 1-9, 
These results together with (81) imply that 

1-3| sin} @—o)| << 1-38—4 |a—G?—1 at (o) at, 


. 


so that 
3\t 
1-3| sin $(y—o)|+ 3 Jara] +4} (75°) —d} <1-38 


As the three terms on the left are non-negative, this inequality implies the 
inequalities (78), (79) and (80). 
Lemma 16. Jf ¢= 62+ 1 then 
PAST. ws. . (82) 
eee Sd, ew se ee (83) 
and 


pe) 22 82. Meet yy (84) 


a 
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Proof. If a= 62 +1, it follows from Lemma 14 that 


\t 
59) att . (85) 
for n = 4 and.5. 


By Lemma 12 and the numerical values of w and w we have 
326° <4 (y+) < 327°, 
so that 127° <4(y+o0)—4y+2a< 129°, 
—13° <4 (xyt+o)—Syp+2a<—11°. 
Thus 2 sin {4 (y + w)—4y} > 2 sin 129° > 1-5, 
2 sin {4 (y + @)—5 yp} <—2 sin 11° << —0°-38. 
These inequalities together with (85) imply that 


— 2sin {4 (y+o)—ny} sin 3 (x—@) 
<1:36—4le— G7 — 13 


z 
0-38 |sin 4(y—)| << 1-36—4 |a—8?—1| — 4}(542) —dt. 
so that 
esa cind ol le +4} (244) af <1-30. 


As the three terms on the left are non-negative, this inequality implies 
the inequalities (82), (83) and (84). 


Lemma 17, 
O%4 al raed, SG ee eee 
and hot PPh 13d... a) eee 


Proof. The inequality (86) follows immediately from Lemmas 15 and 
16. By those lemmas we also have 


tsincd, (y—co)| <u 3270.) 4) 08s. a ee ed 
and 5+6 

eee do 68> oo ee ee 
Thus 


e+ 1 
5 + 6\3 5 + 6\! ey 
|e (FEY +02) 
' 5+6 i 5. ht 
a egies alae 


7:86 
<2. 66+ (gry) a+) Mit gece 
< 2:66 + 1°35 (++ 
<7 6, 


) 
-n=|é-(oh) 
| 


ey 
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Now, as in Lemma 10, 
| B—2—1 |? = |deix_Rel» 2 
= {d’?+ R?—2dR cos (w—y)} 
=(d—R)? + 4dR sin? } (w—y) 
< 496? + 4dR (3° 7)?62 
< 196*-+ 4R(R +76) (3: 7)78? 
< 169 6?, 
Consequently 
|B—¢?—1| < 138. 
7. In this section we prove that a= 62+ 1 and B= 2+ 1. We write 
coe SON ae 
Pa <8 Be” 


{i ——— 


(90) 


and if | | >0 we write 
Pee as NIT ACN ance (ON) 
where 0 > 0 and 7 is real. 


Lemma 18. For every rational integer n, we have 
| (a"(—6*"— 1) (6 (—9?)"—1) (6 (—9"—1)| > 1-4... (92) 
Proof. It is clear that the number 


ee ee 1) pe NOE Drees 
+i] e+] [ee ae i) 


is an integer of k(@) for all rational integers n. So, by (17) with € = &n, 


ier x al CG PP Saal x i-eey 
(« +1 —1) ( GR 1) FEE rae 


for all rational integers n. Thus, by (90) 


eg ee a |e 


so that, using (16) and (56), 


&,= 


a1 


| (a’— (—0?)") (B’— (—92)") (6’— (—9)") | 
> (1—2) |(1—a’) (1— f’) (1 6’) 


ae (0? + 1) (G?-+ 1) (9? +1)| 
[ap B| 


(94) 
==> = sa 4 2 F 2 1l— 

for all rational integers n. Replacing n by —n and using the identity 

@0=1 we obtain (92) from (94). 
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Lemma 19, 
OO? es i ae ae eee ee 
Proof, We have, by (90) 
Pe 
Cee eee aeotl ae 
So, by Lemma 17, 
? FESO 7°80 2 
[ot So ee aaOe,) 4 eee ea) 
and | 
136 136 136 
[Bl< =x Bot == {002 "5; Faas) 


(Go i 130 Reis ee 
Suppose, if possible, that | a’ | > 0. Then since 
BO, 


it is clear that there will be two consecutive non-negative rational integral 
values of n satisfying 


6-4 <|a’(—@)"| <1. 

So we can choose a non-negative rational integer n satisfying 
6-*< a’ (—O7)" <1. 

For this integer we have 


a ey = 1a. 


and 
| (B’ (—?)"—1) (6’ (—?)"—1)| < (1 + | A’ |)? << (1 + 108). 
Thus 
| (a’ (—6?)"—1) (6’ (—?)"—1) (B’ (—9?)"—1) | 
< (1—@-4) (1 + 106) ~~ < (98) 
<i 20-* oey, 
Since 
225< 0-4, 


the inequalities (92) and (98) are incompatible for this rational integer n. 
This contradiction proves that a’ = 0. Now (95) follows from (90). 


Lemma 20. 
Beagrie ene 
Proof. Suppose, if possible, that p’ +H 0. Write 
B= 2y—w = 999206... = = 2 = 251100) 
so that 


—p? = @-1 eft, 
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Then it ides from Lemmas 18 and 19 and (91) that for every rational 
integer n, 


1—36<|(6’ (—¢*)-"—1) (6’ (—¢?)-»—-1)| 
= | (2 O% ef U2 _}) (0 Or e~tU—nw)_—_)|>. 5 . (101) 
= 1+ @? 6"—20 6" cos (A—n pw). 
Now, since @ >0, there will be three consecutive rational integers n 
Satisfying 
6-3 cos 86° < 0 6" < cos 86°. 
Further, it follows, from (100), that for at least one of these three con- 
secutive integers 
cos (A—n pu) > cos 86°. 
For this integer n, 
206" cos (A—nu) — 0? 67" 
= 09" {2 cos (A—nu) — 0 67} 
> 8 cos 86° {2 cos 86°—cos 86°} 
= 0 cos? 86° > 0- 001 > 44, 
contrary to (101). This contradiction proves that f’ = 0. Now (99) 
follows from (90). 


8. Lemma 21. If a, b, b are of the form (9), where £4, 4, n, are con- 
jugate integers and o, t,t are conjugate units of the fields k{@), k(Q), 
k () then there is no integer & of k(9) with conjugates y, 7 satisfying 


| (€—a) (7—b) (@—6)| <3, 
but there are an infinite number of integers & of k(@), with conjugates 
n, 4 satisfying 
\(€—a)(y7—D) G—b)|=—4.. . . . . . (102) 
Proof. It is clearly sufficient to prove this lemma in the special case 


when 
pete l eer 
eee Pt a 


Then 


| (€ —a) (n— 5) Gj —b)| = | ea) ent (-a73)| 


ee 7 }) 
—| N@+1) 


=$|N(E(? + 1)—1)]. 
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Now for any integer & of k(9), 
N(E(@2+1)=N@+INE=SNE), 
so that either 
N(E(62+1))=0, or |N(E(6?+1))| 25. 
‘So there is no integer é of k(@) satisfying 
&(6?+1)=1. 
Consequently 
|(E—a) (9—b) @—B)|=4 |NE@+N—-I| Ss 
for every integer € of k(0). This proves the first part of the lemma. 


But, when £— én, where n is any rational integer and &n is the cor- 
responding integer of k(@) given by (93), 


| (En—a) ( eae —b)| 
= = (25 zy a pe 
GP 1) Nal ea +1 
= iN ee Vie 


This completes the proof of the lemma and the theorem. 


University College, London. 


Mathematics. — SoMMERFELD’s Polynomial Method in the Quantum 


‘Theory. By A. Rusinowicz. (Communicated by Prof. H. A, 
KRAMERS. ) 


(Communicated at the meeting of February 26, 1949.) 


§ 1. Scope of the paper. A. SOMMERFELD 1) has devised a method 
which enables us to treat many of the most simple but also most important 
eigenvalue problems in a very easy manner. His polynomial method starts 
with a splitting of the eigenfunction f of the originally given differential 
equation into two factors E and P 


The factor E takes care of the convergence of the normalization integral 
and of the fulfilment of the boundary conditions. The factor P however 
is assumed to be a polynomial, so that it does not disturb these properties 
of E. We suppose moreover that P is a solution of a differential equation 
of the second order with a recurrence formula containing only two terms, 
that is one of the form 


2 
PTA OB Nx 2 YA, 1 Be 3) D0, oi (2) 


(Az + By x*) x? 
dx 


A; and B; 2) are constants and h is a positive integer. 

In the present paper we will indicate the eigenvalue problems which 
can be treated by SOMMERFELDs polynomial method. In § 2 we suppose 
that the polynomial P in (1) is a solution of (2). Then a relation con- 
necting the coefficients of the differential equation (2) and the original 
one determines the general form of the potential V which appears in the 
original , differential equation. But the function V obtained in such a 
way contains the eigenvalue parameter j and the coefficients A; and B; 
which depend generally on /. If V has a real physical meaning it can 
however not depend on /. That means that there are some relations 
between 2 and the coefficients A; and B;. Further relations follow if we 
require that E guarantees the fulfilment of boundary conditions. All these 
equations determine not only uniquely the potential function V but settle 
also the eigenvalues. 

From this point of view we treat in the following sections some special 
problems completely. In § 3 we deal with the radial functions of the 
spherical symmetric case. In § 4 we discuss the differential equation of 
associated Legendre functions. In § 5 we start from the same differential 
equation as in § 4, but apply the polynomial method after a linear trans- 
formation of the independent variable. In this way we obtain the differential 
equation of JACOBI polynomials. 
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§ 2. The connection between the original and the polynomial differ- 
ential equations. Going over from the “original” differential equation 


f+ 2ef + bf=0 | ea ee (1) 
with the aid of (1.1) 3) to the differential equation of the polynomials 
P'352qP' 4 aP=0. 2G 2 ee ee 


we get between the coefficients of both the differential equations the 
relations 


a==te, | ping) Sgt a ee 
se hs E’ 
boa t2ag tb. <-l. a el oe A (3b) 


Eliminating E we obtain a relation between the coefficients of both the 
differential equations (1) and (2) 


at +t af — ba? +a =p. Se 4 ee 
We denote this expression in the following considerations by S. 

We assume that the coefficients of (1) are real numbers and that (1) 
is the differential equation of an eigenvalue problem. (1) is then selfadjoint 
and has the form 
d d 
Aes —qf>-defo=0. 2 & 2 5 
4 denotes the eigenvalue parameter and o(x) the density function. It appears 
in the integral 


Xg 
Ji feds < eh tne ag eel 


which decides in case of discrete eigenvalues about normalization and ortho- 


gonality of eigenfunctions. x, and xg are the boundaries of the fundamental 
interval. 


Comparing (1) and (5) we obtain 


ip 

be a (7a) 

—A 

pba 2S. 
: (75) 
According to (4) and (7) we can represent S in the form 

oS) lee 

s=1(}p ~42 +q-ie). NS <r Pity (8) 
so that q becomes 

g=pS + te +4219 9 
re Te oe (9) 


300 


The form of the expression S which appears in (9) is known because we 
suppose that the differential equation for P is given by (1.2) and therefore 


ee ie iy Ba Ao t Bo x" 
x A, + B, x" ~ x7 A,+ B, xt 
so that according to (4) 
S=a? po’ — pa (Art Bi x4) [Ai— A, + (B,—(h +1) By x!) 
x? (A, + B, xi)? 
ah Ay + (By—h dy) ce 
x? (A, + B, x!) 
Our final result is: If the differential equation (5) with given p and q 
is solvable by f = EP, P being a solution of (1.2), q must be of the 
form (9) where S has the form (11). 
We apply this proposition in cases where we can split the SCHROEDINGER 
equation 


(10) 


(11) 


2m 


Ay+x(&—V)y=0. emed Cade pee A Y 


into a number of differential equations of the form (5). Both p and © are 
then completely determined by the coordinates used for the separation of 
the variables and q contains generally an expression given by the potential 
function V. By (9) and (11) are settled the forms of the q's of all these 
differential equations therefore also the form of the potential function V. 

A more exact determination of V we obtain by the demand that the 
coefficients of V can not depend on the eigenvalue parameter 2, if V has 
a real physical meaning. This takes place if the expressions (9) for q do 
not depend on / #). But this can be fulfiled only if the coefficients B; which 
appear in S are functions of J. To find out this dependence we can expand 
the expression 


adr to Fok Fy kale © watlies (18) 


contained in (9) in a power series in powers of x—xy (x arbitrary) 
eventually after multiplication with a function of x. On this occasion we 
also find, that only for distinct values of h the expression (13) can be 
made independent of J. 

Further conclusions as to the A; and B; and so as to the potential V 
we can draw from the boundary conditions for the eigenfunctions. As a 
rule the fundamental interval in the quantum theory in bounded by two 
singular points of the equation (5). In such points the solution has a 
tendency towards becoming infinite. From the mathematical point of view 
it is the task of the boundary conditions to suppress this tendency. 

If we confine ourselves to the discrete eigenvalue spectrum we must 
claim that the integral (6) is convergent for any two eigenfunctions fy 
and fs, since otherwise we can not speak of their normalization or 


orthogonality. 
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Also in case of singular boundary points x; and x we prove the ortho- 
gonality of two eigenfunctions f; and f2 with the eigenvalues 2, and A, 
by the aid of the well known relation 


(4-4) fedx= lim ®(:)— lim O(x) .7. 2 04) 


X> X, 


where 
&0)=p (a) | A) BA — Ane |. «5 


To apply (14) in this case in the usual manner we have to claim not only 
the convergence of the integral (6) appearing in (14) but also the 
existence and equality of both the limits 


lim ©®(x)== lim @(x).. 7. a Se 2 UG) 


X> Xy X > Xe 

Physical arguments can stipulate to claim more, e.g. that (16) is zero or 
the eigenfunctions are periodic. But also mathematical motives can do it if 
e.g. the eigenfunctions have to form a complete set of orthogonal functions. 

Both factors E and P which form f can facilitate or make more difficult 
the fulfilment of the boundary conditions, P is a product of x? with a 
“genuine” polynomial, beginning with a constant and consisting of integer 
powers of x". The exponent o is given here by the determining fundamental 
equation 


o(o—1)A,+20A,+A,=0. .... . (17) 
But also the form of E is determined. According to (3a) we get 
E= exp ({ (a—a) dx). 
From (7a) we find 


const 


*/o 


exp —fadx= 


whereas exp { adx may be evaluated by the use of (10). Carrying out 
this calculation we have to distinguish three cases, according to the dis- 


appereance or non-disappereance of the constants A, and By. We get 
without an insignificant constant 


(I) in case of Ay #0, By 0 
E=p-" x4l4: (A, + Bz xy 5 eed be i 3 pel heal 
(II) in case of Ap = 0, Bo 20 


Eze poi chile Fons ee ee 
(iil}= in*case of Ap 20; Ba == 0 


1B, xh 
E= ph xAldige® 4” gy ee (18c) 
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The role played by the singular and zero points of @, p, P and E in the 
fulfilment of the boundary conditions we shall discuss when considering 
the different special cases. 


To these relations we have to add SOMMERFELDs condition of breaking 
off the power series 


(6 +n) (6 +n—1)B,+2(6+n)B, Bg 08 Fer 9} 
which makes of P a polynomial and determines the eigenvalues of 4 in 
their dependence of an integer n, divisible by h. 

Finally we indicate a very useful property of the polynomial equation 
(2). It does not change its form, if we multiply the polynomial P with a 
given power of x. Putting E = x” we get from (3) 


vy v(vy+) 


vy 
oie beakers iat Seat, teeta 5 3 


We obtain therefore for f = x” P the differential equation (1), in which 
the coefficients a, b are of the form (10) of the coefficients a, p. 

In both the cases (I) and (III) in which Ay ~0 we shall use the 
abbreviations 


§ 3. Spheric symmetric field of force. Splitting off from a Schrédinger 
eigenfunction a spherical harmonic we obtain for the radial function R(x) 
the differential equation 


Gg t|ee—ve—Fidty]R=o. . (1) 
i.e. the differential equation (2.5) with 
p= ee, AaewS, geen Vi xb i(l4 1). 6 oe (2) 
V is here the potential function and / the azimuthal quantum number. 
From (2) and (2.11) we obtain V in the form 


oS 241) 
ee ee gy eee eg Ee eS (3) 
We assume that the fundamental interval is given by O0<x< +0 and 
use for our considerations the factor E first. 


To make the normalizing integral (cp. (2.6) ) 
Renn dee a, 
0 


convergent at its upper limit we have to assume case III i.e. 


Pee sy ss 65) 


According to (2. 18c) and (2) E becomes 
ee ie 9 
E=x4 eh 


By hp 
A 


(6) 
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and the convergence for x > + © requires B}/Az <0. Finally we assume 
Avpe'A; . 6 eee eee 


to unite the x-power from E with the polynomial P. 
For further considerations we use V. According to (2.11), (5) and (7) 
the expression S is given here by 


STH ay + ((h + 1) by by) 8 + BE th 
To free V, Eq. (3), from 6 we have here only both the possibilities 
h=1 or =—2. For h = 1 the potential becomes 


iC=2 Czy 


y= AA G5, = 242 od ee 
x 
where the constants 
1 
c-2= — 7 (40 + H(I-+ 1). < 9 “seen ler a) 
1 bi 
C1, = on (2 b, —— bo) an 6 (9b), Co — & + oy eer me (9c) 


are independent of &. Therefore we obtain for V a Coulomb potential 
superposed by a potential inversely proportional to the square of the 
distance. The coefficients ci in (8) are arbitrary because their dependence 
on Aj, B; does not imply any connection between them. 

To the potential (8) belongs the RYDBERG formula. From (2.19) and 
(5) follows 


bose = 2 by (nie) a Se eee 
by 


and hence, in accordance with (96), c_, sae (n -+o6+1) so that we 
obtain from (9c) in fact the RYDBERG formula 


oad (c_1)? 


es Swe ie ee 


6é=— — 
4 (n+o+1) 
Supposing further, that the Coulomb potential has the right constant 
c_, = —e?Z we obtain in (11) the RYDBERG constant. Eq. (9a) not yet 


used determines o and hence the RYDBERG correction. From (2.17), (7) 
and (9a) we obtain 


a(o+ 1)=—ag=i(l+1l)+xer.. . . . . (12) 
In a pure Coulomb field, there is c_» = 0 and therefore oc =I (for 
o = —l—1 the normalizing integral (4) is not convergent) so that we 


obtain the BALMER formula. 
In case h = 2 we have 


V= 4 oe xtte 


x? 
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where 


1 
Cg — — (ay +E +-1)) Sees pee at's... (13a) 


1 2 1 
qi Oe as -. {13b), Co = & + — (3b, — by) <a (bac) 
The potential V is consequently given by a superposition of an elastic 
potential and of a potential inversely proportional to the square of the 
distance. 


The dependence of & upon the quantum numbers we obtain from (10), 


(136) and (13c) 


6=—22(ntotH+o=2|/(ntorH+a. . (14) 


The positive sign of the square root is determined by (6). 
If to the pure elastic field of force corresponds a frequency w (in 


m are He 
2a sec), we have to put cop oy w2. Like in case h = 1 the constant o is 


given by (12). . 
In case of a pure elastic potential there is c_) = 0 and therefore o = 1. 
For the eigenvalues of the spatial harmonic oscillator we obtain then 


S=(n+l+s)hota. . . ... .~ (15) 
The general case (14) we can conceive now as (15) with a RYDBERG 
correction o. The constant cg is in all the formulae of Legendre functions 
arbitrary and we can put cg — 0, if V is normalized as usual. 


§ 4. The differential equation for associated Legendre functions. To 
have an example of an eigenvalue problem in a finite fundamental interval 
we generalize the equation for associated Legendre functions 


ge (1-299) + (tp) K=0 


dx dx 1 jes 
putting V(x) for m2/(1—x?2) 
ga(0 29) + U-V (a) K=0. - Saree) 


We obtain hence the differential equation (2.5) with 
ee ee ee a tee ee see ~((2) 
From (2.9) and (2) we get therefore 


Vat S(i—x*) +4. Wes ste 
For the fundamental interval we choose —1<x< +1 and use first 
for our considerations the factor FE, given by one of the Eqs. (2.18). 
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The zero points of p, which according to (2) are situated in x = + 1, 
endanger the convergence of the normalizing integral given according to 


(2.6) and (2) by 


+1 . 
[K'K daze & ae ee 
=) 


a 


because p-! appears in the factors E in all the cases (2.18). This danger 
can be eliminated only in case I, Eq. (2. 18a), where 


Act, Bp 0 a ae a 
and where according to (2.18a) and (2) 
1 (BA) 
pe ht (Aa + Ba 2th A 
(1 — x?)'le 
To avoid E disturbing the convergence of the normalization integral (4) in 
the endpoints x = +1 of the fundamental interval and to secure the 


finiteness of the expressions (2.15) the binomial A, + B.x* must be 
divisible by 1—x?. Hence follows that 


A, a B, Sy - We hey © uw bin avon race (6) 


and 


h = even integer. “". ss <2) 3) ee ee ee 
Both our demands restrict also the variability of the exponent of 
A, + Box"; the demand that the expressions (2.15) are finite goes further 
and claims in accordance with (6) that 


Apr SO i 


To transpose x4/4: from E to P we must put 


A, — 0. . . . . . . . . + . (9) 
Now let us use V. In accordance with (5), (6) and (9) we obtain from 
(2250) 
Soom —1+(bs +h+1)x® — ap+(bo—hby) x" 
(1 —x"t)? x? (1 — xh) 
Developing S in a power series we get from the expression S p appearing 


in Eq. (3) 


(10) 


Sp=S(l—x2) = — = ag—M (x2 xh) —(M-+EN) (228-2 2h) 4 


+ (M + 2N) (x3"#-2—23h) +... 
where M = ay + by —(h—1)b,, N = —b, (6, +h). 
If V is independent of J the coefficients of a development of V in a 
power series can also not depend on 4. Hence in accordance with (3) the 
coefficients of a power series development of 


S(l—2) 2050. | ee ee 


(11) 


must be independent of 2. 
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First we assume that h F 2, ie. according to (7): h = 4, 6, 8, Pree 
accordance with (11) and (12) then the values of a, 29 + 1, M,N must be 
constant. But from the simultaneous constancy of ag and ay'+ J it follows 
that also the eigenvalue parameter 2 is invariable. In such a case however 
our problem is not an eigenvalue problem. 

For h=2 we have to demand, that only the three quantities ao, 
ag—M + iand N are constant. From this and the Eqs. (2.17) and (2.19) 
we could determine the dependence of the eigenvalues 2 upon the quantum 
numbers. But we may come to this conclusion in an easier way. 

For h = 2 we obtain from (3) and (10) for V the expression 


A B 
V= , ee + x2 + + . . . . . . (13) 
where the constants are given by 


A=(b, +1? (14a), B=—ay (14b), C=i—(8?+ 5,4). (14c) 


If we put 6; +1—=—m and remark that from (2.18) it follows 
by = (0 +n) (0 + n—1) —2b,(o +n) we finally obtain from (14c) 
A=(o+m-+n)(co+mint+i1)+C. .., gle #445) 


In accordance with (8) and (9), m is here a positive constant. 
From (2.17) and (9) we obtain for o the relation 


Phe 3) oe ag 8 eg 3) 


By SOMMERFELD’s polynomial method we can therefore solve a slightly 
more general differential equation than the one for the associated Legendre 
functions. This last differential equation we obtain by putting B= 0. 
In this case is o = 0 or = 1. Remembering that n is an even integer 
(h = 2!) and hence n+ o an arbitrary positive integer we see that (15) 
represents the well known eigenvalues of the differential equation for the 
associated Legendre functions. 

To demonstrate by an example the simplifications caused by the trans- 
position of the x-power from E to P, we indicate the results arrived at 
without the supposition (8). Instead of (11) we obtain an expression in 
which is substituted 


@)— ai (ai—1) foray, M—ay; (2a: +2b,+h—1) for M, 
N—ay, (a: +2b6;-+hA) for N. 


But this does not alter the conclusion that h = 2. 
We obtain V in the same way from (13) but have to substitute 


A’ =(ai+6;,+1)?, B’= B+a; (ai—1), C=C 


for A, i, ©, 
Putting a, + 6; + 1 = —m so that we have again A = m? we obtain 
for the eigenvalues the expression 
A=(o+a;+m-+n)(o+a,tm+n+1)+C’ .. (15’) 


24 
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where o + a, is given by the equation Ree 
(ota,)(otaj—1I)=B. *. . we s. (16’) 


The equations (15’) and (16’) we obtain from (15) and (16) if we write 
o + a, instead of o. But this changes only the notation. 


§ 5. The differential equation of JACOBI polynomials. If SOMMER- 
FELD’s polynomial method is not applicable to a certain differential equation, 
we can try to give the latter a new form by a transformation of the in- 
dependent variable and then to apply this method. We expect to succeed in 
this way from the fact, that in SOMMERFELD’s polynomial equation (1.2) 
the zero point plays .a distinguished role which is after a transformation 
taken over by another point of the fundamental interval. That means: If 
we replace in a “given” differential equation the independent variable by a 
new one and regard such an obtained equation as the “‘original’’ differential 
equation (2.1) or (2.5) we can generally solve the “given” differential 
equation by the polynomial method for other potentials V, as by direct 
application of this method to the “given” differential equation. 

To verify this statement we use the differéntial equation of the associated 
Legendre functions (4.1) i.e. 


dit caaeteedh sichgh | 
Alias 1 Gyr) FE—V(x)) K=0 Se 
where we have denoted the independent variable by x’. Substituting here 
by 

MSS LS AP a oon ee 


the new independent variable x, we obtain the differential equation 


d 


7 x(2—x) 9% 4 G@— Vay) K=0 


which we will consider as the. “original” equation of SOMMERFELD’s poly- 
nomial method. It has the form (2.5) with 


p=x(2—2), =l,. q=eVv i wate & 3) 
so that according to (2.9) and (3) the potential V has the form 


1 
aaah pee fin tod feted ac © USP os le ri) aoe ae (4) 


In the variable x’ the fundamental interval is bounded by + 1, in x it is 
therefore given by 0<x<2, 


For further conclusions we use first the factor E. We have to choose it 
in the form (2. 18a) to guarantee finiteness of the normalizing integral 


2 
LK R det 0 eee en 
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Otherwise p = x(2—x) endangers the convergence of (5) at its upper 
limit. Therefore we have to put A, 0, By + 0 and obtain 


r (BA) 
B= xA7! (Ar t+ Bo xt)t Bi 
a — @=x 


To guarantee the convergence of (5) for x= 2 we must suppose that 
the expression A, + Box" is divisible by 2—x, i.e. that 


1 
- B, — Qh A> Ps Be Rca or, Gy Ne (6) 
Further we have to assume according to the higher demands of (2.15) that 
l B, A, 1 wei’ 2 B, A, 1 : 
aft (BA) Tai 2(2-A)ob a 
Finally the removal of the x-power from E to P gives 
A, — 4 A:. Ses Ty oe Re eis se oh ow ees (8) 


To use V for further considerations we remark that according to (6) 
and (8) the expression S is given by ; 
gs——U+my’)(l+ny") a 1+py't 
Z 16 y?(1—y")? 4g7{l—9") 
where y = x/2 and 
| 
m=2'*15,,. n==—2*t b, —2(h +1), p= Moh hy) 
0 


Developing S in a power series in y we obtain for the expression 
Sx(2—x) appearing in Eq. (4) for V 


Sx(2—x)=4Sy(1— 9) =— (ao +2) (5-1) + My1—y!) + 


+ (M+ N) (y?*#-!—y") + (M+2N) (y34#-!— 34) +... 
where 
M=+—a—2 (b,—hb;), N=4 [2h—1-4+2"+" b, (2*+1 b+ 2h+2)]. 


To fix the value of h we can now use the demand, that the coefficients 
of Sx(2—x) +4 (cp. (2.13) and (4)) are independent of 4. For 
h=2,3,4... we have to claim that aj + 1/4, a + 1/4 +4, M, N are 
constant, so that 4 would be constant. 

For h = 1 the expressions 

Gg-+-4, a9 4+-4+M+1,N 
only have to be constant. But this means that 


@) = const, b; = const (9a), 4=—M + const = 2b) + const (9b) 
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4 can now depend on quantum numbers because bo is not constant now. 
To consider the case h = 1 in detail we remark that according to (2.11 i 
(4), (6) and (8) the potential V has the form 


x—2 ee 
where 


A=—8by (b,-+1)—2 (10a), B=—2ap (10b), C=4—2 (bp +b: +263). (10c) 
The relations (9a) follow also from (10a) and (105) and the relation 
(9b) follows from (10c). 

Reintroducing by (2) again the variable x’, we obtain V in the form 

_ A B _ x'(A+B)+A—B 
VS esi ot ace x’2?—1 

But (1) represents with this V the differential equation of JACOBI poly- 
nomials. It is therefore situated at the limit of the applicability of SOMMER- 
FELD’s polynomial method. 

In the quantum theory of a spinning symmetrical top we have to do with 
this equation with 


A=—4(t—v)?, Bo t(rtrP-. se 
where + and 7’ are positive or negative integers. 
Using that we have according to (2.19) and (6) 


(o +n) (o + n—1) — 4(0 + n)by—2bp = 0, 
we obtain from (10) the eigenvalues 
A= (o + n—2b,) (o + n—2b,—1) + C, 
According to (2.17), (8) and (10) we have to calculate o from 
os —ag=— 4.8.) Sa See 


Supposing especially the case (11) we get from (10) according to (7): 


+C, 


b, = —3 et and from (12) in accordance with the fact that o> 0 
(otherwise we would obtain for x’ = —1 an inadmissible singularity): 
feo ier A 


. Hence we get the well known result 


Asn ena 1) eG 


/ / 
where 1t* = anal is the larger of both the integers |r| 


and | 7’ |. 
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1, A comprehensive treatment was given by A, SOMMERFELD in Atombau und Spektral- 
linien, Vol. II, Braunschweig 1939, cp. p. 716. 
2. SOMMERFELD denotes our coefficients 2A, and 2B, by A; and By. 
3. (a, b) means equation b of séction a, 
4. Compare however (3.1) where q depends on 1, 


Biochemistry. — Elastic-viscous oleate systems containing KCI. IV. The 
flow properties as a function of the shearing stress at 15° and 
constant KCl concentration. By H. G. BUNGENBERG DE JONG and 
H. J. VAN DEN BERG. 


(Communicated at the meeting of February 26, 1949.) 


1. Introduction. 


An orientating investigation!) on the viscous behaviour of KCl con- 
taining oleate sols at 25° C, was performed with a viscometer which 
worked with constant hydrostatic head. The maximum shearing stress 
(at the wall of the capillary) was therefore nearly constant (of the order 
of 40 dynes/cm2) 2). As PHILIPPOFF 3) has shown that under certain 
conditions Na oleate sols may show ‘“‘structural viscosity’, it seemed safe 
to adopt for a further investigation of the viscous behaviour of the KCl 
containing oleate systems methods which allow of a wide variation of the 
shearing stress. In Part I, II and III of this series 4) the elastic properties 
of these systems have been studied at 15° C (at which temperature the 
damping of elastic oscillations is in general much smaller than at 25° C) 
and at a favourable KCl concentration (at or very near to the minimum 
damping at the given temperature). We chose therefore the same con- 
ditions for the present work and performed all experiments with ‘““Na- 
oleinicum, medicinale pur.pulv. Merck’, which is present in the KCl 
(+ KOH) medium in a relatively low concentration. 

The first impression these systems make on the observer is that they are 
gels. The definition of gels often includes the presence of a yield value. Our 
first aim will therefore be to investigate if a yield value can be detected, 
which involves the study of the flow phenomena at low shearing stresses. 


2. Measurements at very low shearing stresses. 


The viscometer used (see fig. 1) is in principle that of MICHAUD 5), 
whose characteristic feature is that slight level differences are obtained by 
immersion of a glass rod (in our case 0.55 cm diameter) into one of the 


1) H. G. BUNGENBERG DE JONG and G. W. H. M. VAN ALPHEN, Proc. Kon. Ned. 
Akad. v. Wetensch., Amsterdam, 50, 849, 1011 (1947). 

2) Dimensions of the capillary: R= 0.0425cm, 1 = 24cm; constant hydrostatical 
head A = 43.1 cm. ; 
3) W. PHILIPPOFF, Viskositat der Kolloide, Dresden (1942). Cf. p. 116, 122, 279—280 

and fig. 61 (p. 128) and fig. 162 (p. 280). 

4) H. G. BUNGENBERG DE JONG and H. J. VAN DEN BERG, Proc. Kon. Ned. Akad. v. 
Wetensch., Amsterdam, Part. I: 51, 1197 (1948); Part II: 52, 15 (1949); Part III: 52, 99 
(1949). 

5) F. MICHAUD, Ann. de Phys. 19, 63 (1923). 
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wide tubes (we used 3.6 cm diameter) and that the movement in the 
capillary of the system investigated, is observed in the microscope. 

For the viscometer used a level difference of 0.0239 cm is obtained by 
screwing the glass rod downwards or upwards over a distance of one 
centimeter. With the aid of the formula P = R.p/2.1 we obtain from the 
dimensions of the capillary (R = 0.0425 cm; 1 = 24 cm) and from the 


eS 


Fig. 1. 


density of the soap system (@ = 1.066) a value of 0.0169 dynes/cm2 for P, 
the shearing stress at the wall of the capillary for one centimeter immersion 
of the rod. 

Before starting a series of measurements the stopcock in the horizontal 
tube connecting the wide vertical tubes is closed and a certain quantity of 
the oleate system is brought into the left hand wide tube. With the aid of 
suction the capillary and the horizontal tubes connected with it are next 
completely filled with the oleate system. The apparatus is then put into 
its definite position, the capillary plus attached horizontal tubes 6) being 
immersed in a tray (standing on the stage of the microscope and provided 
with a glass window at the bottom) through which flows water of constant 
temperature 7). Now the glass rod is moved into such a position that after 


6) Only these parts of the apparatus need be at the desired temperature, because 
the very high viscosity of the oleate systems practically excludes all convection currents. 
*) Running tapwater was used, which practically did not change in temperature during 
the measurements (constant within 0,2° C) 


. 
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completely filling the apparatus it will be immersed some 5 cm-into the 
oleate system. This original position of the glass rod will be called position 
zero”. 

The levelling of the oleate system and the drainage of the walls of the 
wide tubes above the levels (which during the manipulations have become ~ 
wetted by the oleate system) is then obtained by opening the stopcock 
(wide bore) for a sufficiently long time (we took 2—3 hours). 

One is not sure, however, that after closing the stopcock exact levelling 
has really been accomplished (e.g. on account of very high viscosity of the 
system or on account of any presence of a yield value). 

But the position of the glass rod which corresponds to real levelling 
follows from the measurements themselves 8). 

The latter consist in measuring the velocity of displacement (v) of 
suspended particles in the axis of the capillary at various positions (p) of 
the glass rod. These positions are varied into both directions (reckoned 
from position “zero”’), so that a number of values p are obtained in which 
the system in the capillary flows to the left and another number in which 
it flows to the right. Providing these p and v values with + or — signs and 
putting them in a graph, a symmetrical figure must result, the symmetrical 
point of which lies on the abscis axis (position of the glass rod) and 
represents that position of the rod at which real levelling exists. Various 
forms of such symmetrical figures can be imagined, from which a few are 
given under a, b, c and d in fig. 2. The vertical, dotted line represents the 


Fig. 2. 


position of the rod at real levelling (which will in general differ slightly 
from position “zero’’). ° 

In a and b a yield value is present, in c and d it is absent (or smaller 
than the experimental error). In c the velocity of flow is proportional to 


8) This position must remain constant during a series of measurements. This will only 
reasonably be the case if the stopcock after filling the apparatus has been left open for a 
‘sufficiently long time to allow the drainage of the walls. Otherwise a gradual shift 
of this position will occur during a series of measurements, 
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the shearing stress, in a it is a linear function of the latter and in b and d 
this function is more complicated. 

Two oleate systems were investigated with the technique described 
above. They were made by mixing 50 and 125 cc respectively of a stock 
solution (10 gr. Na oleate + 500 cc HgO + 50 cc KOH 2 N) with 
75 cc HXO+ 75 cc KCl 3.8 N and with 75 cc KCl 3.8 N respectively. 
Practically the systems differ only in the final concentration of the oleate 
(0.45 and 1.14 gr per 100 cc). The electrolyte concentration (1.43 N KCI + 
+ 0.05 N KOH or 1.43 N KCl + 0.11 N KOH) is such that the damping 
of the elastic oscillations is near to its minimum value. The mixtures had 
now to be provided with particles, that would indicate the velocity of flow. 

We will discuss this problem together with our experiments on the 
0.45 % oleate system. In preparing the mixture it was thoroughly shaken 
and had to be left to itself for a considerable time to become approximately 
free from air, though some small bubbles still happened to be present. It was - 
now tried if such entrapped air bubbles could be used as “indicators’’ of the 
velocity of flow in the axis of the capillary. An ideal indicator should have 
the following properties: 


a. it shall not alter the properties of the oleate system; 


_b. the dimension of the particles (bubbles or drops) of the indicator 
shall be small in comparison with that of the capillary; 


c. its particles, etc. shall neither rise nor sink during the measurement; 


d. a sufficient number of particles should always be present in the 
microscopic image, so that one can select one in the axis of the capillary 
for each single measurement. 


For entrapped air bubbles only a. is fulfilled. Very small bubbles are 
absent or relatively rare (contrary to b. and d.), so that one must select 
larger ones, which, however, are all rising upwards (contrary to c.). Soon 
after the beginning of a series of measurements no air bubbles are present 
any longer in the axis of the capillary and all lie now pressed at the upper 
part of the glasswall of the capillary. . 

The use of air bubbles is therefore not possible for indicating the velocity 
of flow in the axis of the capillary. Nevertheless it was decided to start 
a series of measurements on such an air bubble. There were no indications 
of the presence of a three-phase contact (air — soapsystem — glass), 
so that it was supposed that the air bubble could freely move along the 
glass boundary. If this assumption is correct, we would get information 
on the mean velocity of flow of the oleate system in a zone?) close to 
the wall as function of the shearing stress and the results could be compared 


®) This zone of course depends on the relative diameter of bubble and radius of the 
capillary. It is further necessary that the bubble shall not alter its diameter during the 
whole series of experiments. In our case this requisite was nearly fulfilled, 
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with measurements in the axis of the capillary which were to be performed 
with a more appropriate indicator later on. 

The results are given in fig. 3 in which v is expressed in scale-divisions 
of the eye piece micrometer per 100 seconds. The order in which the 


Fig. 3. 


successive measurements have been taken is given by the numbers next 
to the experimental points. The figure obtained is of the type represented 
in fig. 2a, which would mean that there is a simple type of functional 
linear relation between rate of flow and shearing stress in the investigated 
range considered (0—0.07 dynes/cm?) as well as a yield value. The latter 
expressed in dynes/cm? can be calculated from half the horizontal distance 
of the two parallel straight lines through the experimental points, bearing 
in mind that every cm change in the position of the glass-rod corresponds 
with 0.0169 dynes/cm?. The horizontal distance being approximately 
0.35 cm change in position of the rod, we obtain a yield value in the order 
of only 0.003 dynes/cm?. 

But it may be seriously doubted if a yield value of the oleate system 
itself is really indicated by fig. 3. It is very well conceivable, that the 
assumption of wholly free movability of the air bubble along the glasswall 
(as we did) is not true and that a certain difference in pressure between 
the two sides of the air drop must be exceeded before the air bubble itself 
“yields” from its original position (pressed against the glasswall). 

The measurements of the velocity of flow in the axis of the capillary 
were performed by using paraffine oil as an indicator. We added two drops 
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of paraffine oil, which slowly dripped from a 5.5 mm diameter glass rod 
held vertically into 200 cc of the 0.45 % oleate system. This was shaken 
vigourously in order to emulsify the paraffine oil. The oleate system being 
full of entrapped air was then left standing till the next day, to become free 
from air bubbles. 

The system was turbid in consequence of the paraffine oil having 
emulsified into very small drops, which were extremely useful for indicating 
the rate of flow in the axis of the capillary. Indeed, paraffine oil is the 
nearest approach to the ideal indicator, fulfilling the requirements men- 
tioned sub b., c. and d. on page 366. 

As to exigence a. (according to preliminary experiments) the influence 
on the period of the rotational oscillation was found to be very slight or 
negligeable. The velocity of flow was measured with a stopwatch by 
checking the time necessary for a very small drop (diameter + 0.2 scale 
divisions) to move over a number of scale divisions of the micrometer, 
and by calculating from this the number of scale divisions per 100 sec. See 
table I (upper part) and fig. 4, from which it appears that: 


a. There is no indication of the presence of a yield value; 
b. the velocity of flow in the axis of the capillary is proportional to 
the shearing stress at the glasswall 1°). 


Comparing fig. 4 with fig. 3 (in which the same arbitrary units for 
the velocity of flow have been used) we perceive from the far smaller 
inclination of the straight lines in the latter figure, that the mean velocity 
in a zone next to the wall (in which zone the air bubble moved) is much 
smaller than in the axis of the capillary. Unfortunately we did not measure 
the dimension of the air bubble and of the capillary (the diameter of the 
bubble may have been in the order of 1/5 of the capillary diameter), so that 
no quantitative calculations can be made to examine whether this difference 
in rate of flow is compatible with laminar flow. Except for any future 
indications to the contrary it seems probable that in the range of small 
maximum shearing stresses (0 — approximately 0.07 dynes/cm2) the oleate 
system, although it shows marked elastic properties, behaves as a New- 
tonian liquid. 

We therefore conclude that our oleate systems, although they make the 
impression of being gels, can be described better by the term: “elastic 


fluids”. 


10) Fig. 4 deviates slightly from the simple scheme of fig. 1c, in as much the experimental 
points do not lie exactly on one straight line. Instead of it the points lie on two straight 
lines, with a slightly different inclination, meeting each another at the same point of 
the horizontal level at v =0. In connection with the order of the measurements 
(alternatingly a point with a positive v and a point with a negative v, and further 
choosing these points in such a way that their absolute values become increasingly smaller) 
the shape of fig. 3 can be explained by a slight gradual shifting to the left as has been 


discussed in note 8. The conclusions a, and b, in the text are of course not endangered 
by this circumstance. 


+] 
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TABLE 1. 
Velocity of flow (v, in scale divisions per 100 sec) in the axis of the capillary. 


System _| Immersion depth (p) of the glassrod in cm, relative to the original position 


iavestigated |, . 1, a 1). 1). 1 )n 2 Doo 
+4 | +3 | +2 |1+]+0.5/-0.5] -1 | -2 | —3 [| —4 
0.45%o oleate 6.85 (5.36 | 3.69 1.94| 1.04 0.65 1.55 3.32 —5,22 | —6.96 
system at (7.06 (5.45 | 3.61 | 1.90| 1.09|—0.69\—1.52,-3.27 |—5.18 | —6.96 
15.4° | (3.55 | 1.86) 1.03 |—0.61 —1.50 —3.27 
(see fig. 4) ) ) (1.88 ) ) | 
meen: 6.96 5.41 3.62 1.90) 1.05 —0.65 —1.52—3.29 |—5.20 | —6.96 


: 


—0.202 |—0.276 | —0.351 
—0.204 |—0.259 | —0.348 


1. 149/p oleate | 0.357 | 0.248 | 0.164 | 
. —0.195|—0.277 | —-0.349 
: 
: 


system at |0.350 0.241 0.158 
14.4° | 0.351 0.242| 0.161 
(see fig. 5) ) 


mean: | 0-351 (0.244 0.161 


: 
' 

: 
: 

: 

; 
: : 
: 
; 


| 


—0.200 | 0.271 . —0.349 


The results obtained with the 1.14 % oleate system on the whole are 
the same. Because of the much higher viscosity of this more concentrated 
oleate system the drainage factor becomes of great importance now. (We 
waited 3 hours for levelling with open stopcock). The influence of this 
factor becomes plainly visible during the measurements, for after each 
displacement of the glassrod it will now take a long time until the new 
position of the levels is exactly reached in consequence of the slow drainage 
of the walls of the tubes and of the glassrod after the latter has been 
moved up, and of slowly renewed positions-of the menisci against these 
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glass-surfaces aftere the glassrod has been lowered). The experimental 
results show this influence of the drainage factor (see table II and fig. 5 }. 
The apparently irregular positions of the experimental points around the 
straight line drawn can be accounted for from the above if we take the 
order of the separate measurements into consideration (indicated by the 
figure next to each experimental point). 

We perceive, for instance, that the greatest deviation of a point from 
the straight line occurs if the glassrod has been displaced over a great 
distance just before the measurement. See point 3 and point 6, which both 
show too high a velocity (as can be expected from the above). Taking 
these disturbing influences of slow drainage and of slowly renewed wetting 
into consideration, we feel justified in concluding that also in the case of 
the 1.14 % oleate system no yield value is present and that the velocity 
of flow is proportional to the maximum shearing stress at the wall. In the 
domain of very low shearing stresses (0O—0.07 dynes/cm?) this oleate 
system, since it shows very marked elastic properties and at the same time 
behaves as a Newtonian liquid, should not be described as a gel but as 
an elastic fluid. 


3. Measurements at higher shearing stresses. 


An investigation of the flow behaviour in a large range of much higher 
shearing stresses (order of 1—1000 dynes/cm2) showed that the simple 
flow behaviour we met at very low shearing stresses (section 2) is lost. 

In principle we followed the method described by PHILIPPOFF 11) and 
calculated from the experimental data the quantities V (= 4 Q/zR3), 
the mean velocity of flow expressed in cc/sec., and P (= R.p/2L), the 
shearing stress at the wall of the capillary expressed in dyne/cm?. 

A survey of the flow behaviour is then obtained by drawing a curve 
— the flow rate curve — through the experimental points in a diagram in 
which log V is represented as a function of log P. 

But for the differences mentioned below, the same device was adopted 
for observing the velocity of displacement of a drop of petroleum in a 
calibrated tube placed horizontally. Differences: a. only one capillary was 
used, and as a consequence a number of tubes of different diameters for 
the drops of petroleum; b. the viscometer was calibrated with a New- 
tonian liquid of known viscosity, from which measurement the radius R 
of the capillary (its length being known = 11.8 cm) was calculated (which 
value R = 0.0542 cm was needed for the calculation of the quantities V . 
and P mentioned above). 

The use of only one capillary, though not in every respect a happy 
choice 12), made it possible to exclude, that the peculiar shape of the flow 


11), W. PHILIPPOFF, Kolloid Z., 75, 155 (1936). 

12) It is of course not to be expected that with one capillary it will be possible to 
measure the rate of flow accurately in the large tract of P-values used. 

At very high P-values and low viscosity the rate of flow of the oleate system is so great 
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rate curve, obtained with oleate systems containing KCl, was somehow 
connected with the succesive use of capillaries of different diameters. The 
curves obtained in b. made it further possible to control that in our visco- 
meter the Newtonian calibrating liquids showed, indeed, a quite normal 
behaviour. 

Fig. 6 gives in a log V vs log P diagram the flow rate curves of the two 
calibration liquids used (60 % sucrose at 15° and a mineral oil purchased 
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Fig. 6. 


from and calibrated by the laboratory of the ‘‘Bataafsche Petroleum Maat- 
schappij’, Amsterdam 13) and the flow rate curve of an approximately 
0.6 % oleate system (containing 1.43 N KCl + 0.05 N KOH) at 15° C, 


that the correction for the kinetic energy (HAGENBACH) had to be made (see note 14). 
At very low P-values the rate of flow of the oleate system is extremely small and a 
capillary tube had to be used for the drop of petroleum. The displacement of this latter 
drop (read off on the scale divisions of an eye piece micrometer) was observed in the 
telescope of a cathetometer. Certain difficulties, which seriously interfered with the accurate 
measurements of the rate of flow of very high viscous oleate systems, will be discussed 
in note 16. 

13) The mineral oil had a kinematic viscosity of 4,424 Stokes at 25°. According to 
F. J. BATES and associates, Polarimetry, Saccharimetry and the Sugars, circular of the 
National Bureau of Standards, C440, Washington 1942, see table p. 671, the 60% 
sucrose solution must have a viscosity of 80,1 centipoises at 15°. The sucrose used was 
recrystallized and from the flow data obtained with the two calibration liquids, followed 
a viscosity of 79,6 + 0.7 centipoises of the 60% sucrose solution which agrees well with 


the value given above. 
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The curves (dotted ones) of the two calibration liquids, as is to be 
expected, are straight lines with a slope of practically 45°. 

The curve for the oleate system, although presenting a similar course at 
low values of log P, considerably deviates with further increase of 
log P14). The shape of this curve, as we will see in the next communication 
of this series, is characteristic for the oleate systems containing KCl, which 
show marked elastic phenomena, although its position in the log V vs log P 
diagram still depends on a. number of variables. 

The typical flow behaviour is perhaps more conveniently represented by 
fig. 7, in which the abscissae are once more log P values, the ordinates, 
however, giving the values of log P—log V. The latter quantity has the 
meaning of the logarithm of a viscosity 7 (expressed in poises) at the given 
value of log P. These curves will therefore be called viscosity-shear stress 
curves. 

We see from fig. 7 curve A that our 0.6 % oleate system behaves as a 
Newtonian liquid (log P—log V independent of log P) at low values of 


log. P- log V 


Figs 7s 


the shearing stress only. Obviously we have here the continuation of the 
simple flow behaviour we met in section 2, at much lower log P values. 

The viscosity of the 0.6 % oleate system in fig. 7 in this range is of 
the order of 20 poises (log P—log V = +1.30). However, with 
increasing log P values the viscosity decreases extremely strongly. At the 
highest P values investigated, its viscosity (of the order of 4 centipoises 
only) is still on its way downward. 

It seems probable that at still higher log P values the viscosity-shear 


14) Application of the correction of HAGENBACH had a visible influence on the course 
of the flow rate curve in fig. 6 at the three highest P-values used. These corrected values lie 
on the drawn curve, those that have not been corrected on the dotted curve 
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Stress curve will reach asymptotically an end value, which, according to 
PHILIPPOFF, may be denoted by 7. If we denote the highest horizontal 
level by 7, the pecularity of our oleate systems containing KCI consists in 
the occurrence of an additional level 15) situated between Ho and y, . This 
level may be denoted by m (i = intermediate). 

If it is permitted to interpret the viscosity-shear stress curve of the 0.6 % 
oleate system from the point of view of its structure, we obtain the following 
picture: At very low shearing stresses (0O—0.07 dynes/cm2, see section 2) 
there is a structure in the oleate fluid which allows a steady flow propor- 
tional to the shearing stress. This is still possible at slightly higher shearing 
stresses (the 1 level in fig. 7) but above a certain value (in fig. 7 in the 
order of 8 dynes/cm?) this structure breaks down in two steps. In the 
first step the “viscosity” decreases from approximately 20 poises to 
approximately 0.7 poises, in the second step, beginning in tigs “7ecat 
approximately 100 dynes/cm?, it decreases to values which are only a few 
times larger than that of the corresponding KCI solution (order of 0.01 
poise) without the oleate. It is further important to remark, that the two 
steep curve branches and the intermediate level of the viscosity curve are 
well defined. One may switch over to other P values situated far apart 
and into both directions without influencing the position of the flow rate 
curve. This means that if we spoke above of a structure in the oleate fluid, 
this structure is not one which breaks down irreversibly. On the contrary, 
the structure postulated is one which in accordance with the prevailing 
shearing stresses, practically immediately breaks down to or rebuilds itself 
to certain equilibrium states. 

In this connection it seems of importance, that soaps belong to ‘“‘Associa- 
tion Colloids” and it suggests that the very high viscosity at the 7 level 
of rather diluted oleate systems is connected with large scale associations 
of the oleate molecules. Such associations in principle represent equilibria. 

It might be supposed that these equilibria are influenced by sufficiently 
large shearing stresses and that the above mentioned “breaking down or 
rebuilding the structure’ is connected with it. 


Fig. 7, curve B, gives the results of measurements with an approximately 
1.2 % oleate system at 15° C. The black dots represent measurements 
with the same viscometer used for the 0.6 % oleate system. At the low 
log P values the rate of flow of this far more viscous system becomes so 
small, that as a result of disturbances connected with the small radius of 
the capillary 16), it became impossible to obtain reliable measurements. 


145) A similar type of flow rate curve was found by PHILIPPOFF in 13,84 % Na-oleate 
solutions, containing NaOH and 1,62—7% m-Kresol at 20°. See his book (quoted in 
note 2) on p. 279—280, fig. 162. 

16) The entrapped air between the meniscus of the oleate system and the drop of 
petroleum in the capillary makes a sensitive air thermometer. Slight fluctuations in the 
temperature of the thermostate reflect themselves in an irregular displacement of the 


374 


That here too a 7 level (analogous to that of the 0.6 % oleate system) 
is actually present 17), follows from measurements (white dots) with a 
viscometer provided with a capillary of wider bore (R = 0.1062 cm). 

Curve B shows the same general characteristics as curve A so that we 
need not discuss its form again. Increase of the oleate concentration has 
obviously two consequences, a. the values of the 7 levels become higher 
(for curve B 7) = + 52 poises, ni = + 1.5 poises, 9, not reached, lies 
below 0.08 poises), b. the values of the (max.) shearing stress at which 
“the structure breaks down” are greater (first step: order of 20 dynes/cm?; 
second step: order of 300 dynes/cm?) 18). 

Comparing the curves A and B more closely we perceive, however, a 
difference in the position of the intermediate » “level”. In B the experi- 
mental points do not lie on a horizontal line, but on one that is slightly 
sloping upwards to the right. 

We are here dealing with “shear rate thickening’, the viscosity increasing 
(+ 1.3 > 1.7 poises) with increasing shearing stress. This might indicate, 
that the “residual structure’ left after the first step downwards (from the 
49 niveau) is modified at further increase of the shearing stress in such a 
way that it offers more resistance to flow, before definitely breaking up 
(fall of the curve towards the 7, level). . 

This shear rate thickening at the intermediate “level’’ is obviously 
facilitated by an adequate concentration of the oleate, which point will be 
of importance for a theory on the internal state of our oleate systems. 

In the next communication, in which another viscometer will be used, 
we will meet with more examples of this shear rate thickening and once 


petroleum drop (standstills with time and occasional temporary changes in the direction 
of displacement occurred). This irregular movement of the drop does not reflect a property 
of the oleate system; indeed in section 2, where we used still lower shearing stresses, but 
observed the flow of the oleate system itself in the capillary, we always found completely 
steady flow. It is still possible to investigate highly viscous oleate systems with the 
petroleum drop method, but one must for that end follow the displacement of the drop for 
a long time. Yet such a determination may be wrong, if apart from the relatively fast 
fluctuations in temperature, there happens to occur a slow and slight displacement (e.g. a 
few 0.01°) of the mean temperature of the thermosatte during the single measurement. 
Therefore one ought to repeat such measurements many times, which would take many 
hours for the location of a single point on the 7 level. It is easier to use a viscometer 
having a capillary of wider bore, which reduces the influence of the said disturbances 
on the measurements considerably (the quantity of fluid flowing through the capillary 
per second being proportional to the fourth power of the radius). The viscometer used 


to measure the 7 level of the 1,29 system is, however, not suited to measure the whole 
flow curve with it. . 

17) From the results of section 2 the occurrence of an 7 level at low values of the 
shearing stresses was already certain (proportionality between rate of flow and shearing 
stress). 

18) The approximately 0,6 and 1,2% oleate systems used for these experiments are 
not wholly comparable in age and composition. It is therefore not allowed to compare 
them quantitatively, e.g. as regards the functional relation between yo and oleate 
concentration. 
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more we will find it clearly present at higher oleate concentrations (1.2 and 
1.8 %) and absent at lower ones (0.6 and 0.3 %). This absence might of 
course be only apparent, as the method used is not accurate enough to 
detect a very slight shear rate thickening. 


4. Certain particularities of the oleate systems containing KCl explained 
by the shape of their flow curves. 


If we have a relatively large vessel (e.g. 1L) partially filled with an 
1.2 % oleate system and try to pour out its contents, we have not the least 
difficulty in doing so. Certainly the fluid does not impress us as being 
extremely viscous. This is caused by sufficient shearing stresses being 
set up, so that the system does no longer behave as a fluid of + 50 poises 
viscosity, but as one of materially lower viscosity. When the fluid is 
poured out, our working point is no longer on the 7 level, but much lower. 

In Part I (see section 12) we described the appearance of a dimple in 
the surface of the oleate system soon after excitation of the rotational 
oscillation. It is followed by a slight elevation of the surface, which may 
last some seconds (see fig. 6g and A in Part I). Now we can understand 
why this elevation may last several seconds. During the vigorous motion 
attending the formation of the dimple and its disappearance, the oleate 
fluid does not behave as a very high viscous system (too large shearing 
stresses). 

A certain quantity of the oleate system is then brought above the 
horizontal level. As the large shearing stresses have meanwhile disappeared, 
a certain quantity of oleate system with the viscosity corresponding to the 
No level, lies a fraction of a millimeter higher than the surrounding surface. 
So the shearing stresses which try to level out this elevation are very small, 
and therefore even its edge may remain visible for a couple of seconds, 
before it rounds off. 


5. Summary. 

1. The viscous behaviour of a few oleate systems containing KCl which 
showed marked elastic properties, has been investigated at 15°, with the 
techniques given by MICHAUD in the range of very small shearing stresses 
and by PHILIPPOFF in the range of larger shearing stresses. 

2. In the domain of very small shearing stresses (shearing stresses at 
the wall of the capillary varying from 0—0.07 dynes/cm?) the oleate system 
showed a steady rate of flow proportional to the shearing stress. 

3. A yield value could not be demonstrated. Our elastic oleate systems, 
which at first sight may make the impression of gels, can for this reason 
be better characterised as elastic fluids. 


4, In the domain of larger shearing stresses (order 1—1000 dynes/cm?) 
flow rate curves were obtained, which show a very strong decrease in the 


viscosity (more than 500 times) in a two step process. 
25 
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5. After the first step downwards (from the 7 level) the viscosity 
retains over a certain range of shearing stresses, nearly the same value 
(the 0.6 % oleate system) or it increases slightly (‘shear rate thickening; 
1.2 % oleate system) before the second step downwards (towards the 
N- level) sets in. 


6. The sequence in which the shearing stress is varied does not change 
the position of the flow rate curve. If a structure is postulated in the 
original oleate system, which is broken down in a two step process, it must 
be one which, in accordance with the prevailing shearing stresses, allows 
a rapid and completely reversible rebuilding. In this connection some 
suggestions are made regarding the role of large scale associations of the 
oleate molecules. 

7. A few particularities of the elastic oleate systems containing KCI 
described previously find an explanation from the shape of the flow rate 
curves. 


Department of Medical Chemistry. 
University of Leiden. 


Biochemistry. — Elastic viscous oleate systems containing KCl. V1). 
Viscous and elastic behaviour compared. By H. G. BUNGENBERG 
DE JonG, H. J. VAN DEN BERG and L. J. DE HEER. 


(Communicated at the meeting of March 26, 1949.) 


Introduction. 


According to MAXWELL’s well known formula 7 = G X 4, the product of 
the shear modulus G (in dynes/cm2) and the relaxation time ] (in sec.) 
must have the meaning of a viscosity coefficient (in poises). 

In this communication measurements of the elastic and of the viscous 
behaviour have been performed side by side at 15° to compare the product 
G X 4 with the viscosity levels yo, n; and 1» Which according to the results 
obtained in Part IV of this series, occur in the viscosity shearing stress 
diagram of markedly elastic fluids containing KCl. 


2) M ethods used. 


Previous to measuring the flow behaviour each oleate system was inves- 
tigated as to its elastic properties, for which purpose we used the rotational 
oscillation in completely filled spherical vessels of different radii. For 
particulars see Parts I, II and III of this series. The method used for 
obtaining a survey of the flow behaviour was in principle the same as in 
Part IV, in which only one capillary and a number of tubes of different 
diameters were used to observe the rate of displacement of -a drop of 
petroleum. 

In connection with this choice it was no longer necessary to keep the 
original shape of PHILIPPOFF’s viscometer. The capillary and the other 
tubes, all of Jena glass, were therefore sealed together, and instead of the 
original spherical reservoirs of relatively small capacity, we used relatively 
wide cylindrical reservoirs (see fig. 1). 

The difference in level of the oleate system (or of the calibration 
liquids) could be easily read off by means of a cathetometer at the begin- 
ning and the end of each separate measurement. This mean hydrostatic 
pressure added to or subtracted from the air pressure applied (read off on 
a H,O or Hg-manometer) gave the pressure over the capillary. From this 
corrected pressure we calculated the value P, the (maximum) shearing 
stress at the wall of the capillary, for which purpose (P = aH the 
length of the capillary (L = 11,8 cm) and its radius must be known. This 


1) Part I has appeared in these Proceeding 51, 1197 (1948), Parts II, III and IV 
in these Proceedings 52, 15, 99, 363 (1949). 
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radius (R = 0.0716 cm) was calculated from the calibration of the visco- 
meter with the same Newtonian liquids of known viscosity as were used 
in Part IV. 

The other quantity necessary for the plotting of flow rate curves, V, the 
mean rate of flow, was calculated from the rate of displacement of the 
drop of petroleum by using the formula V = 4 Q/xR3. 


Bigzel. 


This time a wider radius of the capillary was chosen than in Part IV, 
as this diminishes the irregularities in the displacement of the drop of 
petroleum (cf. Part IV section 3, note 16), which do occur in the range 
of low shearing stresses with the very high viscous 1.2 and 1.8 % 
oleate systems. These irregularities could be further diminished by using 
in this range exclusively level differences between the oleate system in the 
cylindrical tubes and by closing the apparatus by a bent connecting tube 
as depicted in fig. 1. The drop of petroleum now separates two volumes 
of entrapped air (viz. between its own menisci and the menisci of the 
oleate system in the left and right reservoirs). Slight fluctuations in the 
temperature of the thermostate now effect slight variations in pressure 
on either side of the drop of petroleum which compensate one another 
partially, Although this device gave an improvement, nevertheless irre- 
gularities were still present, so that the measurements at low shearing 
stresses and very high viscosities showed a relatively bad reproducibility 
and the mean of a number of separate measurements may still embody 
a considerable error, 

To check the conclusion in section 3, based on the results obtained with 
this viscometer, we will use in section 4 a technique by which the quantity 
of fluid flowing through the capillary is measured directly. 
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3) Comparison of the values G. 2 with No mi and 7,. 


As in all the preceding parts of this series, we will also restrict our- 
selves in this section to a KC] (+ KOH) concentration, which at 15°C 
coincides or lies very near to the minimum damping of the elastic oscil- 


lations for the preparation of Na-oleate used (Na-oleinicum med. pur. 
pulver. Merck). 


TABLE I. 


Measurements of the rotational oscillation with 1.8, 1.2, 0.6 and 0.3% oleate systems 
(containing 1.43 N KCl -+ 0.18 N KOH at 15°C). 


Na oleate} R 
g/100cc | (cm) 


G 
(dynes/cm?) 


Gs<a 


(poises) 


7.46 | 50.0 deat 105193) 5741 99.8 


1.8 4.92 | 57.8 3.57 | 0.137} 3.57 | 2.61 ( mean | 99.7( mean 


4.12 | 65.3} 2.97 | 0.124] 2.97 | 2.40( 2.63 | 101.0( 100.2 oe 
7.46| — 5.40 | 0.201| 5.40 | 2.69 100.2 
7.46 | 35.3 | 8.44 | 0.348] 8.43 | 2.42 ANA 

1.2 4.92 | 39.5 | 5.60 | 0.246] 5.60 ise s0.of mean 97 
4.12 | 41.6 | 4.55 | 0.196] 4.55 | 2.33)" 43.0) 41.6 
7.46 | 18.8 | 18.09 | 0.755] 17.96 ~ 9.06 

0.6 5.01 | 19.0 | 12.35 | 0.757] 12.26 wa 8.77{ meas Ls 
4.12 | 19.1 | 10.03 | 0.763) 9.96 eet 8.98) 8.94 
5.65 | 7.7 | 35.48 | 1.105) 34.94 a2 237 
5.03 | 8.0 | 31.53 | 1.105 | 31.05 = CE see 

0.3 | 4.12] 8.4 | 24.48 | 1.112| 24.11 os 1,53) 4 49 
3.18 | 9.0} 20.75 | 1.105| 20.44 = : 
1.70 | 9.4} 10.03 | 1.109| 9.88 = 1.55 


In a first series of measurements we investigated four oleate systems 
of different concentrations (1.8, 1.2, 0.6 and 0.3 g per 100 cc 1.43 N 
KCl + 0.18 N KOH). The results of the elastic measurements (rotational 
oscillation at 15° C) have been collected in Table I and represented in 


rs 
3 %oleate /” 


03 % oleate 


o——— 0 —— o> oo — 


° 


é /0 


a 
ve y er 
2 


06% oleate 


a5 


380 


fig. 2 and 3. In accordance with the results obtained in Parts II and III 
of this series, we find that for oleate systems containing 0.3 and 0.6 g 
per 100 cc (see fig. 2B) A is independent of the radius of the sphere and 


mes 


sec 


12% oleate 
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18 % oleate 
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TABLE II. 


Flow measurements at 15° C of the oleate systems of Table I. 


Te dleale 


ra 


1.8 g. oleate per 100 cc. 


1.2 g. oleate per 100 cc. 


0.6 g.-oleate per 100 cc, 


0.3 g. oleate per 100 cc. 


log V | log P| log P/V | log V | log P| log P/V | log’V | log P | log P/V | log V | logP | log P/V 
O42 0.15) 2501 OF56=22) 0:29 1574 O97 2) OSS iad rs 0.99—2/0.50—1] 0.52 
0.53—2] 0.50 | 1.97 0.75—2| 0.56 | 1.82 0: 10—1), 0.36 | 1.26 0.33—1/0.81—1| 0.49 
O20 ==2)) O76) || 2.05 O59 1 1181259 022510. 5151526 0.55—1/0.00 0.45 
OIG ——2 ele OOM 2. iil 1275 1273) 10. 9810s fe O82) pala tO O5St— 1017 0.35 
OF 41— 1) 1237 1)" 1.96 2,03 1.88 | 0.86—1 |0.99 1.05 | 0.06 0.35.7 10.58 0.03 
Oa G1 1.57 1581 2°35 2.26 | 0.91—1 | 1.53 1.21 | 0.67—1 |0.69 0.47 0.79—1 
0.38 1.78 | 1.40 252 2.53 20201 1.86 1.51 | 0.66—1 10.83. 10.53 0:70—1 
1.99 2.05 | 0.07 eT eal 2,66) 095-1 2,44 1.76 | 0.62—1 |1.24 |0.68 Ol45—1 
XGA 2.29 | 0.03 3.05 Dl DN One 4a i240 1.92) |0.52-—9F | 1960! “10.84 0.24—1 
2.43 Zao le O,.09 3.76 S07 aOgo kno 2-06: | 0.532 il 88 Vb. Oo O;12=1 
PARI Pt sey || (Oh es 4.24 3,31 | 0.07—1* {2.89 223) (0,59 —) eee BLO Ois 1 
2.68 2.86 | 0.18 3.46 2.09 | OSs 2r020 Ge 0, 051% 
2.69 2.88 | 0.20 4.00 289) 0589=2* 1-3. 15 eae oG 0..81—2* 
3.61 Son Omae——l . 4.66 3236) 0.7025 3 s50ee esos 0.73—2* 
4.18 3.43: | 0,.25—1* = oylloh <aPAR As} 0, 60E=25 


log 79 = 2.02 (mean of first 
five values of log P/V) 


log 7; =0.03 0.20 
LOGE Hing <a Oe = 


log 79 = 1.78 (mean of first 
two values of log P/V) 


log 7; =0.86—-1-> 0.01 
log) Hen < 07207 —1 


log 79 = 1.23 (mean of first 
three values of log P/V) 


log ; =0.65—1 (mean) 
LOGia One? 


yo = 105 poises 
7; =1.17 1.6 poises 
N «» <.0.18 poises 


After application of the 


*” 


* 


no = 60 poises 
n; =0.7 1.0 poises 
" o@ <.0.12 poises 


correction of HAGENBACH. 


No = 17 poises 
7; =0.45 poises 
 @ <.0.05 poises 


log 79 = 0.49 (mean of first 
three values of log P/V) 


log 1; =0.13—1 (mean) 
log 7~< 0.60—2 


Ny = 3.1 poises 
7; =0.11 poises 
" @ <.0.04 poises 


ts ¢ 
logV 1 G01 poise 


+4 


+s 


+2 


t/L’ 


xy 
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that at higher oleate concentrations (1.2 and 1.8 g per 100 cc) A is pro- 
portional to this radius (see fig. 3B). 

In these two latter cases only the relaxation time J and consequently a 
value 7 = G- 2 can be calculated. Though for the purpose of comparing 
G-A with mo, i and 7, it would suffice to measure the viscous behaviour 
of the 1.2 and 1.8 % oleate systems only, we have investigated the 0.3 
and 0.6 % systems as well, in order to gain some insight into the flow 
behaviour as a function of the oleate concentration. The results have been 
given in Table II and represented in fig. 4A and 7A. 


a - 60% sucrose 
Be i a 15° 
fo ft of 
Lif 
, oe: tia P -calibration oil 
i Ve Sipe Bi 25° 
2 ; f Ve 


i . rs 4 a 
7 4 d C4 12% oleate 
- 18% oleate 


2% » YW 


=06% +» 
ae ae -2L A 


: blank 

« 2.00063 mol/l. Hexylale 
» 000188 : 

+ 0.00527 


baAw> 


MY A 
aff é 


A second series of measurements was started from the 1.2 % oleate 
system (which served as a blank) to which small amounts of n.hexyl- 
alcohol were added. The results of the elastic measurements (rotational 
oscillation at 15° C) have been collected in Table III and represented in 
fig. 5 and 6. 

It was found that the addition of hexylalcohol in the two concentrations 
which still permitted elastic measurements, does not change the pro- 
portionality between A and the radius of the spherical vessel. Therefore 4 
values and consequently 7 = G-4 values could be calculated. For the 
third concentration of hexylalcohol elastic measurements were no longer 
possible. It was nevertheless interesting to include also this latter oleate 
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TABLE III. 


Measurements of the rotational oscillation with 1.2% oleate systems (1.43.N KCI + 
0.18 N KOH) containing n. Hexylalcohol at 15°C. (For blank see Table 1.) 


a 


n. Hexylalcohol | R ids Tee A G Gx 
n aie 2 7 : 
mol/L (cm) faces nee (sec.) (dynes/cm?) | (poises) 
7.46) 33.6) 8.12 | 0.459 8.09 | 1.76 44.6 
5.04] 34.7] 5.45 |0.337| 5.44 | 1.61 ( mean | 45.1 mean 
0.00063 4,12] 38.1] 4.44 |0.288| 4.44 |1.54\ 1.64 | 45.2 43.5| 71.2 
3.18] 39.4) 3.70 |0.224| 3.70 | 1.65 38 .9(?) 
. 7.46) 10.2] 7.68 | 1.528] 7.46 |0.49 52.5 
0.00188 4,92|13.7| 4.95 | 1.026] 4.89 0.48 one 53.1 rr 25.9 
4.12/ 17.2) 4.22 |0.804| 4.19 |0.52 i DOat ° 


system in our investigations of the viscous behaviour. The results of these 
measurements have been collected in Table IV and represented in fig. 4B 
and 7B. 


Fig. 6. 


In discussing the results of the above two series of measurements we 
will first say a few words on the viscous behaviour. The flow rate curves 
of fig. 4A and 4B (log V as a function of log P) and the viscosity- 
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TABLE Iv. 


Flow measurements at 15° C of 1.2% oleate systems containing various concentrations 
of n Hexylalcohol.. 


0.00063 mol/L 
hexylalcohol 


0.00188 mol/L 
hexylalcohol 


0.00327 mol/L 
hexylalcohol 


vi i|i 
0.86—2 | 0.53 | 1.67 0.03—1 | 0.41 | 1.38 0.26—1] 
0.15—1 | 0.87 | 1.72 0.38—1 | 0.68 | 1.31 0.74—1 0.58 
0.62—1 | 1.11 | 1.50 0.57—1 | 0.81 | 1.24 0.17 0.42 
0.52 1.40 | 0.88 0.91—1 1.09 | 1.19 0.55 0.31 
1.94 1.72 | 0.78—1 | 1.47 144) 7 0-°97—T | 1.36 0.86—1 
2.25 1.92 | 0.67—1 | 2.24 1266) |) 0.421) 2514 ORS 7a 
2.43 2,12 | 0,569—1 | 2.82 1 1595 | O14 1272 05031 
2.68 2.42 | 0.74—1 | 3.01 Zo22) | OC2I—= 3 3225 0.82—2* 
2.97 Fhe HOF 5-1 | .3235 2.98 |0522—1 173.65 0.74—2* 
3.49 2.93 0.45—1*| 3.69 2.89 | 0.20—1*| 4.01 OL67—=2* 
4.19 3.26 | 0.07—1*| 4.19 3.21 | 0.01—1*| 4.46 0761 == 2% 
| 4.60 Bras Oro —2 
; 
log 1 = 1.70 (mean of first log y= 1.35 (mean of first | log No >0.71 
two values of log P/V) two values of log P/V) 
log 7; =0.66—170.75—1 | logy; =0.14—170.22—1 ; level degenerated 
log 1 <0.07—1 | log 1 <0.78—2 log No <0.61—2 
No = 50 poises | Ny = 22 poises 1) > 5 poises 
4; = 0.46 + 0.56 poises | 7; =0.14+0.17 poises 7; = level degenerated 
1» <_0.12 poises 1 «@ <.0.06 poises 1 <0.04 poises 


* After application of the correction of HAGENBACH. 


shear stress curves of fig.7 A and 7B (log V—log Pas a function of log P) 
show with one exception (the lowest curve in fig. 7B) the same characte- 
ristics viz. jo level, a tendency to reach an 7, level and the presence of an 
intermediate “‘level’’, 7, the experimental points on it lying either on a 
practically horizontal line or on lines that are sloping upwards to the right: 
shear rate thickening (as the corresponding curves in Part IV of this 
series, cf. fig. 7 in that publication). These characteristics were already 
discussed in Part IV, so that we need not discuss them again. 

It is remarkable that the curve deviating from this typical curve form, 
(at its lowest in fig. 7 B) belongs to an oleate system, the elastic properties 
of which could not be measured on account of the exceedingly great 
damping of the oscillations (at most only two turning points were present). 
The viscosity-shear stress curve characteristic of the markedly elastic 
oleate systems (marked on account of smaller damping) has now 
degenerated to one in which an 7: level is no longer present and it already 
strongly resembles the curve form found in many cases of non Newtonian 


viscosity, in which there are only two levels ( and 7,,). 


384 


We might further draw attention to a similarity on the one hand and a 
difference on the other between the influences of lowering the oleate. 
concentration and of adding n. hexylalcohol at a constant oleate concen- 
tration. Comparing fig. 7A with fig. 7B we see in both cases that the 
viscosity-shear stress curve as a whole is displaced downwards. In lowering 


5 A.18% oleate 
(=e Bs/2h% = 
+/ sei : C 206% ” 
Dis05%. 
SS 
ee 


12% oleate 


A . blank 

B - 2.00063 molfl Herylale. 
C. 000g =, x 
D. 0.00327 : 


eT, 0 / ra ) ¥ 


the oleate concentration, see fig. 7 A, the viscosity-shear stress curve is at 
the same time displaced to the left in a highly marked way. When adding 
hexylalcohol at a constant oleate concentration this displacement is absent 
however. - 
We have now come to the comparison of the values G- 4 with yo, 7: and 
N. and refer to Table V, in which we have collected the numerical values 
of these quantities for the four oleate systems in which this comparison is 
possible. | 
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We perceive from this Table that the product G. 4 does not coincide 
with any of the three 7 levels observed. The product is extremely different 
from 77,,, still very different from 7:, but it is of the same order of magnitude 
as 4o- 

TABLE V. 


Comparisson of G X 4 with My 2; and 7 — 
enters eee 


Composition 
Ga 2 n; Ui) 
Oleate Hexylalcohol F : ; 
(g. per 100ce) (mol/L) | (poises) (poises) (poises) (poises) 


4) The discrepancy between the values found for no and for G- i. 


The same order of magnitude of the values found for 7) and G- 4 gives 
rise to the question if these values are really equal and if only experimental 
errors-cause the discrepancy formed. 

The difficulties we met with in measuring high 7 values at low shearing 
stresses, and the discussion of these difficulties?) made us much more 
suspicious regarding the reliability of the average yj values than regarding 
that of the average G or / values. 

We decided therefore to compare a new G-/ with 70, the 79 value no 
longer being obtained with the aid of the petroleum drop method, but 
by directly measuring the quantity of the oleate fluid flown through the 
capillary of the viscometer. The latter was of the type as depicted in fig. 1 
and the measurements consisted in reading at stated intervals the position 
of the levels of the oleate fluid in the equally wide reservoirs, (which 
were in direct communication with the air) with the aid of a cathetometer. 
A larger radius of the capillary 3) (R = 0.1062 cm, L = 12.1 cm) was 
taken and a smaller diameter of the cylindrical reservoirs (= 1.773 cm) 
to reduce the time which the measurements still required (the series of 
seven given in the table below took two whole days). 

From the effective hydrostatic height 4) | A h/ln i) the density of 

2 
the oleate fluid, the change in position of the levels during the time elapsed, 
the dimensions of the capillary and the diameter of the reservoirs the 
viscosity was calculated with the aid of POISEUILLE’s formula. 

This method gave satisfactory results only if, after filling the apparatus 
and bringing about an initial level difference one does not start with the 


2) See section 2 and in Part IV section 3, note 16. 

3) The radius of the capillary was calculated from flow measurements with the same 
calibration liquids as were used in Part IV. 

4) See E. HATSCHEK, Die Viskositat der Fliissigkeiten, Dresden 1929, 
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measurements immediately, but checks the flow of the oleate fluid for 
several hours by means of a counter air pressure (to give the walls of the 
cylindrical reservoirs time to drain) 5). The counter pressure is then 
removed and some hours later one begins with the measurements. 

The apparatus functions well if after the intermission between two level 
readings, the level has risen in one reservoir quite as much as it has 
fallen in the other. 

As we had no more Na oleinicum “Merck” in stock and it could not be 
purchased either, we used Na oleate, neutral powder, “Baker” 6) for the 
renewed comparison of G-A and yo at 15° C. The composition of the 
system investigated was 1.2 g oleate per 100 cc (1.25 N KCI + 0.05 N 
KOH) 7). Table VI contains the 9 values obtained by the above method 


TABLE VI. 


Measurements of the viscosity of 1.2% oleate system (oleate “Baker” in 1.25 NKCI+ 
0.05 N KOH) at 15° C. 


Sar I Sn rs Mie ac St a Ns rh a te oe 
Mean P (dynes/cm?) | 32.10 | 27.97 | 23.98 | 23.26 |...47282. | 34556 
7 (poises) | 68.2) |:76.0 “| 74.9>ep 77.0" [76.68 [omer 


7, mean = 76.0 poises 
at a number of (mean) shearing stresses.at the wall of the capillary 8). 
One perceives that compared with the large fluctuations of the log mo 
values in the Tables I and II the fluctuation of mo is small now, so that 
we may trust its mean value much better. Table VII gives the results 


TABLE VII. 


Measurements of the rotational oscillation of 1.2% oleate system (Oleate “Baker” in 


1.25 N KCl + 0.05 N KOH) at 15°C. 


il 
10X> 


(dynes/cm?) |(poises) 


60.3 B13 1.118 0.112 2.79 46.6 
63.5 Sal 1.198 0,181 2.88 ( mean | 47.1 ( mean 
68.0 7.40 1.291 0.255 2.90 \ 2.87 | 48.1 \ 46.6 
61.240.2/5.35+0.03} 1.202+0.0034 | 0.184+0.003. | 2.91 44.7 


* Measurement several days later during the flow experiments. From this single 


measurement would follow G.4= 130 poises. One is however not sure to conclude 
from a single measurement to a real decrease of G.A with time. 


5) Cf. in Part IV section 2, notes 8 and 10 from which follows the necessity of taking 
a long time for this drainage. 

8) We are glad to express our profound gratitude to the Rockefeller Foundation, 
which kindly put at our disposal a large quantity of KCl and of the preparation mentioned 
in the text in order to enable us to continue our researches on oleate systems. 

7) The KCl (+ KOH) concentration at which the damping of the elastic oscillations 
is a minimum lies lower than for the oleate preparation of MERCK. It was recently 
ascertained that it lay even a little lower than the concentration mentioned in the text, 
which lies still close to this minimum. 


®) In the execution of these measurements we were assisted by W. W. H. WEIJZEN 
and W. A. LOEVEN, to whom we also wish to express our thanks here. 


GA 


134 
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obtained in the investigation on the rotation oscillation (see also fig. 8). 
When we compare now the mean No value (= 76 poises) with the mean 
G-1 value (= 134 poises) we still come to the same result as in the pre- 
ceding section: 7 is not equal to, but smaller than G- J, though both are 
of the same order of magnitude. 


Fig. 8. 


We could leave it at this general conclusion, but it is interesting that 
there appears to be a system in the disagreement between the two values. 
See the following survey. (Table VIII), giving the five cases in which 7 
and G-j could be compared, arranged in the order of decreasing 1. 


TABLE VIIL. 

; 1 : Ng (mean) G - 4 (mean) ; 
Concentration Oleate | n. Hexylalc. Mol/L ron etces 9 /G-a 
1.8/9 (MERCK) _ 105 0.40 
1.2% (BAKER) _ 76 0.57 
1.2/9 (MERCK) mae, 60 0.62 
1.2 9/9 (MERCK) 0.0006 50 0.70 
1.2% (MERCK) 0.0018 22 0.85 


We see from column 5 that the disagreement between 7) and G-4 
decreases with a decreasing value of mo. Cf. fig. 9, which gives yo/G- A 
as a function of 7) and which shows that the experimental points fairly 
coincide with a curve drawn through y)/G.4—=1 at y= zero. This 
might mean that MAxwELL’s relation L = G- {A in principle holds good, 
but that certain unknown systematic errors 9) in the experimental methods 
we used for determining G or / or %, distort this equality and the more so 
as the absolute values of these quantities are higher. But it might also mean 
that the equation representing the connection between 7, G and 4 for our 


9) The fear of any occurence of a systematic error in the mean 7 values obtained in 
section 3 has been considerably diminished by the fact that the experimental point obtained 
in this section with a quite different viscosimetric technique lies on the same curve as 
the other experimental points (see fig. 9). 
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systems is really of a more complicated nature and approximates to 
n = G- 2 at low values of these quantities. 

Be that as it may, the same order of magnitude found for G- 4 and 
for 49 seems to indicate that the mechanism of viscous flow in the range 
of small shearing stresses is mainly due to relaxation of elastic stresses. 


Summary. 


1) The viscous and the elastic behaviour of a number of oleate systems 
containing KCl have been investigated at 15° C side by side to compare 
the viscosity coefficient calculated from elastic measurements (y= G- 4, 
MAXWELL) with the viscosity coefficients 7, 7: and , manifesting them- 
selves as levels in the viscosity-shearing stress diagram of the oleate 
systems. 


2) It has been found, that both 7, and 7: are far smaller and of a quite 
different order than G- A, that 79, however, the viscosity coefficient at low 
shearing stresses, is of the same order of magnitude as G- 2d. 


3) The percentual difference between %) and G-d (1 has always 
turned out to be smaller than G-A) decreases in a marked way with 
decreasing absolute value of these quantities. 


4). This suggests that either still unknown systematic errors, which 
increase percentually with increasing values of G, 4 or 7, bring about this 
difference; or that a more complicated relation between G, 4 and yn is 
operative for the oleate fluid, which at small values of these quantities 
may practically be simplified to 7 = G. Ad. 


Department of Medical Chemistry, 
University of Leyden. 


Geology. — Tectonics of the Mt. Aigoual pluton in the southeastern 
Cevennes, France. Part I. By D. DE Waarb. (Communicated by 
Prof. H. A. BRouWER.) 


(Communicated at the meeting of March 26, 1949.) 


1. Introduction. 


In the summer of 1947 detailed tectonic field-work was carried out in 
the southeastern Cevennes of France. The purpose of the investigation was 
to tackle structural problems in the little-known granite massif near the 
Mt. Aigoual. This massif is shown in the southwerstern corner of sheet 
Alais of the French geological map, 1: 80,000, as a red dot of granite 
stretching away to the north in slates. Mapping has been done in and 
around this granite dot (fig. 1) in a purely structural sense. The well- 
exposed country in this part of the Cevennes facilitated detailed study, 
though mapping proved to be hampered sometimes by glacial deposits, 
hillside waste and afforestation. Excellent exposures are provided by a 
large number of small rivers. 

During the survey, the area of the red dot of granite on the geological 
map appeared to comprise a small batholith surrounded by a complicated 
pattern of granitic, porphyritic and dark coloured dikes in the slaty country 
rock. The large scale map at the back shows the geological units of the 
mapped area. 

The author's thanks are due to Professor W. NIEUWENKAMP for helpful 
discussions and criticism in field and laboratory, to Mr. R. C. HEIM for 
valuable data and good companionship during the field-work and to Miss 
D. E. WispEN, M.Sc., Southampton for assistance with the English. 

The writer acknowledges the receipt of a Z.W.O. grant (Dutch Orga- 
nization of Pure Sciences) which made possible some comparative studies 
in adjacent granite massifs of the Cevennes. 


2. Synopsis of literature. 


The Central Plateau is situated in the middle of France as a large isle 
of mainly metamorphic and crystalline rocks with a local cover of young 
volcanic material on its peneplained surface, surrounded by Mesozoic and 
Tertiary sediments (fig. 1). 

The age of the metamorphic rocks has long been uncertain. The 
“Archean age” theory has been rejected since the discovery of Cambrian 
and Ordovician fossils. Precambrian rocks are still mentioned however. 
According to most authors, folding and metamorphism of the Palaeozoic 
geosynclinal strata has taken place in the Sudetic phase of the variscan 
orogeny. The tectonic units with metamorphic imbricate and nappe struc- 
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tures have been analysed especially by the studies of DEMAY (e.g. 1931b, 
1934), summarized by VON GAERTNER (1937 ie 
Also from this part of France granitization phenomena are extensively 
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Fig. 1.. Map of the southeastern part of the Central Massif showing the geographical 
position of the mapped area between slates (white) and adjoining granite massifs (dotted). 
The surrounding mesozoic strata are coarsely dotted. 


described by French authors (e.g. DEMAY, 1935, 1942, RAGUIN, 1930, 1946, 
Roqu_Es, 1941). According to these publications the granite plutons in the 
Central Massif may be distinguished in ‘pretectonic’, ‘syntectonic anatec- 
tic’, ‘synkinetic intrusive’ and ‘posttectonic intrusive’ with respect to the 
Middle Carbonian orogeny. There seems to be no apparent relation 
between anatectic and intrusive granites. No rule could be laid down for 
the mode of occurrence, or for the texture of the posttectonic intrusive 
granites. Sometimes they are porphyritic, but syntectonic plutons too may 
have a similar texture (ROQUES, 1941, RAGUIN, 1946). 

The Mt. Aigoual granite pluton, subject of this paper, is briefly 
mentioned in several French publications (e.g. BERGERON, 1889, 1904, 
Demay, 193la). According to all authors it is a posttectonic intrusion, 
cutting the schistosity of the slates. A narrow metamorphic aureole is 
further evidence of its intrusive character. The texture of the granite is 
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described as homogeneous and porphyritic, with large crystals of orthoclase. 
Because of similarity in these properties, connections are postulated with 
the adjacent St. Guiral pluton. The Mt. Aigoual pluton is mentioned by 
Demay (1931a) as an apophyse of the much larger St. Guiral massif. 

The first mapping of this area was achieved by FABRE and CAYEUX 
(1901) in sheet Alais of the geological map of France and for the second 
edition by THIERY (1923). Of the numerous granitic dikes north and south 
of the granite body only a few have been mapped separately in this sheet. 
Those nearer to the intrusive body are lumped together in one large granite 
outcrop extending southward as far as the St. Guiral massif. 

The morphological study of BAULIG (1928) of the Central Plateau, as 
well as his maps, show many faults and fault-blocks. Two of these faults, 
continuing in one are drawn through the Mt. Aigoual granite. In the 
mapped area, however, little evidence for these faults could be found. 

According to VON GAERTNER (1937) and based on DeMmay (1931a), 
the Mt. Aigoual region forms part of the Orthocevennes complex, a pennine 
nappe-like tectonic unit of the Central Massif. The granite of the Mt. 
Aigoual, together with the southern adjacent massif, is mentioned and 
drawn in section by VON GAERTNER as a late orogenic (posttectonic) 
intrusive sheet, flatlying in the cleavage plane of the slates. 


3. Geology. 

In the centre of the mapped area a small elliptical pluton of granite 
crops out surrounded by slates. The granite body, 5 km long and 2,5 km 
wide is very homogeneous; no endomorphism could be observed. 

The granite is a light coloured, black and white rock of coarse grained 
texture with 10—15 cm long idiomorphic phenocrysts of orthoclase. Main 
components are quartz, orthoclase, plagioclase and biotite. Thanks to the 
granite-porphyritic texture, planar and linear flow structures could be 
measured nearly everywhere in the pluton. Joint systems occur usually in 
a variety of directions; only in few cases have slickensides been formed. 
Aplite, lamprophyre and quartz porphyry dikes have been observed inside 
the granite massif. 

Nearly everywhere round the granite body results of contact meta- 
morphism are visible. The contact plane is always sharp; its strike and 
dip vary locally. 

The country rock farther away from the contact is in general slaty. 
According to DeMAy (1931a) a Cambrian age of the slates would be 
probable. The slates are not subdivided on sheet Alais of the geological 
map of France; on sheet Séverac, some kilometers west of the mapped 
area, a probable Potsdamian age is indicated. 

Normally the slate is veined with white quartz. Though differing in 
quantity, usually numerous thin veins and lenses as well as small dikes 
to 40 cm wide occur. Intensely folded veins have been observed. They 


do not penetrate in joint systems or faults and very seldom in cleavage 


planes. 26 
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In several parts the slaty rock is a well-developed phyllite. Locally 
quartzite layers parallel to the schistosity of the slates are observed; else- 
where microscopical folds in bedded slates cross the cleavage planes. 
Spotted slates too are found in some places. 

Usually the slates possess one or two systems of ribs on their cleavage 
planes. In the south of the mapped area they have minor folds with wave 
lengths between 10 and 100 cm. In this region they often pass into folded 
gneiss. Directions of ribs and fold axes as well as cleavage planes and 
joint systems in the country rock have been measured for structural 
purposes. 

Dikes cut through the country rock round the pluton. Granitic dikes 
in large numbers present themselves as a fringe-like extension of the 
pluton. This edging has a limited extension at the northern and eastern 
side. Dikes are found much farther towards the south, but they are almost 
absent on the western pluton border. 

Lamprophyre dikes occur in groups, mostly outside the edging of granitic 
dikes. They are often intruded along the cleavage planes of the slates. 
Their occurrence within the granite and in porphyry dikes prove them to 
be younger. More basic dike rocks have been observed. 

Quartz porphyry dikes of light coloured material with phenocrysts of 
quartz and also of both quartz and orthoclase are observed between granite 
dikes and inside the granite pluton. They cross granitic dikes and are 
crossed by lamprophyre dikes. 

Faults and faultzones, with and without drag, occur frequently in the 
mapped region. Faults cutting dikes have not been found. In one case 
drag is observed in the country rock at the contact with a granitic dike. 
Quartz veins always are disturbed by faulting; quartz does not occur 
within the fault fissures, only crushed in fault breccia. 

Phe geological history in the mapped region according to the observed 
phenomena may be summarized in the following succession. 

Sedimentation of strata probably in Cambrian time. Folding and meta- 
morphism accompanied by exudation of quartz (according to most authors 
in the Middle Carbonian, Sudetic phase of the variscan orogeny). Develop- 
ment of schistosity and joint systems. Both segregation of quartz and 
jointing will be products of the same cause, the latter somewhat younger 
than the former. Intrusion of granite accompanied by contact meta- 
morphism and intrusion of quartz porphyry, lamprophyre, aplite and 
pegmatite. 


4. Petrology and microtectonics. 


The microscopic compositions and microtectonic features of the mapped 
rocks will be mentioned here in brief. The petrology of the Mt. Aigoual 
area is discussed in a separate paper by HEIM (1949). 

The granite of the pluton is a coarse-grained, light-coloured rock with 
up to 15 cm large twinned automorphic phenocrysts of white K~-feldspar. 
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Main components are white K-feldspar, plagioclase, dark grey grains of 
quartz and many small black biotite crystals dispersed through the rock 
and ina dark rim around the large feldspars. Thin sections show Carlsbad- 
twinned orthoclase, much albite- and pericline-twinned oligoclase, brown 
biotite, quartz and occasionally hornblende. According to the mineralogical 
composition the rock may be classified in the granodioritic subdivision of 
the granites sensu lato. Quartz shows in general slight undulatory 
extinction. 

In the field little difference could be observed between the granite in 
the pluton and in the dikes. The latter may show however some more 
variation in size of the smaller minerals between the large feldspars. Quartz 
too may thus have an automorphic habitus, giving the rock frequently a 
granite-porphyritic texture. In that case the quartz phenocrysts are corroded 
and rounded. Evidence of stronger dynamo metamorphism in the granite 
dikes is shown by wide-spread and strongly undulatory extinction of 
quartz and by bent and frayed biotite crystals. Sometimes recrystallised 
nonundulatory, suture-grained parts within quartz crystals do occur. The 
microscopic texture as a rule is lacking in orientation. 

The quartz porphyry dikes contain a light-coloured, dense rock with 
many automorphic, dark-coloured quartz grains and irregularly distributed 
feldspar phenocrysts. The dikes are sharply jointed in small regular blocks. 
The rock is composed of quartz, orthoclase, biotite and albite phenocrysts 
in a dense, partly spherulitic or micrographic groundmass of mainly quartz 
and orthoclase. The quartz phenocrysts are largely corroded and rounded 
and sometimes show undulatory extinction. 

Dikes of a dark quartz porphyry have been mapped within the granite 
body, near the contact of the pluton, cutting through granitic dikes and 
forming bordering zones of the latter. It is a dark rock of biotite and 
amphibole with quartz and plagioclase phenocrysts and a few K-feldspar 
phenocrysts. Under the microscope this rock proved to be a tonalite 
porphyry. Large rounded, corroded and undulatory quartz phenocrysts, 
zonal, Carlsbad-, albite- and pericline-twinned, corroded andesine, brown 
biotite, hornblende and occasionally orthoclase phenocrysts are enclosed 
in a fine-grained groundmass, often with fluidal texture, of andesine, 
biotite, hornblende and some quartz. 

As mentioned above different types of basic rocks have been observed. 
They show difference in darkness, in quantity of feldspar and in size and 
orientation of biotite. Usually they are dark coloured, fine-grained, glitter- 
ing rocks, largely composed of biotite and feldspar. In dike outcrops they 
often show rounded blocks caused by weathering along joints or they 
may be weathered all through into yellow-brown, spotted, sandy material. 
In some dikes the rock has an oriented texture, due to relatively large 
biotite crystals, being arranged parallel to the dike contacts. 

Most dikes have kersantite-like composition; one consists of greenish- 


grey porphyrite. 
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The aplite dikes normally contain quartz, orthoclase and some biotite; 
no exceptional minerals have been found in the aplites, nor in the few 
pegmatite or runite dikes and veins. 

The country rock is usually microfolded slate with quartz veins, phyllitic 
with sericite and recrystallised quartz, blastophyllitic with garnet and 
kyanite, and schistose to almost gneissic with sericite, muscovite, biotite, 
chlorite, quartz and some albite. In the narrow contact zone a fine grained 
gneiss is locally developed with orthoclase and plagioclase. Andalusite, 
sillimanite and cordierite also have been observed near the contact. 


5. Granite contacts. 


Many exposures allowed a detailed observation of the pluton borders. 
The western contact forms a nearly straight line in the plane of the 
structure map. Slates border the granite along a sharply defined contact 
plane which dips, at an angle of 70° to 55°, under the slate. This contact 
has all the aspects of a set of joints, a pre-~ or syn-intrusive joint system 
without doubt, because of the existing contact metamorphism in the slate. 
A narrow zone of slightly metamorphosed slates and phyllites borders the 
contact with penetrating veins of aplite and pegmatite. No offshoots such 
as granitic dikes are to be found on this side of the pluton. 

Joint-like contact planes have been observed nearly everywhere around 
the pluton, The northern border is complicated by the branching off of 
many dikes. These offshoots, splitting apart euhedral blocks of slate, 
could be observed in detail. It has not been feasible however to reproduce 
these complications in the structure map; the curved lines merely represent 
a general outline of the granite body in these places. The northern contact 
planes show great variety of orientation; steep discordant and low angled 
often concordant contacts have been measured. Blocks of slates, surrounded 
by granite dikes could be observed. 

Contact metamorphism in the north and northeast of the pluton may 
locally be stronger than near the western border, though more than one 
meter of the country rock has usually not been changed megascopically. 
The slates are transformed in these places in a fine-grained gneiss often 
veined or brecciated with aplite and occasionally with pegmatite. The 
slates and metamorphosed country rock near the contact show variations 
in strike and dip. Dips proved to be usually steeper and strikes are 
systematically bent to northeasterly direction. 

Parts of the eastern and southeastern contacts are hidden by hillside 
waste. Where visible the same joint-like contact planes with shallow 
metamorphism, acid veins and offshoots of granitic dikes occur. Here too, 
strikes and dips of the slates change near the contact. The contact plane 
usually dips under the slates with varying angles. 


On the map there are two gaps in the southern border. In the eastern 
one the contact loses itself under debris. In the western gap there has 
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been no occasion for mapping further south; it seems quite possible 
however that across this track a bottleneck may join up the Mt. Aigoual 
pluton with the St. Guiral massif. The southern contact too has little 
contact metamorphism. Gneiss and acid veins are observed; mostly there 
is hardly any change beyond a hardening of the country rock. 

Though the contact may be clean and sharp as in the western border, 
complicated contacts have been found. In a river exposure of the northern 
contact the following details have been observed. The granite near 
the contact, being quite normal with dark inclusions and at most a little 
less large orthoclase crystals is bordered by a small tonalite porphyry dike 
of 20 cm followed by a quartz porphyry dike of 30 om thickness with 
inclusions. Then, along the contact there is a kind of conglomerate zone, 
25 cm thick with rounded pieces of granite and metamorphosed slate in a 
fine-grained crystalline matrix. On the outer side of the contact and 
parallel to it a second tonalite porphyry dike cuts the slightly metamor- 
phosed and hardened slate. In the south an exposure is found with the 
succession of normal granite, a 3 m zone of somewhat darker granite, 
a 50 cm zone of granite, normal in colour but more crowded with large 
orthoclase crystals than usual, a 75 cm dike of tonalite porphyry in contact 
with hardened slate of which a slice of 5 m is cut off by a quartz porphyry 
dike running parallel with the contact. 

Still less of metamorphism has been found near the contacts of the 
granite or granite porphyry dikes, though local differences occur. 

The contacts of the granite massif, usually cross-cutting and occasionally 
concordant with the slaty cleavage and the slight but never missing 
contact metamorphism points undoubtly to an intrusive character of the 
granite. As a whole the pluton seems to be bordered by joint systems. 
They must have been developed before or during the intrusion. During 
the intrusion movement along these joint systems has taken place. The 
intrusion has been accompanied by an uplift of the country rock. This is 
proved by strike and dip alterations near the contact. Post-intrusive move- 
ments along contact planes cannot have been of much importance because 
of the uninterrupted zone of contact metamorphism. 

The contact zone described above with rounded granite and gneissic 
pebbles in a fine-grained crystalline mass, points to intrusive movements 
along this plane by which parts of both sides of the contact were crushed, 
rounded and embedded in a fluidal mass, causing an intrusion bordering 
mylonite. 


6. Outline of tectonic phenomena. 

The main point of the investigations has been to collect tectonic data 
in order to work out the movements and origin of the granite pluton. 
A summary is given in the structure map at the back. The detailed mapping 
in a relatively small area, with heights between 800 and 1400 meters, made 
it desirable to ‘have all data transposed in one horizontal plane of 1000 
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meter. Topographic distortion is eliminated in this way. The structure map 
shows as it were a hypothetical peneplain at 1000 m in this area. Most 
measurements have been made near this 1000 m plane. 

A fairly close network of tectonic data has been collected in the granite 
massif. Special attention was given to fluidal phenomena and fracture 
systems in order to be able to reconstruct the intrusive movements. Slicken- 
sides and directions of dikes were likewise of much importance with regards 
to these movements. Marginal observations on contact planes and inclusions 
complete the network of structural data. 

In the country rock many important tectonic features could be measured. 
Strike and dip of cleavage planes of the slates are disturbed near the 
contact as mentioned in the preceding paragraph. Joint systems, ribs and 
minor folds in the slates give details of pre-intrusive tectonics. Dike measure- 
ments will possibly reveal relations between pre-intrusive tectonics and 
intrusive movements. Faults with drag phenomena may be connected with 
block. movements during the intrusion. 

These tectonic data, plotted in equiareal diagrams, will be discussed 
systematically in the following paragraphs. 


7. Flow structures. 


Flow structures in the Mt. Aigoual granite are mainly marked by 
oriented orthoclase crystals and, in its marginal parts, also by small 
inclusions. 

The orthoclase crystals are mostly tabular after (010) and elongated 
in the direction of the c-axis. Their orientation may be linear with 
parallelism of the c-axes or planar with parallelism of the (010) planes. 
Usually they are twins after the Carlsbad law, showing interrupted (001) 
cleavage on their narrow sides. 

Primary flow phenomena in crystalline rocks may be divided into linear 
and planar structures. Ordinarily one of the two is predominant or exclusive 
in a massif. Only when there are suitable minerals can both these flow 
structures be seen. Feldspar is one of these, because of its elongated 
tabular shape. In the Mt. Aigoual granite planar structures are favoured 
for its feldspars are distinctly more tabular than columnar in shape. 

In spite of the possibility of recognizing flow structures by measuring 
the directions of the phenocrysts, their orientation at a first glance seems 
to be chaotic in many places. Only a careful investigation of several 
surfaces of different orientation of the granite reveals a predominant 
arrangement, A judicious choice of the order in which the surfaces are 
investigated, e.g. horizontal planes first, followed by surfaces at right 
angles to it assists in deciding between linear and planar orientation. In the 
latter case a surface parallel to this plane of flow was selected in order to 
reveal any tendency to a linear arrangement within this plane. As a rule 
only after a prolonged scrutiny and with many measurements of generally 
existing, local directions on many different planes of the granite, could 
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the predominant orientation be detected. Thus the field observations do 
not consist of measurements of single feldspar crystals which could only 
give reliable information if a large number (e.g. 100) were treated at each 
locality, but of measuring predominant directions on fairly large surfaces 
of granite outcrops. 

Parallelism of minerals is caused by differential motion in the intrusive 
mass. Consolidation preserved a record of the ultimate movements. Linear 
parallelism is originated by preponderant one-dimensional elongation of 
the mass. Stretching in one direction results in a parallel orientation of 
the longest axes, the c-axis of the feldspars, in that direction. Their (010) 
planes have arbitrary positions rotated about the c-axes (fig. 2b). {If 
differential motion in the mass effects equal stretching in all directions 
within (shortening normal to) a plane, planar parallelism will develop. 


ro 
PLANAR AINEAR 


Fig. 2. Ideal arrangement of feldspar phenocrysts in primary oriented granite. & (010) 

planes of feldspars in planar arrarigement. b. c-axes of feldspar crystals in linear 

orientation. c. combination of linear flow structure by parallelism of c-axes and planar 
flow structure by platty arrangement of (010) planes. 
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Feldspars will tend to have (010) parallel to this plane, while their 
elongated c-axes show haphazard orientation (fig. 2a). 

In the intermediate case, when planar structures are developed together 
with somewhat stronger stretching in a certain direction within this plane, 
both (010) planes are oriented in planar parallelism and c-axes in linear 
parallelism in the same plane (fig. 2c). Only locally are linear and planar 
orientation found associated; usually one of the two is predominant and 
the other weak or absent. 

Because differential motion is strongest near the contacts of intrusions 
flow structures will be developed very well in their marginal parts. Con- 
versely strongly oriented structures in the central part of the intrusion 
indicate that considerable differential movements also in these parts of 
the mass have existed. Usually these structures are planar. If linear 
parallelism exists too, strong one-dimensional stretching has been predo- 
minant. In parts of the intrusion where no planar structures are developed 
a weak linear orientation may be formed by distinct expansion of the 
pluton in a certain direction. 

The movements of flow are very well determined if both planar and 
linear flow structures could be measured. Planar parallelism shows the 
plane along which flow has taken place and linear parallelism the direction 
of flow in that plane. 


8. Intrusive flow. 


The general aspect of the registered and interpolated flow data in the 
structure map of the pluton suggests comparison with a pan of boiling 
porridge. In at least three sections of this relatively small pluton upward 
culmination of flow is recorded. 

This reconstruction is based on about 65 flow determinations. Con- 
struction of flow planes and lines is done by interpolation of the measured 
data in the least complicated way. These lines record in a nice manner 
the last movements of the intrusive mass before ultimate consolidation. 
The given interpolation is a reconstruction in the simplest way; in reality 
the flow structures may be much more complicated in detail. 

West of the centre of the pluton a half concentric doming is made 
visible by flow phenomena. It has varying but mostly rather steep flow 
structures, planes as well as lines. The planar structure has conformable 
contacts with the westerly border, dips to the south against another, 
incompletely-known doming and continues to the E.N.E. To the north, 
structures are more complicated but presumably still conformable with 
the pluton contacts, continuing into dike offshoots. Conformable structures 
are also well-developed in granite dikes. Linear parallelism in north and 
west cross flow plane strikes nearly perpendicularly, indicating probably 
a strong upward flow. Near the top of the doming the lineation grows 
parallel with planar structures. 

In the north-east of the pluton the continuation or easterly flank of the 
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doming is found. As far as observed it has identical conformable structures. 
Flow lines here may show stronger convergence in upward flow, also from 
northerly directions. 

Flow structures in the south-east are different. They seem to be con-~ 
formable to parts of the southern contact but disconformable to the easterly 
wall rock, representing a mass with parallel oriented planar structures, 
complicated by vertical and changing dip directions. They suggest a lacking 
or hidden part of the mass still S.E. of the southern pluton border. 

These structures might have been caused by motion of a mass ascending 
from south-easterly direction; a northward motion which produces at the 
same time a drawing-out of the northern doming or elongation of the 
northern part of the pluton in S.W.—N.E. direction. Such an ‘elongation 
may have originated the flow lines in that direction, which are bent in the 
central part of the pluton. 

In the south-western part of the pluton a fragmentarily-known small 
doming is visible. The same parallel and steep flow lines as in the south- 
eastern section have been found in the south-western appendix or bottle- 
neck, 

All together the flow structures in the pluton suggest an upward moving 
mass, converging from all directions, which was strongly affected by the 
ascending mass from the south, causing the elongation and bending out 
of the northern doming. 


9. Fracturing of granite. 

In the pluton joint systems are found in amazing quantity. Normally 
the Mt. Aigoual granite is fractured by about six different systems of 
joints in each locality. Up to twelve systems have been measured in some 
places. By fracturing the granite falls apart in different sized and irregular 
blocks, frequently parallelepipe and column shaped. In fresh rock, joint 
planes are hardly visible; in a few cases faint marks of dislocation are 
observed. The explanation of the existence of fracturing may be contraction 
during cooling of the granite, or continued motion within the pluton after 
consolidation of the upper shell, or a combination of both causes. In either 
case, however, the resulting fracture systems cannot be quite arbitrary, 
there certainly will be a relation with the fabric of the rock and the shape 
of the massif. 

The connection of flow structures and some joint systems is advocated 
by Hans Cioos, BaLk and others. They make a distinction between e.g. 
cross joints perpendicular to flow lines and longitudinal joints parallel to 
them. Many joint systems could not be placed systematically by these 
authors. 

In the case in question it proved always to be possible to find some of 
the six or twelve joint systems in regular association with flow structures. 
But there is nothing specific which marks the thus found ‘cross’ or 
“longitudinal” joints; apart from their relations to flow structures nothing 
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may distinguish them from the other fractures. Most joints however cannot 
be explained with reference to these directions of use; they seem to be 
useless from tectonic point of view. In this case it would be a rather 
arbitrary application, to select and use “cross” joints, etc., according to 
recipe. 

In the method here applied all measured joint data are used in a 
statistical diagram. This is done in order to bring out clearly the possible 
regularities of jointing in the pluton as a whole. Though it has proved, 
so far, not to be possible to arrange all existing joint systems according to 
structural or other principles, this of course in no way excludes the 
possibility of a systematic arrangement. There may be suggested here a 
systematic coincidence between the fracturing of the rock and its mechanical 
properties in different directions. Any alteration in equilibrium of the mass, 
e.g. during cooling, contraction, renewed motion etc., will thus cause 
systematic joint systems. They cannot all be explained locally, but they 
will relate as a whole to the fabric of rock and massif. 

Measurements of joint systems have been made at regular intervals 
throughout the pluton outcrop..An average diagram of fracturing of the 
pluton may reveal special regularities in the structure of the massif. 
Projections of poles of all joint planes are plotted in an equal area net, 
thus giving a petrofabric analysis of the whole pluton as far as its joints 
are concerned. Fig. 3 shows the statistical pole diagram of 80 measured 
joint planes. 


Fig. 3. Pole diagram of 80 joints throughout the pluton. Contours 5—4—24 % 
(Southern hemisphere as in fig. 5, 6, 7, 9 and 10.) 


The first striking regularity is the evident symmetry of the diagram. 
‘The N 120 E. central line of the diagram shows the same density distribu- 
tion at both sides. The direction of this symmetry plane harmonizes with 
shape and directions of flow structures of the pluton except for the south- 
western corner. Probably this symmetry in the diagram is connected with 
the predominant motion in the direction of the line of symmetry and with 


the mentioned probability of supply of granite mass from a south-easterly 
direction. ; 
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By further analysis of the diagram two girdles may be observed. They 
resemble B-tectonites or a cross girdle diagram in petrofabrics, one NW— 
SE and the other horizontal, crossing each other at right angles. This 
indicates an orientation of nearly all measured joints either parallel to a 
horizontal NE—SW axis or in a more or less vertical arrangement. 

Because joints largely will be related to, or oriented as a result of the 
fabric of the rock, in this case the only visible variable viz. flow structures, 
this diagram may show an average orientation of joints associated with the 
average orientation of flow. A glance at the structure map shows a predo- 
minance of flow lines in NE—SW direction which may be associated here 
with domination of cross joints in the densest parts of the horizontal girdle 
viz. in the NE and SW of the diagram. 

In the same way a relation does not seem accidental between the densest 
part of the vertical girdle in the middle and south-east of the middle of 
the diagram with the ascent of mass from a south-easterly direction. This 
densest part of the diagram indicates the nearly everywhere observed 
“bedding planes” or primary flat-lying joints throughout the pluton. 


10. Dikes and faults in the pluton. 


Besides numerous ordinary fractures described in the preceding para- 
graph, dikes and faults have been observed in the pluton. In this paragraph 
joints filled with aplite, lamprophyre and porphyry and faults with well- 
developed slickensides or with visible dislocation will be discussed. More 
than fractures they may be expected to furnish information about motion 
in the rigid or semirigid upper part of the pluton. These data are plotted 
in the structure map. 

The discussion of these phenomena is hampered however, by their limited 
number. The pluton body proved to be relatively poor in aplite dikes. 
Three of them may be called of cross joint origin, two are along primary 
flat-lying joints, two others show affinity to longitudinal joints and the 
few remaining cannot be labelled according to flow structures and shape 
of the massif. It will be clear that use of this nomenclature has no practical 
value here for tectonic purposes because of the small quantity of data. 
The few measurements of lamprophyre and porphyry dikes cannot be © 
used either. 

Faults usually with striae could be observed clearly in the north-western 
sector of the pluton. They are plotted in tectonogram fig. 4. The diagram 
is composed of equiangular projections of: a fault zone striking N 37 E, 
dipping steeply W;; the sharp bisectrix N 50 E, 7 S of a system of diagonal 
joints or faults, crossing each other nearly perpendicular (84°); horizontal 
striae in a fault zone striking N 45 E; and three sets of striae on different 
fault planes in one locality, striking in a north-western direction and 
dipping S. The last three fault systems with striae are likely to represent 
the result of the same direction of motion recorded on different planes 
in this locality; combined they may give the direction of the original 
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motion, viz. S (1.2.3.) in the diagram, striking N 38 E and dipping 25 SW. 
These plotted data show remarkable conformity; they have uniform 
directions between N 35 E and N 50-E, with gentle dips, mostly to the 


Fig. 4. Tectonogram of faulting and slickensiding in the north-eastern sector of the 
pluton, Equiangular projections of striae (S), the sharp bissectrix of diagonal joints (B) 
and a fault zone without striae (F). (Projections in the northern hemisphere.) 


SW. Thus, motion is recorded in the rigid upper shell of the granite in 
SW—NE direction, slightly upward to the north-east. 

In an earlier phase of the consolidation, motion in the same direction 
has taken place in the area north-west of the centre of the pluton. Two 
aplite dikes are several times broken and irregularly dislocated by sets of 
parallel steep faults, striking about N 40 E. Bent parts of the aplite dikes 
and fixing of the fault planes indicate a semirigid condition of the granite 
in which SW—NE motion has occurred. 

The SW—NE motion thus found in the northern part of the massif 
may be interpreted as a result of differential expansion of the pluton in 
that direction during the last phases of the intrusion. Similar indications 
have been found by means of flow structures in earlier intrusive phases. 
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Anatomy. — The digital formula in relation to age, sex and constitutional 


type. I. By J. Huizinca. (Communicated by Prof. M. W. 
'W OERDEMAN. ) 


(Communicated at the meeting of March 26, 1949.) 


I. Introduction. 


If we ask various persons to lay the pronated hand on a flat surface 
in such a way that the longitudinal axis of the hand is a prolongation of 
that of the forearm we shall find that: 


The third finger ends the most distally. 
The first finger ends the least distally. 
The fifth finger follows the first in this respect. 
4. Sometimes the second finger ends more distally than the fourth, 
sometimes the converse is true and sometimes these fingers are equal 


in length. 


oy 


We may represent a given case as follows: 
TSiV > iS V1 


This Woop Jones (1944) terms the ‘digital formula’. 

If we confine ourselves to man, we find that the interindividual difference 
in digital formula consists in the varying relation between the fingers 
II and IV. 

Although it is far from our intention to recommend memoirs such as 
those of CASANOVA as a source of scientific information, we feel justified 
in quoting CASANOVA’s description of his conversation with the painter 
RAFAEL MENGS by way of introduction to our problem, the more so as 
it is our experience that the same dispute about the digital formula can be 
provoked in any group of people at the present day. 

CASANOVA writes as follows (1871): 


“Je me souviens qu'un jour je pris la liberté de lui faire observer, en 
voyant un de ses tableaux, que la main d’une certaine figure me paraissait 
manquée. En effet, le quatriéme doigt était plus court que le second. 

— Voila une plaisante observation, me dit-il, voyez ma main! et il 
l’étendit. 

— Voyez la mienne, répondis je, je suis convaincu qu'elle ne différe 
pas de celle des autres enfants d’Adam. 

— De qui donc me faites vous descendre? répliqua-t-il. 

— Ma foi! lui dis-je aprés avoir examiné sa dextre, je ne sais a quelle 
espéce vous rattachez, mais vous n’appartenez pas 4 la mienne. 

— Alors votre espéce n’est pas I‘humaine, car la forme manuelle de 
l'homme et de la femme est bien celle que voila. 
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— Je parie 100 pistoles que vous vous trompez, lui dis-je. Furieux 
de mon défi, il jette palette et pinceaux, sonne ses gens, et leur fait a 
tous exhiber leurs mains; sa colére fut grandi quand il reconnut que chez 
tous le doigt annulaire était plus long que I’index. Cependant il voulut 
bien sentir le ridicule de sa conduite et termina la scene par cette 
plaisanterie: 

— Je suis charmé du moins d’étre unique en mon genre sur un certain 
point.” 


In 1875 EckER, was the first anatomist-anthropologist to bring up the 
problem of the individually-differing prominence of the fingers. Since 
then numerous publications have appeared; these may be. classified as 
follows on the basis of certain general principles: 


Prominence and sexual dimorphism 

Prominence and age 

Prominence in the light of typology 

Prominence in the light of racial differences 
Prominence in the light of problems of evolution 
Prominence differences between the two hands. 
Discussion of the causation of differences in prominence. 


SR Bee rae 


No single author, however, has dealt with all these aspects at once; 
this is partly due to the lack of due insight into anthropologo-phenomeno- 
logical problems. In some cases the contradictory nature of the statements 
made can be ascribed to the fact that the groups studied were non- 
comparable. For instance, data furnished by the study of a group of 
females aged 4 to 71 years, in the absence of any previous investigation 
of differences according to age, may well lead to premature conclusions 
about the phenomenon in women, 

Although the population of the Netherlands can certainly not be 
regarded as racially homogenous (in addition to definite Nordic, Alpine 
and Mediterranean characteristics we also find Baltic and Dinaric 
features), analysis of our material from this point of view is vitiated by 
so many uncertainties of race-diagnosis that its results are not worth 
reporting. 

The amount of anthropoid material available was so small that we do 
not feel justified in including it in our study. 

To give some idea of the scope of the problem of prominence, a discussion 
of the points 1 to 7 mentioned above will precede the description of our 
own observations. 

WOLOTZKOI (1924) devised a useful nomenclature and we shall follow 
him in speaking of hands of the radial type (Rd.) when the index finger 
extends more distally than the ring finger, and of hands of the ulnar type 
(UlIn.) when the converse holds. The results of our own investigations 
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lead us to use the term transitional type (T.) for hands in which the index 
and ring fingers extend equally far distally. This nomenclature will also 
be used in discussion of the work of other authors. 


II. Survey of the literature. 


1. Prominence and sexual dimorphism. 


ECKER, who ‘was the first anthropologist to publish an article on this 
‘oscillating character in the hand of men’ (1875), found 24 examples of 
the ulnar type and one of the T. type in 25 outlines of the hands of 
American negroes aged 19 to 65 years. In his group of negresses (age 4 
to 71 years) he found 15 Uln., 6 Rd. and one T. 

With a total absence of criticism as to the age-composition of his groups 
(note the women) he concludes that there is an unmistakable sex 
difference: Rd. occurs more in women than in men. From a group of 
Europeans (composition unknown) he drew the same conclusions, 
although with some reservations. 

The data reported by MANTEGAZZA (1877) (ages not stated) make it 
possible to calculate 75 % Uln. for 258 Italian women and 92 % Uln. for 
336 Italian men. In these groups of individuals examined by him, only Rd, 
or Uln. is found. Thus the Uln. type predominates in both sexes, as found 
by EcCKER (1875). 

PFITZNER (1893) found just the contrary for skeleton hands of adult 
Alsatians: 70 % of 175 male hands were of the Rd. type and 79 % of 90 
female hands, while RUGGLES (1930) and BAKER (1888) found in white 
Americans that more Uln. types on the whole occurred among men and 
more Rd. types among women. 

WEISSENBERG (1895) is more inclined to agree with MANTEGAZZA; he 
found more Rd. in women than in men but the the percentage of Rd. was 
always below 50 %, so that Uln. predominated in both sexes. 

The next publication on sex differences in the digital formula did not 
appear until 1924 (WoLOTZKO!I): in adult Russians Rd. was found in 
62 % of 190 men and 77 % of 159 women. For adult Russian Jews the 
figures are 59 % Rd. in 29 men and 62 % Rd. in 58 women. These figures 
agree with those of PFITZNER (1893). 

WOLOTZKOI then draws the inaccurate conclusion ‘that the radial form 
is a special property of the female hand’. 

RUGGLES (1930) concluded from a study of 402 male and 218 female 
‘white adults’ that ‘In white adults the ring finger in males is generally 
longer than the index finger and in females the reverse is found’, a 
conclusion which he (wrongly) believes to be identical with those of 
PFITZNER, ECKER and MANTEGAZZA and contradictory to those of SCHULTZ 
(1924) and Woop Jones, neither of whom, however, (the latter at any 
rate not in his book published in 1944) makes any mention of sexual 
dimorphism in a comparative sense. As we have not been able to obtain 
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a copy of the publication of BAKER (1888), we are doubtful as to what 
RUGGLES describes as BAKER’s results (see above). 

WECHSLER (1939) found more Uln. types in men than in women. 

A number of the older anatomists (e.g. GEGENBAUER (1885), KOLLMANN 
(1886) believed Rd. to occur more in women than in men, although they 
did not give quantitative expression to this. The more frequent occurrence 
of Rd. in women together with the greater beauty of form (in the opinion 
of many) of the female hand, led various investigators to study the way 
in which artists depict the hands of their models. ECKER (1875) makes the 
following pronouncement: ... ‘wherever a great artist has endeavoured, 
whether instinctively or consciously, to depict a hand of perfect beauty ...... 
he certainly never makes the index finger appreciably shorter than the ring 
finger as this formation definitely gives the stamp of a lower type’. 

WEISSENBERG (1895) did not confirm this in his study of Egyptian 
and Assyrian art. However, we shall confine ourselves to the mere outline 
of this aspect of the problem. Summing-up we may remark. that the 
literature fails to provide us with unequivocal information on the sex 
differences in the digital formula. We shall return to this question in 
connection with our own investigations. 


2. Prominence and age. 


Much less has been written about the connection between relative 
length of fingers and age than about the difference between the sexes 
in this respect. We have already seen how ECKER (1875) put females 
aged 4 to 71 years in a single group labelled ‘women’ and then came to 
the conclusion that there was an unmistakable difference between the sexes. 

WEISSENBERG (1895) classified his 574 male Jews according to age 
as well as sex. For the right hand he gives: 


45 .50/o 57 18 .90/g 29 25.0% 
Uln. 34 51.5%) | 222 73.8% 86 70.2% 
x 2 3. 00/p 22 7 .30/g 8 4.80 


From this it follows that boys from 5 to 10 years of age show the Rd type 
more frequently (45 %) than older boys (25 % for age about 20 yr.). 
Although the proportion of Rd types is higher between the ages of 5 and 
10, WEISSENBERG’s data show it to remain still below 50 %. He also 
remarks ‘both types of hand may be found even in new-born infants’. 

WOLOTZKO!I (1924) arranged his Russian and Jewish men and women 
in age-groups as proposed by STRATZ (1903). In the periods from 1 to 4 
and 8 to 10 years growth in breadth is regarded as predominating over 
that in height (first and second filling-out periods; turgor primus et turgor 
secundus), while those of 5 to 7 and 11 to 14 years correspond to relatively 
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greater increase in height (first and second periods of extension; proceritas 
prima et secunda). Then follows the maturation period from 15 to 20 years. 

In order to facilitate comparison with our own findings we give 
WOLOTZKO!'s figures in full: 


Russians. 


aS Males | ___s Females 

7 : 1 Se 

number | °/9/gRd. | %/0%/o Uln. | %/o%/o T . number |9/00/o Rd.| %/00/o Uln. | %/o% T 
is | 13 | 7 21 | 2 | 18 | 60 2% | 12 
5-7 | 399 | 81 5 | 4 67 | 76 a1) 3 
S10 1 ss | oe] ho Cs 53 77 2) OF a2 
11—14 | a oe ee? ae 184 63 29 8 
5-20 | 61 | 53 44 3 52. | 67 33 = 
21—older | 190 . 62 344 0CUd| C4 159 77 ll eee 


a, A Males Females 
| number 0/,9/gRd.| 9/6%/g Un. | %/99%g T | number | /o%/o Rd. 9/99/9 Uln. | 9/9%o T 
s—7 | 9 99 | 10 = ie 
8-10 | 10 7O | 30 ~ | 616 
11-14 | 10 20 80 ba | 8 
15—20 13 54 38 | 8 | 20 80 20 | as 
2l—older 29 | #59 ml wp ss | 62 Le hee 


WOLOTZzKOI concludes that the hands of children show a predominance 
of the Rd. type, but that with increasing age the number of Uln. forms 
increases. He also remarks that after the 20th year the converse 
phenomenon appears and the number of Uln. decreases in favour of Rd. 
forms. (This is not the case with the Jewesses, J. H.). 

Although both WEISSENBERG and WOLOTZKO! examined Russian Jews, 
there is an enormous difference between their percentages for the different 
hand-types. WEISSENBERG finds 70% of Uln. forms in adult males, 
whereas WOLOTZKOI gives 27 %. To what extent the cause of these 
discrepancies is to be sought in the number of subjects examined, we need 
not discuss here. It seems justifiable to conclude from both these 
investigations that the Rd. type is more frequent in youth than at a more 
advanced age, although the actual percentage reported by these two 
authors differ widely. 

We shall discuss this conclusion further in connection with our own 


work. 


As far as we know the only publication in which mention is made of the 
hands during intrauterine life is that of MIERZECKI (1946). Without 
specifying the number of subjects examined or their sex, he gives the 


following percentages for negroes (N) and whites (W) 
27 
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w | N w | N wo TUN 


| 3th. month | 4th. month | __ 9th. month 10th. month 
33,5 | 63.2 )) 38 54.6 


w | N 


Uln. 3353 57.1 ay 50.0 
ls (Kes // 427) 69.0 50.0 Dona 36.8 5450 45.4 
Rd, 2.0 = 16.3 — 1321 — 13.6 — 


In whites, thus, it appears that the index finger is longer than the ring 
finger in the third month in only 2% of cases, whereas at birth the 
occurrence of this relation is 14 %. According to MIERZECKI, in negroes 
the hand during intrauterine life is invariably non-Rd. ‘In the development 
of the hand of the negro, as in that of apes, no tendency whatever is seen 
for the length relation between the fingers to change in favour of the 
index finger.’ 

The high percentage of T. in the above table is remarkable. The number 
of cases in which the investigator is unable to make a decision is largely 
dependent on the method of examination used. What looks like a T form 
on simple inspection may be shown by accurate measurements to be Rd. 
or Uln., or vice versa. We may also agree to use the terms Rd. or Uln. 
only in cases where the difference in prominence exceeds a given number 
of millimetres. MIERZECKI is silent on this point. 

In addition to this, the uncritical fashion in which he makes use of 
EckERr’s data makes it impossible for us to have very great faith in his 
percentages, 


3. Prominence in the light of typology. 


In 1875 (long before the time of Kretschmer) ECKER remarked in his 
publication that, where the Rd type occurs in European men, this ‘is found 
more frequently in tall, thin individuals than in those of short, stocky 
build’.. We understand him to mean that this does not hold for women. 

The Viennese investigator ROMICH (1932) describes the distribution 
of the two hand forms among the constitutional types distinguished by him. 

The progressive constitutional type (P.T.) ‘includes the sum of all 
progressive characteristics that are necessary for the static functions and 
that find expression in the transformations of the entire locomotor 
apparatus which are brought about by and adapted to these functions’. 
This part of his definition will, we believe, suffice to show the direction 
in which the characteristics are oriented. 

The conservative constitutional type (C. T.) ‘adapts itself, with respect 
to the locomotor apparatus, for the dynamic function and shows a 
pronounced accentuation of the rudimentary formation’. 

Among 300 adults of both sexes (proportions not stated, J. H.) 
consisting of 150 P.T. and 150 C.T. individuals, 40 % showed Uln., 
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51% Rd. and 9% T. forms. Classifying these according to C. T. and 
P. T. types, however, ROMICH found: 


——————————— 


| Pols Cyr; 


Uln. 609/o 209/o 
2 80%/o 10°%/o 
Rd. 32/0 700/o 


The percentage of Rd. forms for the whole group (51 %) appears thus 
to result from the occurrence of Rd. in 70 % of the conservative and 32 % 
of the progressive groups (these groups being numerically equal). RoMICH 
further states that in his progressive type (ie. that in which the static 
function is to the fore) a narrow, gracile hand with long fingers is usually 
found, while the hand of the conservative type is short and broad with 
short fingers. 

This would mean that the progressive hand corresponds to that of the 
leptosome type, in which, therefore, the UlIn. form must predominate 
(60 %). If we are really justified in using the term leptosome here, it 
seems that this conclusion conflicts with the findings of ECKER already 
mentioned. 

We are, further, inclined to believe that ROMICH’s typology is no great 
acquisition. 


4. Prominence and racial differences. 

We have already mentioned the distribution of the two hand forms 
found by EckER (1875) in American negroes. Although ECKER makes no 
definite statement, SCHAAFFHAUSEN (1884) after comparing the-‘so-called 
savages with civilized human beings’ considers himself justified in 
remarking that it is ‘as ECKER was the first to show, a characteristic of 
culture versus savagery that the index finger increases in length relative 
to the fourth or ring finger’. Although we ourselves have not investigated 
racial differences, we do not wish to cast doubt on the possibility that 
differences may exist, even to a very large extent. But we see in the fact 
that judgements on this problem, pronounced at the end of the 19th 
century, were frequently based on the examination of 2, 3, 4 or 5, 
individuals of little-known races — and often without any attention to 
age or sex — a reason for not attaching undue value to such statements. 

ViRCHOW (as stated by WECHSLER in 1939) asserted in 1898 that 
cultured peoples have the Rd. type while primitive peoples show the less 
elegantly proportioned hand of the Uln. type. ViRCHOW also remarked 
(in our opinion humorously): ‘the tendency to a longer index finger 
happened to be greatest in the chief of the negro tribe studied: in him 
there was no difference between these two fingers’. 

We have already mentioned the results of PFITZNER (1893) with 
Alsatians and those of WoLOTZKOI (1924) with Russians and Jews, as 
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well as those of WEISSENBERG (1895) with Jews. According to SCHULTZ 
(1924) the rule that primates (see below) have the Uln. type holds for 
negroes, 

RUGGLES (1930) states in connection with his American whites: ‘There 
is no indication of any relationship betwee the finger type and (either 
handedness or) eye color’. 

Comparison of the results of racial investigations gives a confused picture. 


5. Prominence and problems of evolution. 

Here again we must mention ECKER who was the first to include 
anthropoids in his investigations. He remarks that in all the apes examined 
— but least in the Gorilla — the Uln. type occurs. The number of apes 
examined was very small. 

SCHAAFFHAUSEN (1884), the student of ‘savages’ already mentioned, 
states in connection with the occurrence of the Rd. type ‘this is seen in 
none of the anthropoid animals; in these the ring finger is invariably the 
longer and the index finger the shorter’. 

HARTMANN (1883) whom WEISSENBERG calls one of the greatest experts 

on the anthropoids, states, however, that Rd. prominence does occur in 
apes. : 
SCHULTZ (1924) remarks: ‘Among all primates, except in a large per- 
centage of white men and perhaps of some other human races, the fourth 
finger surpasses the second in length’. Woop JONES (1944) also describes 
the Uln. form as typical of all ‘monkeys and apes’. The Rd. type is 
‘definitely non-simian and constitutes a characteristic human specialisation’. 
‘He also remarks that, although the Rd. type is found only in man ‘this 
formula is found only in a certain number of cases, that it may be present 
in one hand and not in the other and that it depends upon the greater 
development of the index finger’. 

It is precisely this ‘elongated index’ that we may regard as a ‘distinctly 
human specialisation’. In this connection WoopD JONES points to the 
differentiation of a separate deep index flexor muscle from the musculus 
flexor digitorum profundus vel perforans. “The factor underlying the 
differentiation of this portion is undoubtedly the human specialisation of 
the index.’ 

From the literature it seems justifiable to conclude that, while the 


Rd. type may perhaps occur sometimes in anthropoids, the Uln. type is 
the rule in apes and anthropoids. 


6. Difference in prominence between the two hands. 
Little attention seems on the whole to have been paid to a possible 
difference in prominence between the right and the left hand. According 


to WEISSENBERG (1895) the Rd. type is commoner in the left hand. We 
quote the figures for his 574 Jewish boys and men. 
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RUGGLES (1930) comes to exactly the opposite conclusion, stating that the 
Rd. type predominates in the right hand in both men and women. 


Rd. type (WEISSENBERG 1895). 
rr rrr SSS 


5—10 yr. fy 3-230 yr. “4 21—30 yr. 31 yr. ond older 
Right n=30 | 45.5% | n=57 | 18.9% | n=29 | 23.60%) | n=21 | 25.00% 
Left n=29 | 43.9%/p | n=75 | 24.909 | n=39 | 31.7% | n=27 | 32.1/p 


We have already seen that Woop Jones (1944) says of the Rd. type 
that it may be present in one hand and not in the other’. 

According to RUGGLES (1930) there is no relation between hand type 
and ‘handedness’. 


Rd. type (RUGGLES 1930). 


Men ) Women 
Right | n=56| 28% | n=57| 520% 
Left | n=40| 20% | n=51 | 47% 


As is the case with practically all the aspects of the prominence problem 
which we have discussed here, a study of the literature once again presents 
us with contradictory opinions. 


t 
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WoLorzkol! (1924) determines the ‘mean’ prominence of the two hands 
of an individual by placing the hands, each with its longitudinal axis along 
the prolongation of that of the forearm, in the same plane, with the middle 
fingers tip to tip. If the distance between the tips of the index fingers is 
smaller than that between the the tips of the ring fingers we have a (mean) 
Rd. type. As this method doubles the difference in prominence between II 
and IV and thus shows it more clearly, we also used it for our investigations 
(see photo). 


Physiology. — Phosphate-exchanges in purple sulphur bacteria in con- 
nection with photosynthesis. By E. C. Wassink, J. E. Tyla and 
J. F. G. M. WinTERmans, (From the Laboratory of Plantphysi- 
ological Research, Agricultural University, Wageningen.) (21st 
Communication on Photosynthesis) *). (Communicated by Prof. 
A. J. KLUYVER.) 


(Communicated at the meeting of March 26, 1949.) 


Introduction. 


Since the time VOGLER and his collaborators made the important 
discovery that in cultures of the chemo-autotrophic sulphur bacterium 
Thiobacillus thiooxidans the shift from the energy-producing to the 
energy-consuming phase of metabolism is accompanied by a phosphate 
exchange (1), suggestions as to something analogous in photo-autotrophic 
metabolism have not been lacking (2—7). Already VOGLER himself put 
the question: “...... is it possible to irradiate photosynthetic organisms 
in the absence of COg and to store at least a portion of the radiant energy 
within the cell in a form which can later be used for COgz fixation in the 
dark?” (2). A certain uptake of CO, in the dark by photosynthetic 
organisms had already been demonstrated by sensitive methods (8), so 
that there appeared to be a fair chance that VOGLER’s question would be 
answered in the affirmative. Nevertheless, the attempts undertaken so far 
to furnish direct proof herefor cannot be said to have shown very definite 
results (5, 9). 

These attempts all made use of green cells. Now, in our opinion, there 
are some good reasons to give preference to purple sulphur bacteria for 
these studies. In the first place they are more closely related to the 
organisms VOGLER used. Some strains of purple bacteria even have been 
shown to be capable of an anaerobic, photosynthetic mode of life as well 
as of a heterotrophic, oxydative one. In the second place, in purple 
bacteria, both carbon dioxide and the hydrogen donor can be supplied at 
will. In the third place, previous studies had given some general idea of 
the kinetics of the metabolism of at least one strain (Chromatium, 
strain D), 

We, therefore, decided to look for possible connections of photosyn- 
thesis and phosphate exchange in suspensions of Chromatium, strain D. 


The method of phosphate determination was worked out chiefly by the second author, 
the results reported below were collected by the last mentioned author. - 


Methods. 
The bacteria were grown in the medium, described earlier as “combination 23”, 


*) 20th Comm.: Ann, Rev, Biochem, 17, 559578 (1948). 
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containing 0.24% sodium-l-malate and 0.16% sodium thiosulphate in an inorganic 
stock-solution (10). In the reported experiments, cultures from a small amount of dense 
inoculum were used after 1 day of incubation at about 27° in a light cabinet, They were 
centrifuged, suspended in the medium used in the experiments (see below), and centrifuged 
again. After this, the bacteria were resuspended in the same medium, and the experiments 
were started. 

The general trend of an experiment was very simple. A suspension of bacteria in a 
glass cylinder of about 20cm height and about 50 ml. contents was ventilated with a 
suitable gas stream, mostly consisting of oxygen-free nitrogen with either hydrogen or 
carbon dioxide added, and exposed either to light or to darkness. At the moment of 
changing conditions, or at intervals during an exposition, inorganic phosphate-P in the 
Suspension medium was determined. 

To this purpose, it was deemed advisable to separate the cells from the suspension 
medium as quickly as possible, and under sensibly the same conditions as during the 
experiment. In view of the experiments in light, centrifuging appeared unsuitable in this 
respect. At first, an experimental set-up was made in which the bottom of a glass cylinder 
was replaced by a cylindrical funnel with glass filter plate. Suspension liquid could 
be removed by suction during the fully undisturbed experiment. 

Unfortunately, however, bacteria entered into the pores of the filter plate — of about 
1 ~ width — soon preventing rapid filtering. After a few trials the following simple 
procedure was found satisfactory. Cylindrical glass funnels were made, about 5cm in 
diameter, and 6cm in height, with a nearly flat bottom and a tube of about 1 cm width, 
A stiffly rolled strip of ordinary filter paper, about 2cm broad and 30cm long, was 
pressed into the upper end of the tube. About 5 ml. of suspension were now removed 
from the glass cylinder in which the bacteria were exposed to the experimental conditions, 
introduced into the funnel — which was either illuminated or in darkness — and ventilated 
with the same gas stream as used in the experiment. Then suction was quickly applied, 
upon which the suspension medium rapidly went through clear. 

In order to have a chance upon measurable relative changes in PO, -P, solutions 
of low phosphate contents had to be used; in general about 10mg P/ml was taken. 
Initially, dilute phosphate buffers of pH ~ 8.0, according to CLARK and LUBBS, and to 
SORENSEN were used. Since the bacteria soon became inactive in these dilute media, 
1% NaCl was added. A remaining drawback was the weak buffering capacity of the 
medium. In case of ventilation with gas containing COs, the buffering capacity could be 
increased by addition of sodium bicarbonate, which, however, had the disadvantage of 
preventing subsequent removal of COz — e.g., when applying No-+ He afterwards — 
and, moreover, turned out to have an influence upon the phosphate determination. 

In search for other suitable buffering systems it has to be observed that systems 
containing organic acids, as e.g., acetate, citrate, efc. are less advisable since thése 
substances are likely to be used as sources of hydrogen and/or of COz by the bacteria. 
Finally, a borate mixture was found suitable and was stood very well by the bacteria; 
PALITZCH’ mixtures were used. Ten ml. solution containing 3 ml. M/20 borax, and 
7 ml. M/5 boric acid + M/20 NaCl, pH ~ 8.0 were added to 40 ml. of a phosphate 
mixture containing about 104g P/ml. When ventilated in dark with Ne+ 1% COs, a 
decrease of pH was observed which, after about 3 hours, remained at about 7,15. 
B-concentrations used did not interfere with P-determination. 

Phosphate-P was determined colorimetrically, according to the phospho-molybdic acid 
method first devised by OSMOND (1887). BELL and Doisy (1920) first used this method 
in connection with an organic reductant. Several modifications, differing in the sort of 
reductant, have since been described. We used that of LOWRY and LoPEz (11), with 
ascorbic acid as reductant. Advantages of this method are the operation at moderate pH 
and the use of a rather dilute solution of molybdate. 

Two ml. of suspension medium sucked through the filter (cf. above) were introduced 
into a measuring flask of 25 ml., and some acetate buffer, pH 4.0, was added, Subsequently, 
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Qml. of a solution containing 1% of ammonium molybdate in 0.05. sulfuric acid, and 
2 ml. of a 1% solution of ascorbic acid were added, and filled up to 25 ml. with acetate 
buffer, A stopwatch is started after addition of the reductant, and the blue color, gradually 
developing, is measured after 5 and 10 minutes with a ‘“lumetron, model 400 A” colorimeter, 
using the red filter, and a solution without molybdic acid as reference. In general, the 
readings after 5 and 10 minutes turned out to be the same, so that each moment between 
was suitable for the determination. About 10 readings were taken from each sample, 
yielding one measurement. The phosphate content according: to the reading was computed 
from a calibration curve obtained by submitting solutions of known phosphor contents 
to the same procedure. 

The accuracy obtained may be indicated by the following figures: 20 measurements 
of blank P+ B-solution (before introducing bacteria into it) yielded P = 8.9 + 0.05; 
6 = 0.23 4g P/ml. 

The following experimental details still have to be mentioned. In the earlier experiments 
glass cylinders were used to take up the bacteria in the experiments (cf. above). As a 
rule two of these were used in each experiment, submitted to different conditions (e.g., 
light versus darkness, Nz + He versus Ne + CO2). The suspensions contained about 5 cmm 
bacteria/ml. The cylinders stood in a thermostatic water bath with glass sides, at about 
29°, illuminated from 2 opposite sides with a 100 Watt incandescent lamp, yielding on 
either side an intensity of 8—9000 lux. In later experiments flat glass boxes of about 
20 X10 X 0.6.cm were used instead of the cylinders. , 

During the experiments a flow-meter-controlled gas stream from a bomb was passed 
through the suspensions. Either Ne + 30% He, Ne + 30% He + 5% COs, or 
No + 1% COs, obtained by mixing pure nitrogen with nitrogen containing 5% of carbon 
dioxide, were applied. The gasses were freed from oxygen by passing them over electrically 
heated copper gauze. With the aid of flow meters the gas streams in both cylinders or 
glass boxes were adjusted to equal velocity. At the top the vessels were closed except 
inlet and outlet of gas, the first ending in the suspension near the bottom. Before entering 
into the vessel with the bacteria, the gas passed a washing bottle with water at room 
temperature, to control evaporation of the suspension. 

Unfortunately, so far, we were not in a position to follow simultaneously the gas 
exchange of the cells used, owing to the lack of suitable apparatus herefor. This is planned 
for a subsequent part of this investigation. 


Experimental results and their preliminary discussion. 


The first experiments tended to “translate” VOGLER’s crucial experiments 
in terms of photosynthesis as closely as possible. A suspension of Chro- 
matium, Strain D, in a borate mixture. with about 10 ug PO,—P/ml., 
pH ~ 8.0, under No'+ 30\% Hz» was submitted to a short dark period for 
“adaptation”, and then illuminated. During the illumination period, a 
decrease of phosphate in the external medium is observed. After some 
hours, the light is turned off, and the gas phase replaced by No + 5 % COg. 
Then, an increase in external phosphate is observed. The procedure may 
be repeated with the same result (fig. 1, curve a). It is obvious that these 
observations are concordant with a concept, analogous to the one proposed 
by VOGLER for the chemo-autotrophic sulphur bacteria. In the light, 
hydrogen may be taken up, and fixed in some form, involving phosphate 
uptake. If this form is an “energy rich phosphate”, able to aid in the 
reduction of carbon dioxide in the dark, release of phosphate in a sub- 
sequent dark period may be expected. 
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It may be remarked already now that, in a M/15 phosphate buffer of pH ~ 8.0, 
WASSINK and KUIPER observed a long-lasting uptake of gas from No-+ He, in the 
absence of COs, even in cells submitted to special “starvation” treatments (in an investig- 
ation on the relation between redox potential and photosynthesis [12—14]). So far studies 
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Fig. 1. Phosphate changes in suspension medium with Chromatium, strain D, in relation 
to light and darkness, and under various gas phases. Exp.S of 13, 14.10. 48. 
+ : shift to darkness, + : shift to light. 


as to a “dark pick-up’ of COs have not yet been made. In the course of the mentioned 
investigation some incidental indications were obtained that cells, allowed to take up 
hydrogen in the light for a long time, already quickly thereafter produce considerable 
amounts of COs in the dark. 


Not always, however, the cells showed clear responses as to phosphate 
in the various phases of the above mentioned type of experiment (c/., e.g., 
fig. 1, curve 6). This forced us to a more systematic study of the various 
phenomena involved. 

The next observation made was a release of PO, in the dark (increase 
in the external medium) when freshly harvested, and washed cells were 
suspended as described. This PO,-liberation starts at once, and shows a 
tendency to decrease with time. In N, + 30 % Hg, and in No + 1 % CO, 
it was studied in greater detail and only little difference between these 
two gas phases was found (fig. 2, 3, 4, 5, 6, 8, 10; Table I). This release 

TABLE I. 
Release of POs—P (increase of P in suspension medium) by Chromatium, strain D, 


in darkness under various gas phases. 
Suspension medium: borate buffer, pH ~ 8.0, with ~ 10 4g P/ml. 


Ventilated with . Time | Total PO4-P released Number of observations 
Nz + 30/9 H2 35° eke fof tae ae OP 15 
120’ 1.1 + 0.13 12 
180/ 1.8 + 0.14 5 
285’ Zale tie OLZ5 2 
b N2-+ 19/9 CO2 120’ One EiALe 12 
270' Del. 0555 5 


of PO, probably may be compared with that accompanying “endogenous 
respiration” in VOGLER’s experiments. Notwithstanding the fact that, 
owing to the absence of oxygen, the Thiorhodaceae can not show a 
respiration-proper, some form of energy liberating metabolism is likely to 
continue in darkness, which may well be accompanied by a release of POx. 
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Since neither hydrogen nor carbon dioxide are likely to be active in the 
dark in these bacteria — unless, eventually, after special treatments — it 
is perhaps not very surprising that the rate of release is fairly independent 
of the gas phases applied. 

Next, the behaviour of the cells after admitting light, in Nz + 30 % He, 
was studied in greater detail. The phosphate release is now replaced by 
phosphate uptake (fig. 1, 2, 3, 4, 7b, 8b); after a short dark period, when 
the tendency to release phosphate apparently still was strong, the uptake 
at first was sometimes small or even slightly negative still (cf. fig. 4, 7) 
especially in dense suspensions (low average light intensity). Some typical 
experiments are represented in fig. 2 and 3; in each case, after a short dark 
period, the two halves of the same suspension were submitted oppositely 
to light and darkness, whereas, some 2 hours later, the conditions were 
reversed (fig. 2) or continued (fig. 3) without changing the gas phase. 
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to light and darkness, Flushed with Nz + 30% Hp. E: Exp. of 17.11.48; B: Exp. of 10.11.48. 
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Fig. 3. Phosphate changes in suspension medium with Chromaftium, strain D, in relation 
to light and darkness. Flushed with No + 30% He. Exp. of 9. 11. 48, 
+ : shift to darkness, + : shift to light. 


A number of experiments of the type-fig. 2 is summarized in fig. 4, showing 
the different behaviour in light and darkness very clearly. It may be con- 
cluded too, provisionally, that the rate of phosphate uptake in the light is 
higher after a longer dark period (fig. 4, Table II, a2). The reaction upon 
renewed darkness will be discussed below. 

Also in Ng + 1:% COs, illumination as a rule causes uptake of phos-. 
phate, replacing the release occurring in the dark (fig. 5, a few experiments 
are summarized in fig. 6). Here, too, after a prolonged dark period the 
rate of uptake appeared increased (cf. Table II, b). 
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Fig. 4. Phosphate changes in suspension medium with Chromatium, strain D, in relation 

to light and darkness. Flushed with No + 30% Hoe. Aliquots of the same culture indicated 
by equal characters. Exp.S of various dates, 

+ : shift to darkness, + : shift to light. 
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Fig. 5. Phosphate changes in suspension medium with Chromatium, strain D, in relation 
to light and darkness. Flushed with Noe+1% COs. Exp. of 2. 12.48. 
+ : shift to darkness. + : shift to light. 
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Fig. 6. Phosphate changes in suspension medium with Chromatium, strain D, in relation 
to light and darkness. Flushed with N2+1%COs. Aliquots of the same suspension 
indicated by equal numbers. Exp.S of various dates. 
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TABLE II. 
Uptake of PO4-P (decrease of P in suspension medium) by Chromatium, strain D, 
in illumination under various gas phases. . 
Suspension medium: borate buffer, pH ~ 8.0, with ~ 104g P/ml. 


eT 
ee 


Condition during | Time of 
Condition during | Duration of P cabicean 9 exposure Total sata derare< = 
dark period dark period iffuainaton to light taken up in light | observatio 
Nz + 309/o Ha 30' N2 + 309/o H2 75. 0.29 +0.1 ug/ml. ply) 
135° 0549 4 OS 12 
195" 0.80) 0532" 5 5 
235 1235 tO TOn rs 3 
. same 2 hours same 120’ 0.8042 0.26 6 
240’ | 1.11+0.16 5 
same 3 hours same 120’ 1220 == 0710 a 
same 4.5 hours same 120’ 2.402 OF2 2 
No + 1%CO,) 2 hours | N+ 1% CO, 120’ | 0.31+0.16 1 
b 240’ O20 2057 
same 4.5 hours same 120’ Le Oats Ole | A 
Not 3000 Hp 30° N2 + 300/y Hy 
+ 5/9 CO 150’ OL Ra O14 5 


However, another feature is very obvious, viz., that the uptake in light 
under No + 1 % COg is much smaller than under Ny + 30 % Hy (fig. 8). 
In some cases with CO, it even hardly reaches significantly positive values 
(cf. Table II, 6). The same can be said, comparing Ny + Hg + CO, with 
No + Hs, also here the presence of COzg reduces the phosphate uptake in 
the light practically to zero, and the difference in the general trend of the 
curves is very obvious, notwithstanding some incidental exceptions (fig. 7, 
Table II, c). 

It should be emphasized that most of the curves of a and b in fig. 7 and 8 
have been obtained one by one in the same experiment so that the variable 
“activity” of the cells to take up and give off phosphate cannot have 
interfered with these results. 

In the type of experiment as shown in fig. 2 and 5 (cf. also fig. 4 and 
6) darkness was again given after light in one of the two aliquots of 
suspension. The number of these observations is still too small to draw 
definite conclusions. It is obvious, however, that phosphate uptake stops 
or is converted into release. It would seem that in Ny'+ COg the release is 
somewhat more definitely pronounced than in No + He (cf. fig. 4 and 6). 
The reaction upon darkness after light without changing the gas phase 
throws some doubt upon the meaning of the release of phosphate after 
changing from Hy to CO,. More and quantitative studies, taking into 
account the gas exchange, are required here. In general, however, it may 
be admitted that under all conditions the cells are provided with organic 
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substances, which may be converted under release of phosphate and pro- 
bably of COs, as soon as conditions for active rebuilding of ~ph are 
absent, especially when the source of energy: light, is lacking, 
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Fig. 7. Phosphate changes in suspension medium with Chromatium, strain D. Influences 
of various gas mixtures during illumination and eventually subsequent darkness. Aliquots 
of suspensions were separated after + } hour adaptation in dark with Nz + 30% He. 
(Numbers indicate parallels with the same culture.) Exp.S of various dates. 
+ : shift to darkness. + : shift to light. 
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Fig. 8. Phosphate changes in suspension medium with Chromatium, strain D. Influence 
of various gas mixtures in dark and in light in parallel experiments (indicated by 
characters.) Exp.S of various dates. 
+ : shift to darkness. * : shift to light. 


General Discussion. 


In the foregoing, uptake and release of phosphate have been observed. 
It is obvious, a priori, that neither of the two can last infinitely with the 
same speed, even under as such favorable conditions. It is obvious too, 
therefore, that both uptake and release will show an asymptotic course in 
relation to time. Furtheron, it may be admitted that the rate of reversion 
upon a change of conditions also depends upon the nature and duration 
of the condition before the change. Finally, if indeed there is a relation 
between phosphate metabolism and photosynthesis of the type VOGLER 
found in Thiobacillus thiooxidans, we may expect that in the light, in the 
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absence of carbon dioxide, cells will accumulate phosphate to a higher 
degree than in its presence. In darkness, we may expect phosphate bond 
energy to be released as part of the energy reserve of the cell, used up in 
darkness. 

From these general statements, and founded upon the results of VOGLER 
and UMBREIT, we would expect purple sulphur bacteria to influence the 
phosphate content of a limited amount of medium in which they are 
suspended, in a way schematically represented in fig. 9. 
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Fig. 9, Schematical representation of expected POg-shifts under different conditions. 
See text. 
+ : shift to darkness. + : shift to light. 


Cells brought into darkness from a well-fed condition will start releasing 
phosphate from a level A, which phosphate accumulates in the medium. 
Since neither hydrogen nor carbon dioxide are known to act as energy 
sources in the dark for Chromatium, there is no reason to expect that the 
release of phosphate will depend either on the presence of hydrogen or of 
carbon dioxide. The release will show an asymptotic course (AB). Upon 
illumination, e.g. in Ng + Ho, an uptake of phosphate will occur which 
may be expected to start at higher rate, the more phosphate has been 
released before, so, e.g., the longer the dark period has lasted. The curves 
starting at C, D, E, etc., will be expected to strive towards virtually the 
same level, indicated by H. The probability that this level is, e.g., dependent 
on light intensity, will not be discussed here. Under Ny + CO. we may 
expect a behaviour fundamentally analogous to that under Ny + Hg, only 
pointing to a significantly higher ultimate phosphate content in the external 
medium (level K). 

In fig. 10 we have collected from our total experimental material, the 
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Fig. 10. Summarized representation of actual POy-shifts from table I and II. 
Cf. fig. 9 and text. 
. + : shift to darkness. t : shift to light. 
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data pertaining to the points outlined above. Notwithstanding the fact that 
on various points the material still is somewhat scanty (expressed by the 
large mean errors of some points and the lack of data on shorter intervals) 
it is clear that, in general, fig. 10 shows a very distinct analogy with the 
scheme of fig. 9. This would seem to indicate that, in general, the situation 
answers to the outline given. 

It should, however, not be overlooked that the experiments made so far, 
do not prove that, in photosynthesis of Chromatium, in the absence of 
COs, part of the energy of the light actually is accumulated as phosphate 
bond energy. To prove this would require recovery of these energy rich 
Phosphates from the cells which determinations have not yet been 
attempted. The findings merely are very well in accordance with the 
supposition that indeed light energy is accumulated in this form. 

It may be useful still to emphasize a rather fundamental difference 
between VOGLER’s case and that of the photosynthetic bacteria. In Thio- 
bacillus, the oxidation of the hydrogen donor is achieved by oxygen, and 
is, or at least may be conceived to be fundamentally independent of carbon 
dioxide. This easily leaves room for a storage of energy from the oxidation 
Process in some form which may, or may not, be used for the reduction of 
carbon dioxide. In Chromatium, on the other hand, carbon dioxide itself 
is the ultimate oxidant of the hydrogen donor; for the transfer of hydrogen 
light is essential. This, however, would seem to leave much less room for 
an independent storage of (light) energy at the expense of the hydrogen 
donor, and requires a much more specialized mechanism. The observation, 
alluded to above, that excessive hydrogen uptake may occur without supply 
of carbon dioxide (13, 14), indicates that indeed under certain conditions 
— e.g. pH ~ 8.0 — hydrogen may be assimilated more or less independent 
of carbon dioxide. It is interesting in this respect that so far no indications 
were found for an appreciable CO,-uptake without hydrogen during 
prolonged illuminations. The results reported in the present paper suggest 
that the mentioned uptake of hydrogen — which was not measured in the 
present study — is connected with phosphate uptake. It is perhaps most 
plausible to assume, that phosphate-shifts are connected with the transfer 
of energy from the pigment-protein system (cf. [15, 16, 12]). 

In view of an estimation of the quantitative importance of the phosphate 
shifts observed, it may be of interest to compare their magnitude with the 
presumable phosphate content of the cells. If we suppose a “normal” 
phosphor content in Chromatium-cells, this may be fixed at ~ 5 % of the 
dry weight, as P2O;; the dry weight may be 15 % of the fresh weight. 
This would mean that the bacteria contain 0.75 % of their fresh weight as 
P,O;, or 0.33 % P. The suspension density as a rule was about 5 cmm/nl., 
or about 5 mg/ml. Thus, the bacteria in 1 ml. contain about 17 ug P. On 
exposure to Ny + Hg in the light, an uptake of 1.5 wg P/ml is not rare. 
It thus appears that the bacteria may underly shifts in phosphate amounting 
to about 10 %. In order to achieve that these shifts will definitely surpass 
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the limits of the method of determination, rather dense bacterial suspensions 
are advisable. In order to obtain sufficient illumination inside the suspension, 
in our more recent series we used glass boxes instead of cylinders to 
suspend the bacteria. 


Summary. 

In suspensions of Chromatium, strain D, in borate buffer, pH ~ 8, 
containing some phosphate, shifts in phosphate content of the medium 
were found correlated with shifts from light to darkness, and from a gas 
phase consisting of Ng + Hz to one containing Ny + COnx.- 

The principal observations were as follows. Bacteria taken from the 
culture medium (“comb. 23”, see [10]) release phosphate in darkness, at 
about equal rates under No + He as under No + COs. Shift to light 
results in a marked uptake of phosphate under No + Hy and in a much 
smaller one under No + COg. Shift from light under No. + Hg to darkness 
under No!+ CO, results in a marked release of phosphate. Change from 
light to darkness without change of gas phase results in a certain release, 
too. 

As a whole the data collected so far show a marked similarity with 
VOoGLER’s findings with Thiobacillus thiooxidans, if No + Hg in light is 
parallelized with sulphur oxidation. They thus form a preliminary support 
for the hypothesis that Chromatium in the light is capable of building up 
energy rich phosphate bonds which are broken down in part when CO, 
becomes available and/or light is withdrawn. The phosphate exchanges 
observed so far can amount to 10% of the phosphate content of the 
bacteria. 

Correlations with gas exchange were not made so far. The investigation 
is being continued. 
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Botany. — De hypothese voor de erfelijkheidsformules van de twee 
zuivere lijnen I en II van Phaseolus vulgaris op grond van kruisings- 
proeven. I. By G. P. Frets. (Communicated by Prof. J. BOEKE.) 


(Communicated at the meeting of January 29, 1949.) 


In een vroegere mededeling (Proceed. 50, p. 798, 1947) stelden we 
voor de bonen van de I- en Il-lijn de formules L, L,... Le Lg By, By... 
Bs Bs baby . .. be bg thy thy ... the the en 1,1, ...Igl¢b,b,... bg bg... 
Th, Th, ... Thg Thg thy thy... thgthg op, of eenvoudiger geschreven 
L, Ly By be th, thy en 1, ly by by Thy Thy. De bonen van de I-lijn zijn lang, 
breed en dun, die van de II-lijn zijn kort, iets minder breed en dik. Het 
verschil van de lengten is veel groter dan dat van de breedten en de 
dikten. Bij de formules worden voor de erfelijkheid van de verschillen der 
afmetingen 6 paar genen aangenomen. De grootte-toename door één gen 
is voor alle 3 afmetingen relatief even groot. Steeds zijn één gen voor de 
lengte, één voor de breedte en één voor de dikte tegelijk werkzaam. Om 
het zoveel grotere verschil van de lengten dan van de breedten en de 
dikten te verklaren, wordt aangenomen, dat bij de bonen van de J-lijn er 
6 genen voor de lengte in de homozygote, dominante vorm (als LL) aan- 
wezig zijn, terwijl er slechts 3 genen voor de breedte in de homozygote, 
dominante vorm (als BB) aanwezig zijn. Bij bonen van de Il -lijn zijn er 
slechts 3 genen voor de dikte in de homozygote, dominante vorm (als 
Th Th) aanwezig. Van de 6 genen voor de afmetingen komen dus de 
genen 4—46 bij de breedte en bij de dikte alleen in de homozygote, reces- 
sieve vorm (als bb en thth) voor. Daarom schrijven we de formules ook 
L, Ly By, bs th, thy voor de bonen van de I-lijn en 1, lp by bo Th, the voor 
de bonen van de II-lijn. Volgens deze schrijfwijze hebben we te doen met 
een kruising volgens het tetrahybride schema. Na de kruising ontstaan er 
in de splitsingsgeneraties verschillende niewwe vormen als resultaat van 
genencombinaties op grond van het kruisingsschema zo b.v., voor de lengte 
en de breedte de vormen L, Ly By, bg en Lyly Byby. Er zullen dus, — 
we wezen hier reeds op in onze eerste mededeling (lc.) —, erfelijke 
variaties met een zeer grote lengte en een grote breedte en ook die met een 
niet zeer grote lengte en een grote breedte voorkomen, m.a.w. onder de 
bonen met de grootste breedten (B,by) zullen er zijn met zeer grote 
lengte (L,L2), doch ook met een niet zeer grote lengte (L,1,). De 
grens tussen bonen met een zeer grote en die met een niet zeer grote 
lengte, legden we bij 15.5 mm. Van bonen met een grote breedte is de 
breedte groter dan 8.5 mm. Bonen van de 2e groep, dus met de form. 
L, 1, Byby hebben een hoge LB-index. 

In onze eerste mededeling voegden we hieraan de opmerking toe, dat 
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individueel voortkweken van verschillende varianten de beslissing over de 
juistheid der hypothese zou kunnen brengen (l.c., p. 804). Met het oog op 
dit gezichtspunt bespreken we hier de varianten van de lengte en de 
breedte L, Ly By by en Ly lp By bg van F3-1935. We bedenken daarbij, 
Je dat er duidelijk erfelijkheid is bij de op elkaar volgende generaties, — 
er is een vrij grote positieve correlatie, l.c., p. 800 —, doch er is ook een 
grote niet-erfelijke variabiliteit, blijkende uit de hoge standaard-deviatie en 
variatie-coéfficiént (tab. 1 en 2). De grote niet-erfelijke variabiliteit zal het 
resultaat van ons onderzoek minder evident maken. Ook is het volstrekte 
bewijs der juistheid onzer hypothese niet te leveren. 2e De grond voor onze 
aanname, dat steeds een gen voor de lengte, één voor de breedte en één 
voor de dikte tegelijk werkzaam zijn, vinden we in de positieve correlatie 
tussen de afmetingen voor een willekeurig materiaal en voor de negatieve 
correlatie in groepen bonen met een zelfde gewicht, als uiting van ,,com- 
pensational growth”. Er is bij de groei der bonen een regulerende factor 


TABLE 1. The variation of beanyields of the F3- and the Fy-seedgenerations and of 
the I- and the Il-line of 1935. (Only the variation of the length is given in this table.) 


The length 
Pl. n M+m otm Ww) oF iin Var. r. 
I-line 11 | 39 |16.59-40.24| 1.54+0.18|] 9.3 18.5 |12.0 | 6.5 
3 | 38 |15.19+0.26] 1.61 +0.18| 10.6 | 18.3%] 11.13] 7.2 
9 | 58 |14.334+0.22] 1.59+0.15| 11.1 |19.1 |11.1 | 8.0 
7 | 60 [14.3 +0.14| 1.08+0.1 6 |17.5 |11.3 | 6.2 
IJ-line 22 48) || 1208s 0.09) 0.64 4-0,07 15). ai 1s20 958 15352 
25 8720). 11 502-009) 087 Se 0207 e795) 12.9 8.6 | 4.3 
30 30. 10.49 400103) OF73 2001 7 11.9 8.677) 3.5 
F3-seedgene- ie) 351 159° SOL 24 1539 = OF Fsree se Sea ee sO oes 
ration 38 52 14 58s OcT 7d 210.12 eee Maced oe ee 
44 90. 4/°13'.66 +:0,.09 |.0.84 40.06); G22 15.4 WOc oa eas 
61 122... W12 72s 0.061) (0. 0S 05044) 25.0 15 OO a ais 
97 58. 112.01 +-.0.92)) O94 2 O09) 7 Balls sa LO. aed 
F4-seedgene- | 119 51. 16.32 sO 1818 = Oe tier 2 Leo lS cere 
ration 133 5115.16 = OF11 5)" 0678 Se 00ST S22 OF Sell 18 aan 
150 46 114.2) 0316) Il) OM ror Ges WE See 
178 450 13,266 0711) .0.72 + 0.08) $.4.114.5 1 1) 6 cee 
187 42 NAD O Od 0:63 320.07 10 S22) S35, 1) LOS ese 
198 41 IOs SOc 0.67 +0.07) 6.1 | 1258 bo Vek ye TRIS Yen 
1) Then follows 17.6. 2) Then follows 12. 
The variation of the beanyields of the pure lines I and II, — that is non-hereditary- , 
variation —, is great; it is greater of the I-line than of the II-line. The variation of the 


beanyields of the F3- and the F4-generation, that is hereditary and non-hereditary variation 
is not greater than that of the I-line and of the Il-line (cf. tab. 2, p, 425). 
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TABLE 2. The mean dimensions and weights and the variation of the beans of the I- and 
the II-line in 1935—1937. (Only the lengths and the weights are mentioned in this table) 1), 


The length “a iy 
(| eatee a ile SI AEE ll inc 


The lines o = a m | V 4 . Var. fr. 
1935 Iline | 774 | 15.2 +0.07| 1.78+0.05| 11.7 | 21.62)| 10.1 | 11.5 
1935 | Wine | 802 | 11.1 +0.03/ 0.84+0.02] 7.5] 13.7 | 8.63! 5.1 
1936 I-line | 1617 | 14.2 +0.02) 0.99+0.02/ 7 | 18.1 | 10.2 | 7.9 
1936 | IFline 605 11.16 +0.03 | 0.814+0.02| 7.3/13,3 | 7.3 | 6 
1937 | Ilime | 2383 | 15.2 +0.02) 1.13+0.02! 7.4/ 19.0 | 9.7 | 93 
1937 | I-line | 1292 (11.0 +0.02| 0.80+ 0.02) 7.3| 14.93 6.4 | 8.5 
ESS Eas ee ae Te Sas es ea 

The weight 

1935 Lline 895 |52.5 +0.5 116.4 +0.4 | 31 | 1274] 10 | 117 
1935 | Illine | 799 |42:4 +0.3 | 7.9 40.1 | 19 63. | 15 | 48 
1936 Iline | 1610 |50.7 +0.3 |10.4 +0.2 | 21 89 21 68 
1936 | I-line | 605 |45.1 +0.4 | 8.6 +0.2 | 19 66 10 | 56 
1937 line | 1303 |63.85)+ 0.4 |13.1 +0.3 | 21 | 99 | 22 | 77 
1937 | W-line | 1298 [43.8 +0.2 | 6.5 +0.1 | 15 64 15 +} 49 


1) The F3- and following generations are not included in the table. The averages 
are not comparable here: at the basis of the material is the selection of the initial beans. 

2) Then follows 19.6 mm. 

%) Then follows 13.2 mm. 

4) Then follows 105 cg. i 

5) The mean weight of 1937 is much greater than that of 1935. The mean thickness 
is also very great, whereas that of 1935 is small. 


tussen de afmetingen (1947,l.c.). 3e Onder de uitsplitsingen van de F3- en 
F4-zaadgeneratie zullen ook bonen als de oudervormen zijn. Dergelijke 
uitzonderlijke uitsplitsingen zullen echter phaenotypisch als regel verschil- 
len van de uiterste varianten van de I- en de II-lijn (1947). 

We groepeerden van het materiaal van F3-1935 de bonen met de grote 
breedten volgens breedte-klassen en tekenden de bij de bonen van iedere 
breedte-klasse behorende lengte aan. Hetzelfde deden we voor vergelijk- 
materiaal van de I- en de II-lijn. We gaan dan na, of in de groepen van 
bonen met de grootste breedten meer bonen met een niet zeer grote of 
kleine lengte voorkomen dan in de overeenkomstige groepen van bonen 
van de I-lijn. We vinden het volgende: In tab. 3 zijn van het Fz-materiaal 
van 1935 de bonen met de grootste breedten opgenomen. Ter vergelijking 
zijn bonen van de I-lijn en van de II -lijn er naast geplaatst. We zien 
(tab. 3), dat de grootste breedte van bonen van de II -lijn van 1935 is, 
b = 9.8 mm. De bonen van de II -lijn komen welhaast voor ons onderzoek ° 
niet in aanmerking. Er zijn in het I-materiaal 5 bonen met een breedte van 
b = 13.0—11.5 mm, dit zijn grotere breedten dan bij de F3-bonen van 
1935 voorkomen. Al deze 5 bonen hebben een zeer grote lengte; ook der- 
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TABLE 3. The Fs-seedgeneration of 1935. The number of beans with a length that is 
not very great (15.5 mm and smaller) and with a very great length (15.6 mm and greater) 
for beans with a great breadth (8.6 mm and greater) and comparison beans of the I- and 


the ILline. sm. 1. = number of beans with a length of 15.5 mm and smaller. gr. 1. = number 
of beans with a length of 15.6mm and greater. gr.v. = greatest variation; sm.v. = 
smallest variation. 
1935 
a 
F3-seedgeneration I-line !) 
Br.>| n |Sm.1.|Gr.l.} m, |Gr.v.|Sm.v.]) o Sm. 1.|Gr.1 my Gr v.| Sm.v. 
130 1 0 1 216 
121 1 0 1 179 
17, 1 0 1 193 
als) 2 0 2 184.5] 186 183 
113 1 1 0 Ney #2 0 2 179 191 167 
112 1 0 1 183 5 0 5 180.8| 196 161 
111 1 0 1 164 5 1 47 Olof 148 
110 a 1 3 NCES) BA 9) Ge: 6 0) 6 173.8 183 162 
109 3 1 2 SOEs Oma tak 6 0 6 168.8] 182 158 
108 4 0 4 | 167.8) 172 | 160 9 0 9 170.9} 183 156 
107 ii a ey MYST ier ey I 2 13 164.7| 180 146 
106 2 1 1 149.5) 158 | 141 9 1 8 174.7| 185 155 
W@sy) hy 3 OM i) 1625 3/1835 1465 Hy let 10 166.3) 177 148 
104 | 10 6 AW 15490) 166) tad 18 17 168.8) 186 147 
LOSS 4 se 19 15 154.6] 176 | 137 || 20 2 18 168.8} 196 153 
OR) Sil ake |) aol NSwigtss| 1sxey If tlsis) 19 0 19 165.4} 181 157 
Op” see i 2 LOS S15 Si LSet se Si? 1 11 17a 52) 183 155 
100 | 47 | 38 9 149 174 130" ele 1 16 167.1} 184 148 
OOF Gl) 1 49 S127 | T48eG) 182. 5) 132 16 0 16 166.9] 176 160 
SS WiOS= NWS AN L49R SI 75 136 || 26 Sy Er 165224) 192 142 
OTe ofan eeO 17 147.1) 169 127 1 39 £34 IGSOGI) LS 151 
96 1134 (115 19°) 146.5) 167 132° 120 6 14 160.5| 181 140 


1) We find for the material of the II-line that there are 2 beans with the great breadth 
(ie. the greatest breadth that occurs) b = 9.8mm, the length is resp. 12.0 and 11.8 mm, 
There is one bean with the breadth, b = 9.7 mm; the length is 1 = 13.0mm. There are 3 
beans with the breadth, b = 9.6 mm; the length is resp. 12.6, 12.4 and 12.2 mm. 


gelijke zeer grote lengten ontbreken onder de F3-bonen. Onder de 57 bonen 
van de I-lijn van 1935 met een breedte van 11.3—10.6 zijn er twee met een 
lengte van resp. 14.8—14.6 mm, d.w.z., dat er onder de bonen met de 
grote breedten 11.3—10.6 mm, van de I-lijn, d.i. van bonen met de form. 
L, Ly, B,by van de lengte en breedte, er 2 zijn met de niet zeer grote 
lengte 14.8 en 14.6 mm, als niet-erfelijke varianten hier. Bovendien is er 
één boon met de lengte 1 = 15.5 mm, d.i. de grenswaarde. Er zijn 23 bonen 
in het F3-materiaal van 1935 met de grote breedten 11.3—10.6 mm. Deze 
aantallen bonen van de F3-generatie en van de I-lijn zijn niet direct ver- 
gelijkbaar: het totale aantal bonen van de I-lijn is veel kleiner dan van de 
Fs-generatie. De I-lijn echter is een zuivere lijn en bevat slechts niet- 
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erfelijke variaties, terwijl de F3-generatie een zeer heterogene erfelijke 
samenstelling heeft. 

Onder deze 23 Fs-bonen met de breedten 10.6—11.3 mm, zijn 8 bonen 
met een niet zeer grote lengte (d.w.z. 1= 15.5 mm, en kleiner). Er zijn 
daarbij met een lengte van 14.1—14.4 mm, d.i. met een zo kleine lengte, 
als ze onder bonen van de I-lijn met deze breedten niet voorkomen. Onder 
de 11, resp. 18 bonen van de I-lijn met een breedte van 10.5 en 10.4 mm, 


L B Th W LB th Bth 


Fig. 1. a. Characterogram of 4 p 1b, initial bean of pl. 81, Fe-1934 for pl. 73, F3-1935. 
b. Idem of the averages of the beanyield of pl. 73, Fz-1935. S: = Standardcharacterogram 
of 1934. Se = Standardcharacterogram of 1935. 


vinden we in ieder van deze groepen slechts één boon met een niet zeer 
grote lengte, 1= 14.7 en = 14.8 mm. Onder de 12, resp. 10 bonen van 
de Fz-generatie van 1935 vinden we 3, resp. 6 bonen met een niet zeer 
grote lengte. De kleinste lengte van de Iste groep is 14.6, die van de 2de 
is 14.1 mm. 

In de volgende breedte-klassen b = 10.3 mm e.v. (tab. 3) blijven we in 
de groepen van de I-lijn slechts zeer enkele bonen met een niet zeer grote 
lengte aantreffen. De kleinste lengte, 1 = 14.8 mm, is van een boon met 
de breedte b= 10.0 mm. Van de bonen van de breedte-klassen b = 10.3 
mm, e.v. van de F3-generatie wordt de gemiddelde lengte regelmatig klei- 
ner; het aantal bonen in ieder van deze klassen is vrij groot. Van de bonen 
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met een breedte b = 10.3, b= 10.2 en b= 10.1 mm, is de kleinste lengte 
resp. 113.7, ==13.6) en 13-6 

We vinden dus onder de bonen met een grote breedte, b= 11.3—10.1 
mm, van de F3-generatie meerdere bonen met een niet zeer grote lengte, 
1 = 15.5—13.6 mm. 

Bij bonen met een dergelijke grote breedte b= 11.3—10.1 mm van de 
I-lijn van 1935, vinden we slechts een zeer enkele boon met een niet zeer 
grote lengte. De boon met de kleinste lengte is 1 = 14.6 mm. 

Bij het F3-materiaal zijn bonen met een zo kleine lengte, als ze onder de 
bonen van de I-lijn niet voorkomen. Deze verschillende bevinding laat zich 
verklaren uit de verschillende erfelijke samenstelling van deze groepen van 
bonen wan de F3-generatie en van de I-lijn, iedere groep met dezelfde 
grote breedte. De formule voor de lengte en de breedte is van alle bonen 
van de I-lijn dezelfde, nl. L; Lg By bp. Onder deze bonen van de F3-gene- 
ratie zijn er met de formule L, Ly By, by als van de I-lijn (of tot dit gebied 
behorend), maar ook met de formule L, ly By, bo. 

De geleidelijke afmame van de gemiddelde lengte der bonen van de op 
elkaar volgende breedte-klassen (tab. 3) wijst in de eerste plaats op de 
positieve correlatie van de lengte en de breedte, — we treffen de geleide- 
lijke afname ook aan voor de I-lijn —, en ook op de erfelijkheid door poly- 
mere factoren. L, en Ly, B, en by staan ieder voor een groep polymere 
factoren (Ly—L3 en Ly—Lg; By—Bs en by—be). 

Het aantal bonen van ons F3-materiaal is yoor de breedten 10.0—9.6 
mm, vrij groot, het aantal bonen met een niet zeer grote lengte neemt, 
vergeleken bij dat met een zeer grote lengte, zeer toe. Onder de bonen 
van de I-lijn met dezelfde breedten van 10.0—9.6 mm, zijn naar verhouding 
veel meer bonen met een zeer grote lengte. Dit verschil van bonen van de 
F3-generatie en van die van de I-lijn, berust hierop, dat in de formule van 
de lengte en de breedte van deze bonen met geleidelijk minder grote breed- 
te, de homozygotie der bonen afneemt. Het best voor ons onderzoek lenen 
zich de groepen der grootste breedten en hier vinden we een goede aan- 
wijzing voor het voorkomen, in het F3-materiaal, van bonen met voor de 
lengte en de breedte de formule L, 1, B,b. naast bonen met de formule 
L, Ly By by, d.w.z. onder de bonen met de grootste breedten zijn er met de 
formule L,; Lz By bo en Ly 1, By bo. 

In onze eerste mededeling (1947) groepeerden we de bonen volgens 
afnemende lengte en als aanwijzing voor de bevestiging van onze hypo- 
these, vonden we daar, dat in de groepen met de grootste lengten, de ge- 
middelde breedte dezelfde was; eerst in de latere groepen, daalde met de 
lengte, de gemiddelde breedte der bonen regelmatig. Ook hier de aanwij- 
zing, dat bonen met de grootste lengten (formule L, L,), dezelfde grote 
gemiddelde breedte (B,b.) hebben als bonen met minder grote lengten 
(formule L,1,). In tab. 4 zijn enige voorbeelden van F3-bonen van 1935 
met de grootste breedten en één met zeer grote lengte bijeengebracht en 
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TABLE 4. F3-1935. Some examples of beans with great breadths and not very great 
lengths and comparisonbeans of the I and the II-line. j 


BTh 

Ip 1b 148 107 66 70 72 45 62 

36 Ip 3b 146 107 65 66 73 45 61 
87 6p 2b 141 106 59 60 75 42 56 
68 9p 1b 146 105 77 83 72 53 73 
87 10p 2b 146 105 73 74 72 50 70 
87 6p 6b 141 104 71 74 74 50 67 
87 6p 4b 137 | 103 56 59 75 41 54 
58 5p 2b 139 102 65 67 73 47 64 
82 3p 3b 136 102 69 73 75 51 68 
65 9p 3b 137 101 76 73 74 56a e775 
103 | Ip 2b ) 140 101 65 65 72 46 64 

Comparisonbeans of the I-line (very rare) 

1 4p 3b | 148 | 111 65 78 75 44 69 
4 8p 2b | 146 107 66 70 73 45 62 
20 5p 2b | 148 | 105 61 64 71 41 58 
4 Ip 4b | 147 | 104 52 55 71 35 *| 50 
10 16p 2b 153 | 103 67 72 67 44 65 
Comparisonbeans of the Il-line (very rare) 

28! 5p 4b | 118 | 98 75 58 83 64 7 
28! 12p 2b | 122 | 98 74 59 80 61 76 
22 8p 2b b 490 97 67 55 75 52 69 


enkele zeer zeldzame vergelijkbonen van de I- en de IIJ-lijn. Deze Fs- 
bonen hebben een hoge LB-index, evenals de bonen van cl. 5 met de 
formule 1B Th van ons schema. De dikte is verschillend. Van bonen met 
een kleine dikte is de BTh-index laag, ze komen overeen met de bonen 
van cl. 6 met de formule 1 B th. Deze F3-bonen met een grote dikte hebben 
een matig hoge B Th-index. Ze komen overeen met bonen van cl. la met 
de formule L; lp BTh. Bonen met de formule L; 1, B, by voor de lengte 
en de breedte hebben een grote absolute breedte en een hoge L B-index. 

De hypothese (1947, l.c., blz. 800), waarbij aangenomen wordt, dat het 
grootte-verschil der lengten op de werking van een even groot aantal 
polymere factoren in de dominante homozygote vorin berust als dat der 
breedten en dikten, doch dat de werking van een factor voor het grootte- 
verschil van een lengte-factor beduidend groter is dan die van een factor 
voor de breedte en voor de dikte kan het voorkomen van bonen met een 
grote breedte en een niet zeer grote lengte ook verklaren, maar de ver- 
wachting van dergelijke bonen (dus met de formule voor de lengte en de 
breedte L, lp B, by naast die met de formule LL. By ba) is bij de door 
ons aangenomen hypothese veel groter. 

De verschillende erfelijke samenstelling van Fy-bonen brengt mee, dat 
de F3- bonenopbrengsten onderling zeer zullen verschillen. In tab. 5 zijn 
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enige bonenopbrengsten van de F3-zaadgeneratie met de grootste gemid- 
delde breedten en enkele vergelijk-bonenopbrengsten van de I- en de 
Il-lijn bijeengebracht. We vinden, dat van de bonenopbrengsten met de 
grootste gemiddelde breedten, de gemiddelde lengte van de bonenop- 
brengsten van de I-lijn steeds zeer groot is (1, — 16.6—15.8 mm, eenmaal 


LB Th W LB Lh Brh ‘a 


101 


183 


180- 


170- 


t B Th W LB Lh Bib 


160- 


154 


112 


150+ 


150 


140: 73 


130 + 


Age whe es aS 
132 87 64 49 67 50 74 128 85 


65 


40 


Fig. 2. a. Characterogram of 2 p I b, initial bean of pl. 73, F3-1935 for pl. 309, F4-1936. 
b, Idem of the averages of the beanyield of pl. 309, F4-1936. Ss = Standardcharacterogram 
of 1936, 


is 1, = 151.9 en = 149.3 mm). Onder bonenopbrengsten van de Fs- 
zaadgeneratie zijn er enkele, die, — ofschoon de gemiddelde breedte zeer 
groot is —, een niet zeer grote gemiddelde lengte hebben. De grootste 
gemiddelde breedte van bonenopbrengsten van de Il-lijn is 9.04 mm, d.i. 
kleiner dan de gemiddelde zeer grote breedten, die we hier beschouwen; de 
gemiddelde lengte is 12.08 mm. 


Bij de F3-bonenopbrengsten van 1935 met een grote gemiddelde breedte 
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zijn er enkele, met een kleine gemiddelde lengte, zoals ze bij de bonen- 
opbrengsten van de I-lijn met eenzelfde grote gemiddelde breedte niet 
voorkomen. Deze bevinding wijst er weer op, dat, onder de Fs-bonen- 
opbrengsten met een grote gemiddelde breedte, er met een andere erfelijke 
samenstelling van de lengte en de breedte zijn dan bij de bonenopbrengsten 
van de J-lijn. Van alle bonenopbrengsten van de I-lijn is de formule voor 
de lengte en de breedte L; Ly B,bo, onder die van de F3-zaadgeneratie 
zijner met de formule L,; L, B,; bs enook andere met de formule L, 1, By bo. 


Statistics. — Random sampling frequencies; an implement for rapidly 
constructing large-size artificial samples. By H. C. HAMAKER. 
(Laboratorium voor Wetenschappelijk Onderzoek der N.V. Philips’ 
Gloeilampenfabrieken Eindhoven—Nederland.) (Communicated by 
Prof. H. B. G. CASIMIR.) 


(Communicated at the meeting of March 26, 1949.) 


Artificial sampling is, no doubt, the chief purpose to which tables of 
random sampling numbers have been constructed and applied. To fix our 
argument, let us suppose that we wish to sample the Poisson distribution 
with parameter m — 2, which yields the cumulative frequencies listed 
below. 


TABLE I. 


Cumulative frequencies of the Poisson distribution with m = 2. 


0 1353 :5 .9834 
1 4060 6 .9955 
2 .6767 7 .9989 
3 .8571 8 .9998 
4 +9474 9 1.0000 


Then referring to a table of random sampling numbers arranged in sets 
of 4 digits, we score a 0 for each number lying between 0000 and 1353, 
we score a 1 for numbers between 1353 and 4060, a 2 for those between 
4060 and 6767, and so on; this is a common procedure which does not need 
further explanation. 

If it is our purpose to investigate certain properties of small samples, 
the taking of individual scores will be unavoidable, even though this may 
involve some tedious labour; but if we intend to study the properties of 
large samples, the work involved will soon become prohibitive, and the 
question arises whether it is possible to devise a more rapid and convenient 
method for the constructing of random samples of a large size. Such a 
method will be described below. 

We begin by observing that as we score a 0 for all numbers from 0000 
to 1353, we may score a 0 for all numbers beginning with 00 to 12 without 
any further knowledge of the two digits that follow; likewise all numbers 
beginning with a pair between 14 and 39 will score a 1 regardless of the 
two final digits. It is only for numbers beginning with 13, 40, etc., that 
we require to know the last two digits as well. 

Let us then for a moment neglect the last two digits . We take 1000 
pairs of two-digit random sampling numbers and we score these in a 
10 x 10 array, thereby obtaining the frequencies recorded in table II. 
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The classes and the frequencies in this table will be designated as the 
array-classes and the array-frequencies, to distinguish them from the 
sample-classes and the sample-frequencies relating to the sample that we 
are going to construct from them. 


TABLE II. 


1000 two-digit random sampling numbers scored in a 10 X 10 array, 
—— ES 


=e Second digit 
ane on, US es Ee a a a a a 

Array-frequencies 
0 13 4 9 10 95 
1 4 5 9 7 9 87 
2 12 9 a fil 10 90 
3 9 a 16) it 8 96 
4 8 8 sh) 04 15:1 408 
5 16 | 14 ! 14 9 7h AQ 
6 | 10 “a aad og ee 14 | 100 
7 7 6 | S 4, ts 14 94 
8 4 | 15 | 9 | 15 13, esi? 
9 8 Og oe 8 iy ae oan ie Ya ee 6) } 13.) 106 
Totals| 91 | 95 | 95 | 105,| 94 | 113 | 83 | 103 | 108 | 113 | 1000: 


With the aid of the array of table II we may now rapidly construct 
a random sample of 1000 scores of the population specified by table I in 
the following manner: 

The frequencies in the array-classes 00 to 12 are added giving a score 
of 113 entirely-belonging to sample-class 0. In the array-class 13 we have . 
a score of 7 and to allot these to sample classes 0 and 1 we have to tag 
two more digits behind them, and then to divide the resultant set of 
4-digit numbers into those below and those above 1353 (see table I); 
from a table of random sampling numbers we select a random set of 7 
two-digit figures, for instance 


1S he kr Pa | ia ale PE 


four of these are less than 53, and hence we decide that of the seven 
scores in array-class 13, 4 have to go into sample-class 0, and the other 
3 into sample-class 1. The total score in sample class 0 will then be 
$135: 4 4 =117. 

Next we add the frequencies in the array-classes 14—39, giving a 
total of 248; the score of 8 in array-class 40 is again subdivided by a set 
of 8 two-digit random sampling numbers, 3 going into sample-class 1, 
and 5 into sample-class 2. The total in sample class 1 is now also complete, 
viz: 3 (from array-class 13) + 248 (from array-classes 14—39) + 3 
(from array-class 40) = 254. 

Thus we proceed until we have reached the array-class 99, which has 
to be subdivided into four sets, those below 55, those between 55 and 89, 
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those between 89 and 98, and those above 98, going into the sample-classes 
6, 7, 8 and 9 respectively (compare table I). The sample will then be 
complete. The entire procedure is illustrated in table III. 


TABLE III. 


Illustrating the construction of a random sample of 1000 from the distribution of table I 


by means of the array-frequencies of table II. 
a 


Array Sample Th 
Sample class Array classes Feauracies hequancies eory 
0 0O—12 113 117 135 
4 
7 
7; 
13 N3 
1 14—39) 248 254 271 
3 : 
7 
8 
40 Ns 
2 A266 280 291 271 
6 
7 
6 
67 No 
3 68—84 178 \ 185 180 
7 
7 
9 
85 ce 
4 86—93 85 ; OZ 90 
5 
7 
4 ri 
9 No 
5 9597 40 43 36 
1 
98 6 cf 
N5 ) 
6 ( 9 
4 17 
7 99 13—>7 f 
8 N2 2 


q2 = 6.28; »= 6; P=0.395 


So far this method cannot claim to effect much saving of labour. In 
constructing our sample, however, we started in the left-hand top corner 
of the array of table II, and we proceeded from left to right and from the 
top to the bottom of this array; and a considerable saving of labour can 
be achieved by observing that there is no reason for using the array- 
frequencies precisely in this order. All the array-~frequencies are equivalent 
to one another, and in constructing a sample we may start at any point 
in the array and proceed in any direction; or, more general still, we may 
rearrange the frequencies of table I] at random in any of their 100! 
permutations before we start making up our sample. The simplest method 
to achieve this is to write the array frequencies on 100 chips of card 


board, mix these thoroughly and lay them out in random order while 
constructing the sample. 
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If we use invariably the same set of array-~frequencies (for instance 
those of table II) we may expect our samples to be on the average slightly 
biassed in one way or another, though the variability obtained by 
rearrranging is so large that serious bias is hardly to be expected. Never- 
theless it seems inadvisable to use one set only, and I have therefore 
constructed ten independent arrays (as that in table II), from which one 
may be chosen at random every time we wish to construct a sample; 
thereby the possibility of systematic bias should be ruled out. 


TABLE IV. 
Ten sets of random sampling frequencies for sample size 1000, 
SS 


Se ie aa ee ae 5 6 7 8 9 10 
Array-freq, f Frequencies of the array-frequencies = n 

2 = | 1 re Bema er. i ae heli 

ie 1 — 1 1 = 1 = 1 — 1 

4 2 1 4 1 3 1 1 3 5) 3 

5 3 ae ee a ee a eee! a ae ee ah 

6 4 6 5 8 3 7 9 9 6 3 

7 8 9 Gy diy 43 11 9 10 11 10 8 

8 11 11 8 17 9 a2 7 6 8 6 

9 16 14 15 11 15 9 -12 12 14 14 

10 17 9 17 13 12 19 14 13 10 13 

11 8 10 5 10 12 9 9 8 9 20 

12 9 7 a, Ae 3 1%) 11 10 3 4 

13 11 9 6 6 5 8 10 8 15 5) 

14 5 7 SD ees, 9 6 5 5 2) Hf 

15 1 5 7 6 6 1 2 6 5 4 

16 2 — os 1 3 1 3 4 2 3 

17 1 2 1 sae he Z — 2 1 2 2 

18 = = = = = 1 = ae 1 1 

19 1 2 — = = 1 — — _ — 

20 — — 1 = = = a = = — 

21 = — 1 — _ — — ~ -— — 
=n 100} 100} 100} 100; 100; 100; 100; 100; 100} 100 
nf 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 
S} 8.42 | 11.26| 11.50| 8.62 | 10.50} 8.20 | 9.14 | 10.54] 11.32] 10.70 


If we use the array-frequencies not in the form of table IJ, but as a 
lottery with 100 lots drawn without replacement, the actual array will 
be quite immaterial; it will only be important to know the frequencies of 
the array-frequencies. These I shall call the “random sampling [requencies’’; 
for the ten arrays actually constructed they have been recorded in table IV 
above; set Nr 5 corresponds to table II. 

Also in table V an independent set of random sampling frequencies 
have been provided for the construction of samples of size 500, and in 
table VI a set for sample size 200. 

Samples of other sizes may be obtained by combination; for instance, 
if we require a sample of 1500 we may construct two samples, one of 
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1000 and one of 500 scores, and then add their frequencies; or alternatively 
by laying out one set of 100 random sampling frequencies for sample 
size 1000 and one set for sample size 500 in a 10 X 10 array, and by 
adding the array~frequencies we obtain a set of random sampling 
frequencies for sample size 1500, from which by permutation any number 
of samples of this size may be constructed. If we desire a sample of 
exactly 1537 scores we have only to add 37 more made up in the usual 
way. Thus the sets of frequencies in tables IV, V, and VI should be 
amply sufficient for a great variety of purposes. 
TABLE V. 


Ten sets of random sampling frequencies for sample size 500. 


SL SSS SSS 


Set Nr 1 2 3 4 5 641-7 | 8 be 10 
Array freq. f Frequencies of the array-frequencies = n 

0 1 2)/—|- OT Bee) Mee) hears hoe 2 

il 5 3 2 3 6 — 3 5 1 2° 

2 9 9 5} 13 7 10 10 5 11 8 

3 13 13 15 20 10 19 21 19 2 13 

a 23 17 25 1h at 17 15 15 18 15 

5} 15 13 19 16 16 18 16 aa 17 20 

6 8 19 13 12 16. 13 10 12 20 18 

7 12 12 10 6 10 11 7 12 11 12 

8 8 6 3 5 5 4 10 7 Fi 5 

9 5 3 6 10 5 5 2 1 1 2 

10 1 i! — 2 1 1 3 1 — 1 

11 2 = 2 1 — 1 1 1 2 1 

12 — 2 — 1 1 1 = as = 1 

13 te = Se = 1 — 2 1 — — 
=n 100 | 100: |, 100 |} -100 |} 100) 2100 | 200], 100) 100} 5100 
Snf 500 | 500} 500} 500} 500} 500} -500 |} 500 | 500) 500 
sf 5:34 15,40 | 4.10 | 6.32 | 5.64 | 4.66 | 6.36 | 4.48 193.96 |) 4.74 

TABLE VI. 
Ten sets of random sampling frequencies for sample size 200. 
Set Nr | 1 2 3 4 5 6 7 he 1S one fo 
Array-freq. f Frequencies of array-frequencies = n 

ee ee ee ERE ORT Bea aes a 

0 18 12 16 14 10 14 16 12 14 14 

1 18 30 21 26 27 24 27 24 28 27 

2 30 25 32 31 31 33 29 32 24 30 

3 18 18 18 12 24 13 16 21 pa 12 

4 12 10 7 12 a 11 9 7 7 11 

5 4 4 3 a 3 3 5 3 3 4 

6 — 1 3 1 2 2 1 1 3 1 

if = _ — 1 _ — — — — 1 

8 — — = = — = 1 = — — 
Z 100; 100] 100] 100 | 160-1 °100 1) 100 1100) =1007) 00 
anf 200 | 200] 200] 200} 200) 200] 200] 200| 200| 200 
Sf 1.92 | 1.88 | 2.06 | 2.10~). 1.62 | 1.92 | 1.76 ) 1064.) 2.08 boy ie 
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It may be added that these random sampling frequencies were all 
derived from tables of random sampling numbers constructed as described 
in a previous paper 1). 

To conclude this paper let us briefly discuss a practical application. 
In tables VII and VIII we have represented two sets of ten samples 
each of 1000 units taken from the population specified in table I; the first 
set (table VII) was constructed by using the random sampling frequencies 
of set Nr 1 in table IV in ten different permutations, while the second 
set (table VIII) was obtained by using the ten sets of table IV each once. 
For each sample we have computed the first, second and third moments 
about the origin, and 72, and from these the average moments and the 
standard deviations have been calculated. 


TABLE VII. 


Ten random samples of the Poisson distribution with m = 2 constructed from one set of 
random sampling frequencies in ten different permutations, 


i ee 


Sample 10 sample: constructed from set Nr 1, table IV 
: Theory 
class Frequencies . 
0 127 | 146 | 147 | 148 | 155 | 146 | 145 | 128 | 138 | 129 | 135 
1 282 | 267 . 257 | Wl | 247 | 272 | 270 | 265 | 264 | 277 271 
2 285 | ray eof 290 +f, 209 1 265 253 | 271 |) 276.) 266 Zi 
3 157 | 178 | 178 | 170 | 178 | 190 | 183 | 177 | 176 | 179 181 
4 90 | 88 | 98 2B! 89 89 97 | 110 89 | 100 90 
5 42 | 35 | 40 31 48 29 36 30 38 30 36 
6 14 10 | 17 18 11 11 13 M7, 13 
iz 1 an es 3 1 3 2 6 16 
8 2 1 | oot — = 74 1 = ae 
1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 
m'; | 2.000} 1.966 2.032) 1.946) 2.002) 1.928) 2.996} 2.049} 2.006] 2.016 
wont O01 1) 5.96,) 6.29 | 5.91 | 6.09 | 5.53) 6.05.1 6,23 | 6.04 | 6.06 
Me i22.2 12933 23.6 22eo0 227s) Be O 2234 25.0) (22.0) 122.2 


i aeat 2528 | $22 119.03 |. 9.53 | 6226 | 2.52 | 6.20 | 0.98. | 3.09 


m,= 1.994; 8m, = 0.036 
m',= 6.007; $m', = 9-200 
m'; = 22.02 ; Sm, = 1.19 


In computing z2 the frequencies of the last three rows were pooled, 
which reduces the number of classes to seven and the number of degrees 
of freedom to six; the 95 % and 5% points of the y?-distribution are 
then 1.64 and 12.59 respectively. On comparing with tables VII and VIII 
we see that out of a total of 20 values of y2 two (0.98 and 1.21) lie below 
the 95% point and one (13.03) lies above the 5% point: quite an 
acceptable result. It will be noted that these three somewhat outlying 


1) H. C. HAMAKER, Proc. Kon. Ned. Akad. v. Wetensch., Amsterdam, 52, 145—150 
(1949). 
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values all occur in table VII, and if we restrict our attention to this table, 
we have an event with a chance 1/10 that has occurred 3 times in a series 
of 10 trials: the binomial probability of three or more such events is 
easily computed to be 0.07, which is still above the customary significance 
limit. 

TABLE VILL 


Ten random samples of the Poisson distribution with m=2 constructed from ten 
different sets of random sampling frequencies. 


Pe ee ee 
ca Neild 2 eas. lo 4s ee ee [gy > 1 a 
Theor 
Sample Frequencies sf 
class 
ESS foe i ae ee eee 
0 193 | 131 | 138-4 124 1-143.) 192 11319) 14f 146 aa 135 
i 276 | 249. | 253 | 288 | 265 | 285 | 257-| 267 | 268 | 295 9 2x8 
2 a67 “W270 | 278 | 2731 267 | 278 | 296) 263. 1-247 7 see art 
3 199 | 187 | 173 | 176 | 194 | 162 | 174 | 178 | 187 | 185 | 181 
4 77 | 88 | 98 |. 90 | 844 -91~f—00 1) 91 | 98.) 26 90 
5 40 C95 4 38 1 34) Sar "aa ae ae 36 
6 12M 9 shel ung E 10 | 14.) 14) 16 7 10 
Z 4 1 ff Cal eee 3 5 4 4 4 
8 1 1 1 Pees e 1 2 1 a 16 
9 a ee eee i = = 
1000 | 1000 | 1060 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 
m'; | 2.028] 2.074] 2.057] 1.987) 1.962] 1.974] 2.018] 2.021] 2.003] 1.985 
m', | 6.08 | 6.38 | 6.44 | 5.86 | 5.71 | 5.88 | 6.06 | 6.23 | 6.07 | 5.82 
inig 122.03 23:7 125.1 [211 199.8 1225 (22:3 123.6. e222 ae 
y2 | 5.71 | 5.86 | 4.45 | 2.36 | 4.59 | 3.23 | 4.61 | 2.57 | 6.19 | 4.43 
m', = 2.011; Sm, = 0.024 
m= 6.053; 5.9, 0220 
m’s = 22.24 : Sg = aT 


Moreover, the ten samples of table VII were derived from the set of 
random sampling frequencies Nr 1 in table IV which possesses a standard . 
deviation somewhat below the theoretical value; if anything, we should 
therefore expect the variability between the samples of table VII to be 
less than normal, though it is extremely unlikely that such a tendency 
should already be apparent in so small a series of samples. These 
arguments render it highly probable, however, that the comparatively 
large variations in y2 in table VII are accidental, and are not in any way 
connected with the particular method by which the samples were obtained. 

The average moments in tables VII and VIII lie close to their theoretical 
values, which are 2.00, 6.00, and 22.00 respectively. The standard deviation 
of the moments computed from each of these tables separately are in 
satisfactory agreement, and these standard deviations offer perhaps the 
best illustration of the usefulness of the methods described in this paper; 
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by constructing a set of large-size samples we may rapidly estimate the 
standard error of any statistic computed from them, while avoiding involved 
mathematical arguments which are all too often incomprehensible to the 
practical statistician; and the method also applies to cases not’ amenable 
to mathematical treatment. 

I should like to express my indebtedness to Mr. H. A. C. v. D. LINDEN 
for his careful assistance in carrying out this investigation. 


Eindhoven, 19 December 1947. 
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Crystallography. — Transformation of gnomograms and its application 
to the microchemical identification of crystals. I. By D. W. 
DyksTRA. (Memorandum of the Crystallographic Institute of the 
Rijks-Universiteit at Groningen. ) (Communicated by Prof. J. M. 
BIJVOET.) 


(Communicated at the meeting of February 26, 1949.) 


Summary. The method of construction, dealt with in this article, 
allows: 1° the transition to another plane of projection, crystallographic 
axes of reference and parametral plane remaining the same (rotation, 
Umwalzung), 2° the transition to other axes and parametral plane (trans- 
formation in a narrower sense). 

Here it is not the crystal, that will be rotated (as usual) into another 
position with respect to the plane of drawing, which does not alter its 
position, but the crystal remaining stationary the projecting lines are made 
to cut out another gnomogram in a new plane, which may have any position 
with respect to the original plane of drawing. This plane is revolved to 
coincide with the picture plane, and then the old gnomogram and the new 
one will be found to be related by the principle of central collineation. The 
form of the collineation is a very simple one, so the constructions can be 
effected with only few lines and practically only with a ruler (§ 6, 8, 9 
and 10). 

The most usual types of gnomogram are considered to be particular 
cases of one general gnomogram (§ 1, 2 and 3). 

As an application the identification of a crystal of silverdichromate, 
measured by microscope, will be described (§ 11b). 


§ 1. In the most usual gnomogram the zone lines are equidistant 
parallel lines because of the particular position of the plane of projection, 
the latter being perpendicular to one of the axes of reference. 

In the general form of the gnomogram, however, the plane of projection 
may have any position. 

In practice this may occur inter alia in the under mentioned cases: 

a. In measuring efflorescent or deliquescent crystals time may be 
insufficient to adjust the crystal on the two-circle goniometer. Measurements 
are taken with the crystal in an arbitrary position and this will result in a 
general gnomogram. 

b. With some crystals the planes are arranged in a small number of 
zones, intersecting in a face actually occurring on the crystal. In measuring 
these ‘‘face adjustment’ 1) will be profitable. When plotting the measure~ 


1) M.H. HEy, Min. Mag. 23, 560 (1932), 
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ments one will find a gnomogram with the plane of projection coinciding 
with a face of the crystal, and the latter being in general not perpendicular 
to a possible edge, the form of the gnomogram will be the general one. 

c. When reducing Laue photographs one often meets with gnomograms 
with arbitrary planes of projection 2). 


§ 2. When passing from the particular to the general gnomogram the 
normal pattern of fig. 1a changes into that of fig. 1b. 


This figure 1b is a complete quadrangle with as its angular points the 
“octahedral faces”: (111), (111), (111) and (111) (for the sake of 
clearness its sides are represented by full lines, its diagonals by dotted 
ones). The “cube faces” (100), (010) and (001) are the diagonal points; 
the four zone lines through each “cube point” form a harmonious pencil of 
lines, being two sides and two diagonals of the complete quadrangle. 

In order to construct in this general gnomogram the zone lines cor- 
responding to the equidistant parallel zone lines of the ordinary gnomo- 
gram represented in fig. 2a, the rays (100)—(001) and (100)—(111) in 


Fig.2a Fig. 2b 


2) See e.g. R. W. G. WYCKOFF, “The structure of Crystals” (1924) p. 139, fig. 108. 
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fig. 2b (the points (100), (010), (001) and (111) being given) are inter- 
sected with a line parallel to the ray (100)—(010), and then on this line 
pieces are laid off “equal to the piece cut off by the rays just mentioned 
(fig. 2b). The connection of the division points, found in this way, with 
the point (100) results in the pencil through (100), and then a correspond- 
ing construction gives the one through (010). 


§ 3. a. The first particular case of this general gnomogram mentioned 
here, is the ordinary gnomogram; the line (100)—(010) in this case is the 
line at infinity of the plane, and that converts the pencils into systems of 
equidistant parallel lines. 

b. Another particular case is that of the trigonal notation according 
to MILLER. Here the “dodecahedron zone” (110)—(011)—(101) 3) is the 
line at infinity (fig. 3). 


Ill 


/ 


Fig.3 


c. In the regular system the “triangular gnomograms” are identical 
with those according to MILLER’s notation; in the other systems their form 
is the general one *). 


§ 4. In the usual methods 5) of transforming gnomograms the crystal 
will be made to rotate the plane of drawing remaining stationary (Um- 
walzung), which results in the gnomogram thus constructed having the 
same gnomon circle as the original one. The new place of the pole is found 
by construction or by means of a “gnomonic net”, 


3) This is the zone [111]. 
*) R. L. PARKER, Schw. Min, Petr. Mitt, XVII, 475 (1939). 
P. TERPSTRA, ‘Kristallometrie”, p. 103 seqq. 


°) See eg. V. GOLDSCHMIDT, “Uber Projection und graphische Krystallberechnung” 
(1887) p. 68. 


H. E, BOEKE, “Die gnomonische Projection etc.” (1913) po. 26: 
F, E. WRIGHT, Am, Min. 17, 422 (1932) 


. 
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Let us take WRIGHT's method by way of example. (In § 7 WriGuT's 
construction will be shown to become formally identical with this method 
of ours, by only a slight alteration.) 

In fig. 4 the crystal is to be rotated, the result being the coincidence of 
the pole P’ with the centre of the gnomon circle. First thing to do is to 
construct the angle of rotation P’C,O; this is bisected by the line C,Rp, 
and in R>, the line i, is drawn perpendicular to P’O. At the same distance 
from O, but this time on the other side, the line ig //i; is drawn. These 
lines i; and ig are WRIGHT's “isometric lines”. 


Fig.4 


The rotation makes a point D’ on i, coincide with a point D on ig. 

Next the trace p’ of the plane P’ is constructed; now the point Rp will 
be the angle point of the line p’ (4 A’RpCr = 90°—14 ZOCuP”, 
ZA'CrR; is so too; so A’Cn = A’R,). By rotation p’ goes to infinity, 
making O the new angle point. This means that after the rotation a point 
E’ on p’ has come at the point at infinity of the line drawn through 
O|/R,£’. Thus the point at infinity of the line R»E’ will represent the 
point E corresponding to the point E’. 

The zone line / conjugate to a given zone line I’ through D’ and E’ is 
obtained by connecting D and E, this being the line through D parallel 


to R git 


§ 5. In our method the crystal itself is not rotated; the projecting rays 
from the centre O, however, are made to intersect not only the original 
plane of drawing but another one too, which may have any position with 
respect to the former, provided that its distance to the centre of projection 
— as that of the former plane — is equal to the radius of the gnomon- 
circle. Each pair of corresponding points lying on the same ray from the 
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centre of projection, the figures in both planes are in perspective with one 
another. 

In fig. 5 +’ represents the original plane of drawing in which lies the 
gnomogram to be transformed, and 7+ the new plane of drawing. Both 
planes are tangent to the sphere, having the centre of projection O for 
centre and the gnomon for radius, and their intersection is the line d. The 
plane through O perpendicular to d is the plane of drawing in our figure. 
O’ is the centre of the gnomon circle and P” is the pole of the plane z in 
the original gnomogram. The plane tp, through O parallel to t, intersects 
tv’ by the vanishing trace v’; the plane t’p, through O parallel to 7’, inter- 
sects t in the vanishing axis w. 


Fig. 9 


The distances of the parallel sides in the parallelogram v’Owd being 
equal v’Owd is a rhomb, which makes v’O = v’d and v’d — dw. 

Now tp is revolved round v’ into 7’, and t in the same direction round d, 
likewise into 7’, This revolution makes O. (afterwards called On) come on 
d and the vanishing axis w at wn, v’d being equal to dw,. The centre P of 
the gnomon circle in the new plane 1 comes at M, while O’d = dP — dM. 
Now the figure in the plane 7’ will look like fig. 9 without the point O 
drawn there. 

The revolution has made the old gnomogram in 7’ and the new one in z, 
formerly being in perspective to one another, now to be related by the 
principle of central collineation, the line d being the principal axis of the 
collineation, the lines v’ and wa, on either side of and at an equal distance 


from d, being the secondary axes and the point On on the line d the centre 
of collineation. 


The central collineation has the followin 


g properties, as will be proved 
in § 6: 
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a. Each line I’ of the original gnomogram intersects its conjugate I of 
the new one in a point on the principal axis d. 

b. Each line through On coincides with its conjugate. 

c. A set of parallel lines I’ changes into a pencil / through Q, the point 
of intersection of the secondary axis wz and the line through On parallel 
to I’ (fig. 6a). 


————-} ——} Fb 


#nQ My 
Fin 


d. A pencil ’ through the point P’ on the secondary axis v’ changes 
into a set of parallel lines J, parallel to the line connecting P’ and On 
(fig. 6b). 

e. Conjugate points A’ and A lie on one straight line through On. 


§ 6. In order to prove these properties fig. 5 has been copied stereo- 
metrically in fig. 7. The letters are the same as in figg. 5 and 6, in figg. 5 
and 7 the index n being added after revolution. 
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a. The rays from O, projecting the line 1, describe a plane through O 
and 1;; the intersection of this plane and 7’ is the line I’; conjugate to 1,. 
The intersection D of this plane and the line d is a point common to 
and 1’, and this point D remains stationary during the revolution; so the 
point of intersection of I’, and In too will lie on the principal axis d. 

b. The rays from O projecting the line / through On will lie in a plane 
through O and 1. The parallelogram O,P’OQ formed by the intersections 
of this plane with the planes 1, v’, tp and 1’ will be a rhomb. 

During the revolution P’On remains stationary, the parallelogram hinges 
in its corners, OQ remaining parallel to P’On. When finally O has come 
at On, Q coincides with Qn on P’O, produced, so OnQ (= 1) becomes 
OnQn (= In), In being I’ produced. 

c. The intersections of any set of planes through OQ with 7 will form 
a pencil of lines J, 1, ..., having Q for its vertex, and those with 7’ a set 
of parallel lines ’,/..... After the revolution the pencil will lie in the 
plane v’ with Qn for its vertex, OnQn being parallel to the set I’, I’, .... 

d. The intersections of any set of planes through OP’ with 7 will form 
a pencil of lines l’, I’, ..., having P’ for its vertex, and those with 7 a set 
of parallel lines J, 1, ..., forming, after revolution into 7’, a set of lines 
parallel to P’On. ; ; 

e. The line OnA’ is transformed into a conjugate line OnA, coinciding 
— as shown in section (6) — with O,A’; which is another way of saying 
that On, A’ and A will lie on one straight line. 


§ 7. In this passage will be shown, that only a slight alteration in 
WRIGHT’s construction will bring about formally the samé results as those 
obtained in § 5. ‘ 

In fig. 4 all the zone lines I’, meeting at the point E’ on p’, are shown 
to change into a set of lines 1, parallel to one another because they are all 
parallel to R,E’. 

In order to deduce from the unaccented zone lines the corresponding 
accented ones the rotated crystal is to be turned back through the same 
angle, which results in Q coming at O (fig. 8). q, the trace of the plane Q, 
being located, the construction is an analogous one. 

The zone line J meets ig at D and q at F; so I’ is drawn through D’ and 
F’, that is through D’ // OF* |/ RoF. 

This way of constructing results in a set of accented lines, parallel to 
one another, all of them being parallel to RaF, corresponding to a set of 
unaccented ones, meeting at the point F on q. In order to complete the 
conformity with the method based on the principle of central collineation, 
another modification is made in the drawing represented in fig. 4. The 
full lines in fig. 8 are the same as those in fig. 4 and the letters too are 
placed in the same way. In fig. 4 two gnomograms are drawn with the 
same gnomon circle, one accented, the other one unaccented. 

The unaccented gnomogram is, as a whole, shifted downwards across a 
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distance a = OR, = OR,,; the accented one is translated upwards across 
this distance a. After translation lines and points all have got an asterisk 
and the lines have been changed as follows: ————_— into --------- 
ts a a into -++-+--+---+-.. . Now the gnomon circles will be 
separated and i, and i, will coincide with one another. 


Fig.8 


The way of constructing conjugate lines is now identical with that 
described above in §5; ij i taking the place of the principal axis of the 
collineation (fig. 7), p’* that of the vanishing trace v’ and q* that of the 
vanishing axis wn. Rj (Rj) is the centre of the collineation On. 


§ 8 The transformation dealt with in § 5 can be performed easily, 
when the lines d, v’ and wn and the centre On have been located. 
In most cases these lines and the point On will have to be constructed. 


Three of them will be described here: 


a. A given gnomogram is to be transformed inte a gnomogram on a 
plane indicated by its pole P’ (fig. 9). With regard to an arbitrary plane 
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one has to solve this problem, when wishing to reproduce a face in 
its true form by means of orthogonal parallel projection onto the plane 
of drawing, e.g. for the construction of a model of some crystal or other 
in pasteboard. 

The zone line to the pole P’ (cf. fig. 5) is constructed; this is the 
secondary axis v’. The line bisecting < O’OP” intersects the original plane 
of projection at On; the line through On /|/ v’ will be the principal axis d; 
the other secondary axis wn is the line parallel to d and at an equal distance 
from d as v’. The centre M of the new gnomon circle is found by making 
O’O, = OnM, the radius not being changed. In the gnomogram now 
constructed the pole of face P will coincide with M. 

b. A transformation of a given gnomogram is required, resulting in a 
certain line v’ in the original one changing into the line at infinity of the 
new one, this being the case when a gnomogram in the general form is to 
be transformed into one with equidistant parallel zone lines. 

First P’, the pole to the zone line v’, is located and then the construction 
_ is completed as above (problem (a) ). 

If the zone line v’ should lie beside the drawing paper and so the pole 
close to the centre of the gnomon circle, to begin with, a transformation 
should be performed promoting an arbitrary but farther removed point 
to centre of the gnomon circle, which then should be followed by another 
transformation converting the line v’ of the second gnomogram into line 
at infinity. 

c. For the sake of completeness the case will be dealt with, that the 
line at infinity is to be transformed into a given line wa (figg. 5 and 10). 
In fig. 10 only the line wy and the point O’ are given. 


Fig 10 


Now (fig. 5) OnwnwO will be a parallelogram, Ow and Onwn being 
parallel and equal to one another. The point of intersection S of the 
diagonals bisects the line Own. So from S (fig. 10) the tangent to the 
circle is drawn and the point d results. The point M and the line v’ are 
constructed as above (problem (a)). 


§ 9. In order to change to other axes of reference and another para- 
metral plane the complete quadrangle through the given poles is con- 
structed, as indicated in the first part of this article, and this is transformed 
into an ordinary gnomogram as in § 8b. 
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§ 10. The construction is performed as follows. The line 1 conjugate 
to the given line /’ is to be constructed. 

a. The line Y’ (fig. 11) can be regarded as one out of the pencil through 
B (cf. fig. 66), which will change into a set of parallel lines. The ray BOn 
out of this pencil coincides with its conjugate and the line / conjugate to 


can be drawn through A // BOn. 


Fig. 11 Fig.12 


b. The line l (fig. 12) can be regarded as one of a set of parallel lines 
(cf. fig. 62), which will change into a pencil with its vertex on wn. The 
ray OzC out of this set coincides with its conjugate, which makes C the 
required vertex. The line / conjugate to l’ is drawn through A and C. 

This construction of the point P conjugate to P’ (fig. 13) turns out to be 
guite the same as that of the image formed by a negative lens, On perform- 
ing the function of the optic centre and C that of the focus. Accordingly 
the relation between the positions of conjugate points complies with 


45 1 eee re 
the well known formula for lenses, giving here 0,P~ 0,P°=0,D’ 


Physics, — Photons in extensive Cosmic-Ray-Showers. By J. CLAY. 
(Natuurkundig Laboratorium Amsterdam.) 


(Communicated at the meeting of March 26, 1949.) 


Summary. 


The results are given of 6 and 8-fold coincidence measurements of 
extensive showers with layers of lead above and between the counters, by 
which the absorption of the electrons and the production of electrons in 
the showers could be estimated. It is clear that there are a great deal of _ 
energetic photons in these showers and the amount increases with the 
density of the shower particles. 

By measuring 4- and 6-fold coincidences for surfaces of relative sizes 
from 1 to 6 the value of the coefficient S in the frequency-density relation 
I =I,-N-~S could be found and we come to the conclusion that S is not 
constant but varies from S = 3,0 for very small densities to S = 1,5 for 
high densities, which is in agreement with the frequency-density relation 
for bursts. 


Some time ago in the extensive atmospheric showers we found an 
indication of a great number of photons. We were carrying out an 
experiment for measuring on one hand the absorption of the electrons in 
these showers and on the other hand the percentage of the particles therein 
which penetrate at least 10 cm Pb. To this purpose we were using 4 trains 
of two counter-sets each, which were placed in a quadrangle, each train 
arranged so, that the lower box could be protected by 10 cm Pb. We were 
anxious to measure afresh the number of penetrating particles in relation 
to the total number of particles, as we had found in two separate sets of 
measurements this relation to be 1 to 10. At the Pasadena congress of 
1948 Cocconl protested that he had found a relation which was totally 
different from the one found by us. Therefore we repeated our experiments 
with two different counter arrangements in two different places, in our 
laboratory and in a wooden shed with a thin wooden roof where we had 
already made several measurements before. In both places we again found 
the relation 1 to 10. Vide table. 

We have already stressed that when we observe two particles there are 
two possibilities. It may be two particles belonging to one extensive shower 
and it may be two coherent mesons produced in one and the same process. 
The variation with distance found while measuring 4-fold coincidences 
caused by two particles in two vertical counter-trains at not too large a 
distance from each other appears in either case. We thought however that 
we would be able to distinguish between the two cases when measuring 
2 or 3 particles. And indeed when we suppose the variation of frequency 
of the coincidences with distance to be given by N = Noe-**, we find in 
the case regarding 2 particles at distances of 120 cm and 60 cm a = 0,47 
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and a = 0,37 respectively and in the case of 3 particles at the distances 
mentioned a — 0,62 and a = 0,64 cm-! resp. In the frequency-distance 
curve for the coincidences we find a very distinct difference between the 
variation for slightly larger distances (up till 5 m) for two penetrating 
particles of extensive showers and the variation of the frequency of coin- 
cidence of two coherent mesons which here are measured together 1). 
Above 5 m distance the influence of the coherent mesons is negligible. 
JANoss! 2) is of opinion that also the showers are to be divided in two 
groups, viz. extensive showers and locally more limited ones. 

Two years ago already we came to the conclusion that the extensive 
showers were very rich in photons. We therefore tried to find the pro- 
portion of the photons to the electrons in the following way. We used 
four trains of two sets of counters each, these sets separated by a space 
of 10 cm designed to be filled with lead screens of different thickness. In 
case there is no lead between the counters, the particles, mesons and 
electrons, give coincidences and we can measure at the same time 4 and 
6-fold (2 and 3 particles) and 6 and 8-fold (3 and 4 particles). When 
we apply lead between the counters, the absorption of the soft part of the 
radiation, the electrons, can be measured. From the absorption curve the 
energy distribution may be derived. With 10 cm of lead the electrons are 
practically all absorbed, and we find that the incident radiation contained 
a number of mesons equalling about 1/10 of the number of electrons. This 
ratio between the two kinds of particles is totally different from their 
ratio in the radiation of the atmosphere outside the showers. 

When next we put a sheet of lead above the upper set of counters we 
find an increase. This is on account of the photons in the shower which 
produce in the top layer of lead Compton electrons or electron-pairs. These 
secondaries are able to operate the counters under the layer. The absorption 
of the primary electrons is known from measurements with lead between 
and above the counters. To determine the number of photons we proceed 
in the following way. 

We know that when the active area of m counters, each of an area f, is 
called F and WN is the density per m?, the probability of an n-fold coin- 
cidence is (1—e-”f")" or (1—e-F")". If there are P photons which have 
an opportunity of producing in the top layer an electron able to operate 
the counters under the layer, with an efficiency 4, then the probability of 
an n-fold coincidence in which n, counter-sets are discharged by either 
primary electrons or photon-produced electrons and n—n, by primary 
electrons only, is 

{1 — e~FIN+PH)} 1, (1 — eFN)n-m; 0 ny = 7. 
For each value n, this expression has a special weight gn, and the total 
probability for an n-fold coincidence with lead on top of the counter-sets 
is given by 


Soe Lette (1 eh yom, 


1 
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The experimental results however gave us an indication that this 


lin power of the 
n 
ratio a of n-fold coincidences with and without a layer of lead on top is 
nearly a constant, as appears from the table hereunder: 


formula might be simplified in the following way. The 
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Consequently we may conclude that in the case with lead above the 
counter-trains, the electron density N increases by AP, giving (l—e~F.W**))” 
for the probability of an n-fold coincidence. 

The ratio of n-fold coincidences with and without lead may therefore 
be expressed as 

(1 —e-FIN+iP))n 
. (1 —e-FN)n 

in agreement with the above-mentioned experimental results. 

If F is known, we can compute N from the measurements without lead, 


=="a 


so that e-FN is known too. Substituting then e-F*? =x and eFN=y 
enables us to calculate 1 P in the following way: 


= i 
jay ja=1—gx 
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From earlier experiments we found for counters 0,03<1< 0,06. If we 
now introduce 4 = 0,06 the value may turn out to be too small in this case. 
But we are experimenting at the time to find a more appropriate value. 

The measurements treated in the present paper were made in two series. 
One of them was made in our laboratory under a roof with a concrete 
ceiling of 20 cm. Four vertical trains of two sets of counters each were 
used, so that it was possible to measure 2, 3 and 4 particles mainly in a 
vertical direction. 

Next we experimented with the counter-sets placed side by side and 
6 and 8 coincident particles were measured. These measurements were 
repeated with 1 cm Pb above the counter trains. With this arrangement 
we measured the particles and photons from every direction. 
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Each set contained 3 counters of 30 cm length and 3,5 cm diametre so 
that the active area of a set was 315 cm2. N was found 35 p m2, Using | 
this value, we find AP — 16, giving P = 320 for 4 = 0.06. 

In a second series of measurements 4 and 6-fold coincidences were 
counted with sets of counters of 70 cm length and 3 cm diametre. 1, 2, 3, 
4, 5 and 6 of these counters in turn were used in parallel. As we will see, 
the minimum density of the showers which we are able to measure depends 
on the active surface and is smaller when this area is larger, as we 
had already observed in 1941. In several cases the photons were measured 
too, with lead (1, 3, 5 and 6 counters in parallel). In all these cases the 


frequency of coincidence was measured moreover with a cover of 10 cm Pb 
to the counters. 


In one case we have determined the photon-density of the showers 
using very small surfaces and in doing so we have found a mean density 
of 150 times the density of the total of the particles of the shower. The 
density of the photons in this case was 4000 per m2 for = 0,06. 

For the first mentioned experiments the results are given in the graph 1, 
which shows that the increase by the lead layer is considerable. 

Likewise a very interesting result is that, when 1 cm Pb is placed between 
the upper and lower sets of counters, there is still an increase in the number 
of coincidences in comparison to that in counters without lead. How is this 
possible? This is a consequence of the coherence between electrons and 
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Fig. 1. 3 and 4-fold coincidences open and with different layers of Pb on top (drawn 

curves) and with Pb between the counter-sets (dotted curves). For comparison the 

absorption curve is given for the common soft electron component (1eh)e 
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photons. Due to the absorption of the electrons by the lead sheet between 
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Fig. 2 and 3. The variation of 6 (3p) and 8 (4p) fold coincidences with barometric 
variations. 


the counters a great number of the electrons cannot operate the counters 
under that sheet and this will decrease the number of coincidences. But 
there are photons in coherence with the electrons. And the presence of the 
lead sheet causes thesesphotons to produce secondary electrons. Therefore 
a coincidence may occur caused by a primary electron traversing the upper 
set and by a secondary electron produced by a coherent photon traversing 
the lower set of counters. It may be possible to calculate the effect from 
the respective numbers of electrons and photons but it is difficult to 
estimate the energy spectrum of the electrons which is certainly different 
from that of the electron of the ordinary soft electron component in the 
radiation which is given also in the graph for comparison. 

It is evident that one cannot measure the actual absorption of the pri- 
mary electrons directly by putting lead over the counters nor by putting 
the lead between the counters. 

Another interesting phenomenon is discernible in this series of obser- 
vations, but this complicates the issues considerably. It became manifest 
that barometric fluctuations caused great variations in the frequency of 
the observed showers. This influence is given in the graphs 2 and 3 and 
the coefficient of the variation in frequency is 13 % for 1 cm Hg. Conse- 
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quently we can compare the results of different situations only for identical 
barometric pressure. This big coefficient is further quite in agreement with 
the increase of the frequency of showers with increasing altitude observed 
by HILBerRY in the lower part of the atmosphere and it indicates that the 
production of the showers may lie not very high in the atmosphere. 

We also made a prolonged series of measurements in another building 
with sets arranged in a hexagon, observing 4 and 6-fold coincidences as 
mentioned above. Every set contained six counters and measurements were 
first made with one counter functioning in each box, then two, three ... up 
till six. From counting 4 and 6-fold coincidences in this manner we could 
compute the frequency of showers and their mean density. By putting 
0,5 cm Pb over the counters the photon density could be estimated and by 
placing 10 cm Pb over these boxes the penetrating part was determined. 

We know that with the increase of the active surface the ratio of the 
showers of low density to those of high density increases. It is even a well- 
known fact that with very small surfaces only the densest showers are 
found. There must be a relation between the frequency and the density of 
the showers and in analogy with what was found in the case of the meson 
and the electron spectrum, the most plausible assumption was to suppose 
that the frequency would be proportional to a power of the density. 
I(N)dN = I,N-S dN, We may test by experiments whether the value 
of S is of the same order as the exponent for the energy for the mesons. 
This has also been done by Cocconi, LOVERDO and Tonaiorai 3). 

If N is the density and f the active surface of the counters, the chance 
that a counter system of m counters will be hit is (1-—-e—™f%)" and 
the total probability of the coincidence for an area F = mf will be 


I(F)=Ip [N -S({—e-FN)n dN, 


Number ; 
of Absorber ie ae M 
in cm 
counters 


250 (4000) 


1—10cm Pb 64 
1 0 230 2,0 1,6 0,8 


lcm Pb 230 3,8 a4 26 (410) 
10cm Pb 230 
2 0 460 4,6 3,1 0,68 
3 0 690 9,7 6,3 0,68 
0,5 cm Pb 690 12,1 7,9 
10cm Pb 690 0,6 0,4 5 
4 0 920 13,3 7,9 0,60 
5 0 1150 19,1 9,7 0,51 
6 0 1400 26,0 14,1 0,55 
0,5cm Pb 1400 30,1 17,0 1,5 (24) 
10cm Pb 1400 3,6 1,7 4 


30 
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¥ 
After the substitution of FN = x this becomes ~ 
I(F) =I, FS“ f x-S (1 —e)" dex. 
0 


What we next want to find is the value of S. There are two ways: by 
varying F and by measuring for different values of n. The first has been 
done by Coccon!, LoverDo and TONGIORGI and gives mathematically the 
simplest result, as the integrals for two values of F are the same and 


He) _ (Bi) 


1D aa) 

gives directly the value of S which is constant for different N's, but it is 
more correct to derive the relation from different areas which have different 
observed mean densities. Therefore it is preferable to find S for different 
n’s and use the same surface per set. We have made measurements in both 
ways, as we have observations of the relation of 4 and 6-fold coincidences 
with sets of 1, 2, 3, 4, 5 and 6 counters in every set. 

For the relation of sixfold and fourfold coincidences the values which 
we found were: 


S Theoretical Experimental N (density) 
1,5 0,89 : 0,84 100 
2,0 0,78 0,68 19 
2D 0,63 0,60 16 
3,0 0,58 0,55 9 


The integrals 1(F) were calculated for a series of S and n values by 
Mr H. F. JoNGEN tio whom we are much indebted. 

With these results it is obvious that we may not take the value of S 
constant for all densities, but that with the increase of density the value S 
decreases, and therefore our assumption about S in J(F) is not perfectly 
correct, but this was a first attempt to solve the problem. 

These different values of S are in agreement with our results for the 
track density of bursts1), which is for the low density of the order of 
3,0 and for high densities 1,6. Also the value of proportion of two 
coincidences of different numbers n is different. 

For one counter the value 6/4f 5/4f 4/3f 3/27 

theoretical OFS Ot. O72 Boss 

experimental 0,8 0,6 0,365 701 

This agreement is not too good, but the tendency in the calculations and 
the observations is the same. 


Finally I wish to thank Mr G. KLEIN for his help with part of the 
experiments, 


Amsterdam, Febr. 1949, 
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Biochemistry, — On stearate systems containing methyl-hexylcarbinol with 
viscous and elastic properties comparable to elastic viscous oleate 
systems containing KCl. By H. J. VAN DEN BerG and L. J. DE HEER. 
(Communicated by Prof. H. G. BUNGENBERG DE JONG.) 


(Communicated at the meeting of March 26, 1949.) 


1. Introduction. 


The results obtained in this laboratory 1) with elastic viscous oleate 
systems containing KCl opened up the question if any systems with 
analogous properties also exist with other soaps. A priori one would be 
inclined to answer in the affirmative, as at sufficienlty high temperature 
palmitates and stearates also form elastic systems at KCl concentrations 
which are to low for coacervation. 

Difficulties to be foreseen in the measuring technique at such high 
temperatures withheld us for some time to start with such an investiga- 
tion. 

The elastic systems obtained at high temperatures cannot be inves- 
tigated after having cooled down to room temperature as they are no 
longer stable then (the soaps separate in an insoluble form). Prof. BUN- 
GENBERG DE JONG informed us that the temperature at with the elastic 
stearate system of the composition 1.2 % K-stearate in 0.2 NKOH 
becomes unstable, can considerably be lowered by adding the proper 
amount of a suitable organic substance, e.g. methyl-hexylcarbinol. It was 
therefore decided to perform some measurements on the above stearate 
systems for which we used “K-stearat, doppelt gereinigt, RIEDEL DE 
HAEN”, as this preparation makes it possible to get stable and clear 
systems at even 15°. 


2. Preparation of the elastic system. 

When the system (1.2 gr K-stearate p. 100 cc in 0.2 N KOH) has 
cooled down to room temperature, it is a thick semiliquid somewhat pasty 
system, of a white silky appearance. If at 15° one adds increasing amounts 
of methyl-hexylcarbinol and homogenizes it by vigorous shaking, it is 
found that gradually the white silky appearance becomes less and the 
system changes into a clear markedly elastic system, which on further 
addition of the alcohol grows turbid and finally looses its elastic 
properties completely. It was found in preliminary experiments (rotational 
oscillation in a spherical vessel with a radius of 3.58 cm) that as a 
function of the added amount of methyl-hexylcarbinol, the values 


1) Cf. H. G, BUNGENBERG DE JONG and H. J. VAN DEN BERG, these Proceedings 
part I, 51, 1197 (1948), parts II, III, IV, V, 52, 15, 99, 363, 377 (1949). 
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y (=1/T), 1/A and n increase to a certain maximum and next decrease. 
The maxima, which nearly coincide, are situated in the region of the 
completely clear systems. For the experiments in the following section 
we made a large quantity of elastic system corresponding to the position 
of the above maxima, for which an addition of 0.42 cc methyl-hexyl- 
carbinol was needed per 100 cc. 1.2 % K-stearate in 0.2 N KOH. 


3. Period and logarithmic decrement as function of the radius of the 
sphere at 15.7°. 


The methods used and the calculations of the elastic constants were 
the same as described by H. G. BUNGENBERG DE JONG and H. J. VAN DEN 
BERG in their publication of their investigations on elastic viscous oleate 
systems (see note 1, Part I, II and III). 

The results obtained with the rotational oscillation at 15.7° are 
collected in Table I and depicted in fig. 1. 


Fig; 1. 


We see from fig. 1 lower that the period is proportional to the radius R 
of the sphere. Compared with oleate systems this is quite the same result. 
In oleate systems we meet two different cases as regards the dependence 
of A on R. In 1.2% systems it was found to be proportional to 
R, in 0.6 % systems, however, it was found to be independent of R. We 
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TABLE I. 
Measurement with a 1.2% stearate system (containing 0.2 N KOH and 0.42 cc 


methyl-hexylcarbinol per 100cc) at 15.7°C (rotational oscillation), 


10x A=(T/2-A) 


(sec) (sec) 


G 
(dynes/cm?) 


b,/b3 


9.24 0.168 + 0.005 | 4 

7.46 | 93.8 | 6.54 | 1.124] 0.116 + 0.002 | 5. 64.2 »=GXl= 
mean mean 

5.65 | 86.2 | 5.02 | 1.096] 0.092 + 0.002 | 5.47) 5 4 | 62.5) 639 | 5:12 63.9= 

5.01 | 83.5 | 4.40 | 1.089] 0.085 ++ 0.003 | 5.16\ ~~ | 63.9 327 poises 

4.12 | 74.2 | 3.60 | 1.083] 0.080 + 0.003 | 4.51 64.6 


see from fig. 1 upper that in our 1.2 % stearate system too, A is proportional 
to the radius of the sphere. This means that here too the damping of the 
elastic oscillations must be ascribed to relaxation of the elastic stresses 
with a constant time of relaxation j 2). These values of 4 (calculated from 
T/2.A) are given in column 6 of Table I and fluctuate considerably around 
the mean value 5.12. : 

This is most probably connected with the relatively low damping (far 
lower values of A as were ever found in oleate systems). The lower the 
damping ratio b,/b3; the more experimental errors exert their influence on A 
(which is In 6,/b3). Such great fluctuations are not found in the calculated 
G values (Table I, column 7), which brings out that the b,/b3 measurements 
are far more difficult than the 10 X 7/2 measurements. 


Fig. 2, 


We therefore decided also to perform elastic measurements on a lower 
concentrated system, which we obtained by diluting three volumes of the 
above 1,2 % stearate system with one volume 0.2 N KOH. As it is highly 
probable that the methyl-hexylcarbinol is practically wholly bound to the 
stearate, an extra addition of it was deemed to be unnecessary. 


2) J. M. BURGERS, Proc. Kon, Ned. Akad. v. Wetensch., Amstedram, 51, 1211 1948). 
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TABLE II. 


Measurements with a 0.9% stearate system (containing 0.2 N KOH and 0.32 cc 
methyl-hexylcarbinol per 100 cc) at 15.7° C (rotational oscillation). 


T ( 7 G 
10*K= aw | 
R(cm)| 7 x 2| by/b2 A 2A (dynes/cm?) 
(sec) (sec) 

0.2232 | 4.11 32.6 n=Gxris 
0.1553 }-4,42{ mean | 32.44 RRR 8 age ase 
0.1389 | 4,364 4:30.) 43.7 \ SPro 4 0 oe ee 
0.1138 | 4.317 - 34.6 


The results obtained with this 0.9 % stearate system are collected in 
Table II and depicted in fig. 2, from which we once more see that T~R 
and A ~R, the experimental A points now deviating much less from the 
drawn straight line than in fig. 1 upper. Therefore the calculated 4 values 
(column 6) also show far smaller fluctuations percentually. 

That is why we feel justified in assuming that the irregularities in the 
A and / values met with in the 1.2 % stearate system, were caused by 
the difficulties in the measurement of the damping ratio. 

In Table I we have also calculated the value G X 4, which we will 
need in the next section. 


4. Viscous behaviour of the 1.2 % stearate system at 15.7°. 


For this investigation we used the same method as described by H. G. 
BUNGENBERG DE JONG, H. J. VAN DEN BERG and L. J. DE HEER (see note 1, 
Part. V). It may suffice to give here only the (logarithmic) viscosity- 
shearing stress diagram obtained. See fig 3, in which we notice the same 
characteristic shape as we met in the markedly elastic oleate system, 
though due to the very high viscosity at low shearing stresses it was not 
possible to make sufficiently accurate measurements to locate the mp level. 
The bending off at the left upper end of the curve, however, suggests 
that it must lie approximately as indicated by the dotted line. (log 7 
= 2.08, that is 79 = 120 poises). At the other end of the curve we have 
not yet reached the 7, level, which we may expect to lie at an 7 value 
in the order of 10 centipoises or lower. 

The total fall in viscosity from 7 to 7, by increasing the shearing 
stress is therefore very great, viz. at least to a thousandfold lower value. 

The most characteristic feature of the viscosity curve in fig. 3 is, 
however, the occurrence of an intermediate level. The total fall in visco- 
sity by increasing the shearing stress takes place in two separate steps, 
the first from + 120 poises down to +4 poises, the second from + 4 
poises down to 0.1 poise. (or lower). 

Such viscosity curves with three levels are characteristic of the 
markedly elastic oleate systems. 

Finally we may discuss a further point of resemblance between stearate 
and oleate systems. It was found for the latter that the viscosity coefficient 
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calculated from the elastic measurements (using MAXweELv’s formula 
»=G XA) is of the same order as the experimentally determined 1 
level, though 7 was found to be always smaller than G X 4 3), 

We have drawn in figure 3 the level corresponding with the value 
G X 4, which according to table I was found to be 327 poises (log y is 
therefore 2.51). Assuming the mo level to lie at approximately 120 
poises (log 72,08), as drawn in figure 3, 7) would be approximately 
2,7 times smaller than G X 2, 

For oleate systems the disagreement between 7 and G X / has been 
found to increase with the absolute value of 7 and from this one would 
have expected for an % value of 120 poises a disagreement which comes 
near to the above mentioned value of 2,7 times smaller. (i.e. %o/G X A 
e037). 


+3} bg 7 


+2 


5. Final remarks. 

The investigations described above show that the characteristics of the 
elastic and viscous properties of the soap system used are the same as 
those of oleate systems containing KCl]. They do not claim to give accu- 
rate figures for a pure stearate, however. In spite of the indication 
“doppelt gereinigt’”, the preparation of K-stearate used, appeared in 
further experiments to be anything but pure. BUNGENBERG DE JONG and 
pE HEER, in comparing a number of K-stearates and stearic acids of 

different origin, have found very great differences at to the damping under 


3) Cf. paper quoted in note 1, part V. 

4) P = shearing stress at the wall of the capillary in dynes/cem?; V = mean rate of 
flow in cm/sec. In.figure 3 the ordinate does not indicate log V, but log P—log Vv, 
which stands for the meaning of the logarithm of a viscosity coefficient 7. 
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comparable conditions. Really pure stearate systems (using a prepara- 
tion of Stearic acid “Kahlbaum” with m. p. 69°) cannot be transformed 
into clear elastic systems with methyl-hexylcarbinol at 15°, a0 sEt 
20°, but it does succeed at 40°. At 40° however, RIEDEL DE HAEN’s 
preparation, which we used in this investigation, was no longer measur- 
able because of the now very great damping (n, the number of turning- 
points being only 3) though it had a small damping at 15° (n having 
now values in the order of 60). 


Summary. 


1. It has appeared to be possible to make elastic-viscous systems, with 
analogous properties, with other soaps than oleates. One must slightly 
alter the methods of preparation though, in order to obtain systems which 
are measurable at a conveniently low temperature. 

2. Analogous to the results found for elastic measurements of oleate 
systems of more than 1,1 % Na-oleate, a direct proportionality was found 
between A and R, and T and R. 

3. The viscous behaviour was the same as has been found for oleate 
systems viz. a breakdown of viscosity with increasing shearing stress in 
a two-step process, which seems to be characteristic for these elastic soap 
systems. 

4, There exists a disagreement between the viscosity 79, measured at 
low shearing stress (at the wall of the capillary) and the viscosity predicted 
by MAxwWELL's formula 7 = G X A, which disagreement is found to be of 
the same order of magnitude as could be expected from the behaviour of 
the oleate systems. ; 
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